Binary_c parameters

The following chapter contains all the parameters that the current version of binary_c can handle, along with their descriptions and other properties.

This information was obtained by the following binary_c build:

binary_c git branch: releases/2.2.4 binary_c git revision: 6899:20230528:385edf07a Built on: May 28 2023 13:29:03

Section: stars

stars

Option

Description

metallicity

Description:

This sets the metallicity of the stars, i.e. the amount (by mass) of matter which is not hydrogen or helium. If you are using the BSE algorithm, this must be 1e-4 <= metallicity <= 0.03. See also nucsyn_metallicity and effective_metallicity.

Parameter input type:

Float

Default value:

0.02

effective_metallicity

Description:

This sets effective metallicity of stars as used in routines like the Schneider wind loss. If not set, or set to DEFAULT_TO_METALLICITY (==-1, the default), this is just the same as metallicity. The main difference between effective_metallicity and metallicity is the range of validity: 0 <= effective_metallicity <= 1, while metallicity’s range of validity is limited by the stellar evolution algorithm (so, for BSE, is 1e-4 <= metallicity <= 0.03).

Parameter input type:

Float

Default value:

-1

Macros:

[‘DEFAULT_TO_METALLICITY = -1’, ‘Random variation : double between 0 and 0.03 ‘]

M_1

Description:

The initial mass of star one (in solar units, internally this is star index 0).

Parameter input type:

Float

Default value:

0

M_2

Description:

The initial mass of star two (in solar units, internally this is star index 1).

Parameter input type:

Float

Default value:

0

M_3

Description:

The initial mass of star three (in solar units, internally this is star index 2).

Parameter input type:

Float

Default value:

0

M_4

Description:

The initial mass of star four (in solar units, internally this is star index 3).

Parameter input type:

Float

Default value:

0

vrot1

Description:

The initial equatorial rotational velocity of star one (in km/s, internally this is star index 0). If 0.0, the Hurley et al 2000/2002 prescription is used to set the main-sequence velocity, so for a truly non-rotating star, set to something small (e.g. 0.001). See also vrot2,3,4.

Parameter input type:

Float

Default value:

0

Macros:

[‘VROT_BSE = 0’, ‘VROT_BREAKUP = -1’, ‘VROT_SYNC = -2’, ‘VROT_NON_ROTATING = -3’, ‘Random variation : double between 0 and 300 ‘]

vrot2

Description:

The initial equatorial rotational velocity of star two (in km/s, internally this is star index 1). If 0.0, the Hurley et al 2000/2002 prescription is used to set the main-sequence velocity, so for a truly non-rotating star, set to something small (e.g. 0.001). See also vrot1,3,4.

Parameter input type:

Float

Default value:

0

Macros:

[‘VROT_BSE = 0’, ‘VROT_BREAKUP = -1’, ‘VROT_SYNC = -2’, ‘VROT_NON_ROTATING = -3’, ‘Random variation : double between 0 and 300 ‘]

vrot3

Description:

The initial equatorial rotational velocity of star three (in km/s, internally this is star index 2). If 0.0, the Hurley et al 2000/2002 prescription is used to set the main-sequence velocity, so for a truly non-rotating star, set to something small (e.g. 0.001). See also vrot1,2,4.

Parameter input type:

Float

Default value:

0

Macros:

[‘VROT_BSE = 0’, ‘VROT_BREAKUP = -1’, ‘VROT_SYNC = -2’, ‘VROT_NON_ROTATING = -3’, ‘Random variation : double between 0 and 300 ‘]

vrot4

Description:

The initial equatorial rotational velocity of star four (in km/s, internally this is star index 3). If 0.0, the Hurley et al 2000/2002 prescription is used to set the main-sequence velocity, so for a truly non-rotating star, set to something small (e.g. 0.001). See also vrot1,2,3.

Parameter input type:

Float

Default value:

0

Macros:

[‘VROT_BSE = 0’, ‘VROT_BREAKUP = -1’, ‘VROT_SYNC = -2’, ‘VROT_NON_ROTATING = -3’, ‘Random variation : double between 0 and 300 ‘]

fKerr1

Description:

Sets the angular momentum of star 1 to this fraction of the Kerr maximal angular momentum. Must be 0<=fKerr1<=1. See also fKerr2,3,4.

Parameter input type:

Float

Default value:

0

fKerr2

Description:

Sets the angular momentum of star 2 to this fraction of the Kerr maximal angular momentum. Must be 0<=fKerr2<=1. See also fKerr1,3,4.

Parameter input type:

Float

Default value:

0

fKerr3

Description:

Sets the angular momentum of star 3 to this fraction of the Kerr maximal angular momentum. Must be 0<=fKerr3<=1. See also fKerr1,2,4.

Parameter input type:

Float

Default value:

0

fKerr4

Description:

Sets the angular momentum of star 4 to this fraction of the Kerr maximal angular momentum. Must be 0<=fKerr4<=1.See also fKerr1,2,3.

Parameter input type:

Float

Default value:

0

vrot_multiplier1

Description:

Multiplier for vrot1 or Prot1. Defaults to 1.0. Ignored if vrot1 > 0.0. See also vrot_multiplier2,3,4.

Parameter input type:

Float

Default value:

1

vrot_multiplier2

Description:

Multiplier for vrot2 or Prot2. Defaults to 1.0. Ignored if vrot2 > 0.0. See also vrot_multiplier1,3,4.

Parameter input type:

Float

Default value:

1

vrot_multiplier3

Description:

Multiplier for vrot3 or Prot3. Defaults to 1.0. Ignored if vrot3 > 0.0. See also vrot_multiplier1,2,4.

Parameter input type:

Float

Default value:

1

vrot_multiplier4

Description:

Multiplier for vrot4 or Prot4. Defaults to 1.0. Ignored if vrot4 > 0.0. See also vrot_multiplier1,2,3.

Parameter input type:

Float

Default value:

1

Prot1

Description:

The initial equatorial rotational period of star one (in km/s, internally this is star index 0). See also Prot2,3,4.

Parameter input type:

Float

Default value:

0

Prot2

Description:

The initial equatorial rotational period of star two (in km/s, internally this is star index 1). See also Prot1,3,4.

Parameter input type:

Float

Default value:

0

Prot3

Description:

The initial equatorial rotational period of star three (in days, internally this is star index 2). See also Prot1,2,4.

Parameter input type:

Float

Default value:

0

Prot4

Description:

The initial equatorial rotational period of star four (in days, internally this is star index 3). See also Prot1,2,3.

Parameter input type:

Float

Default value:

0

inclination1

Description:

The initial inclination of star one (in degrees).

Parameter input type:

Float

Default value:

0

inclination2

Description:

The initial inclination of star two (in degrees).

Parameter input type:

Float

Default value:

0

inclination3

Description:

The initial inclination of star three (in degrees).

Parameter input type:

Float

Default value:

0

inclination4

Description:

The initial inclination of star four (in degrees).

Parameter input type:

Float

Default value:

0

B_1

Description:

The initial magnetic field of star one (in Gauss, internally this is star index 0).

Parameter input type:

Float

Default value:

0

B_2

Description:

The initial magnetic field of star two (in Gauss, internally this is star index 1).

Parameter input type:

Float

Default value:

0

B_3

Description:

The initial magnetic field of star three (in Gauss, internally this is star index 2).

Parameter input type:

Float

Default value:

0

B_4

Description:

The initial magnetic field of star four (in Gauss, internally this is star index 3).

Parameter input type:

Float

Default value:

0

B_inclination1

Description:

The initial inclination of the magnetic field of star one (in degrees).

Parameter input type:

Float

Default value:

0

B_inclination2

Description:

The initial inclination of the magnetic field of star two (in degrees).

Parameter input type:

Float

Default value:

0

B_inclination3

Description:

The initial inclination of the magnetic field of star three (in degrees).

Parameter input type:

Float

Default value:

0

B_inclination4

Description:

The initial inclination of the magnetic field of star four (in degrees).

Parameter input type:

Float

Default value:

0

stellar_type_1

Description:

Set the stellar type of star 1 (internal index 0), usually MAIN_SEQUENCE (main sequence). Note that setting the stellar type only works for stars with both age=0 and core_mass=0, i.e. main sequence (hydrogen or helium), white dwarfs, black holes and neutrn stars.

Parameter input type:

Integer

Default value:

0

Macros:

[‘LOW_MASS_MS = 0’, ‘MS = 1’, ‘HG = 2’, ‘GIANT_BRANCH = 3’, ‘CHeB = 4’, ‘EAGB = 5’, ‘TPAGB = 6’, ‘HeMS = 7’, ‘HeHG = 8’, ‘HeGB = 9’, ‘HeWD = 10’, ‘COWD = 11’, ‘ONeWD = 12’, ‘NS = 13’, ‘BH = 14’, ‘MASSLESS_REMNANT = 15’, ‘LOW_MASS_MAIN_SEQUENCE = 0’, ‘MAIN_SEQUENCE = 1’, ‘HERTZSPRUNG_GAP = 2’, ‘FIRST_GIANT_BRANCH = 3’, ‘CORE_HELIUM_BURNING = 4’, ‘EARLY_ASYMPTOTIC_GIANT_BRANCH = 5’, ‘THERMALLY_PULSING_ASYMPTOTIC_GIANT_BRANCH = 6’, ‘NAKED_MAIN_SEQUENCE_HELIUM_STAR = 7’, ‘NAKED_HELIUM_STAR_HERTZSPRUNG_GAP = 8’, ‘NAKED_HELIUM_STAR_GIANT_BRANCH = 9’, ‘HELIUM_WHITE_DWARF = 10’, ‘CARBON_OXYGEN_WHITE_DWARF = 11’, ‘OXYGEN_NEON_WHITE_DWARF = 12’, ‘NEUTRON_STAR = 13’, ‘BLACK_HOLE = 14’, ‘STAR_WITH_NO_MASS = 15’]

stellar_type_2

Description:

Set the stellar type of star 2 (internal index 1), usually MAIN_SEQUENCE (main sequence). Note that setting the stellar type only works for stars with both age=0 and core_mass=0, i.e. main sequence (hydrogen or helium), white dwarfs, black holes and neutrn stars.

Parameter input type:

Integer

Default value:

0

Macros:

[‘LOW_MASS_MS = 0’, ‘MS = 1’, ‘HG = 2’, ‘GIANT_BRANCH = 3’, ‘CHeB = 4’, ‘EAGB = 5’, ‘TPAGB = 6’, ‘HeMS = 7’, ‘HeHG = 8’, ‘HeGB = 9’, ‘HeWD = 10’, ‘COWD = 11’, ‘ONeWD = 12’, ‘NS = 13’, ‘BH = 14’, ‘MASSLESS_REMNANT = 15’, ‘LOW_MASS_MAIN_SEQUENCE = 0’, ‘MAIN_SEQUENCE = 1’, ‘HERTZSPRUNG_GAP = 2’, ‘FIRST_GIANT_BRANCH = 3’, ‘CORE_HELIUM_BURNING = 4’, ‘EARLY_ASYMPTOTIC_GIANT_BRANCH = 5’, ‘THERMALLY_PULSING_ASYMPTOTIC_GIANT_BRANCH = 6’, ‘NAKED_MAIN_SEQUENCE_HELIUM_STAR = 7’, ‘NAKED_HELIUM_STAR_HERTZSPRUNG_GAP = 8’, ‘NAKED_HELIUM_STAR_GIANT_BRANCH = 9’, ‘HELIUM_WHITE_DWARF = 10’, ‘CARBON_OXYGEN_WHITE_DWARF = 11’, ‘OXYGEN_NEON_WHITE_DWARF = 12’, ‘NEUTRON_STAR = 13’, ‘BLACK_HOLE = 14’, ‘STAR_WITH_NO_MASS = 15’]

stellar_type_3

Description:

Set the stellar type of star 3 (internal index 2), usually MAIN_SEQUENCE (main sequence). Note that setting the stellar type only works for stars with both age=0 and core_mass=0, i.e. main sequence (hydrogen or helium), white dwarfs, black holes and neutrn stars.

Parameter input type:

Integer

Default value:

0

Macros:

[‘LOW_MASS_MS = 0’, ‘MS = 1’, ‘HG = 2’, ‘GIANT_BRANCH = 3’, ‘CHeB = 4’, ‘EAGB = 5’, ‘TPAGB = 6’, ‘HeMS = 7’, ‘HeHG = 8’, ‘HeGB = 9’, ‘HeWD = 10’, ‘COWD = 11’, ‘ONeWD = 12’, ‘NS = 13’, ‘BH = 14’, ‘MASSLESS_REMNANT = 15’, ‘LOW_MASS_MAIN_SEQUENCE = 0’, ‘MAIN_SEQUENCE = 1’, ‘HERTZSPRUNG_GAP = 2’, ‘FIRST_GIANT_BRANCH = 3’, ‘CORE_HELIUM_BURNING = 4’, ‘EARLY_ASYMPTOTIC_GIANT_BRANCH = 5’, ‘THERMALLY_PULSING_ASYMPTOTIC_GIANT_BRANCH = 6’, ‘NAKED_MAIN_SEQUENCE_HELIUM_STAR = 7’, ‘NAKED_HELIUM_STAR_HERTZSPRUNG_GAP = 8’, ‘NAKED_HELIUM_STAR_GIANT_BRANCH = 9’, ‘HELIUM_WHITE_DWARF = 10’, ‘CARBON_OXYGEN_WHITE_DWARF = 11’, ‘OXYGEN_NEON_WHITE_DWARF = 12’, ‘NEUTRON_STAR = 13’, ‘BLACK_HOLE = 14’, ‘STAR_WITH_NO_MASS = 15’]

stellar_type_4

Description:

Set the stellar type of star 4 (internal index 3), usually MAIN_SEQUENCE (main sequence). Note that setting the stellar type only works for stars with both age=0 and core_mass=0, i.e. main sequence (hydrogen or helium), white dwarfs, black holes and neutrn stars.

Parameter input type:

Integer

Default value:

0

Macros:

[‘LOW_MASS_MS = 0’, ‘MS = 1’, ‘HG = 2’, ‘GIANT_BRANCH = 3’, ‘CHeB = 4’, ‘EAGB = 5’, ‘TPAGB = 6’, ‘HeMS = 7’, ‘HeHG = 8’, ‘HeGB = 9’, ‘HeWD = 10’, ‘COWD = 11’, ‘ONeWD = 12’, ‘NS = 13’, ‘BH = 14’, ‘MASSLESS_REMNANT = 15’, ‘LOW_MASS_MAIN_SEQUENCE = 0’, ‘MAIN_SEQUENCE = 1’, ‘HERTZSPRUNG_GAP = 2’, ‘FIRST_GIANT_BRANCH = 3’, ‘CORE_HELIUM_BURNING = 4’, ‘EARLY_ASYMPTOTIC_GIANT_BRANCH = 5’, ‘THERMALLY_PULSING_ASYMPTOTIC_GIANT_BRANCH = 6’, ‘NAKED_MAIN_SEQUENCE_HELIUM_STAR = 7’, ‘NAKED_HELIUM_STAR_HERTZSPRUNG_GAP = 8’, ‘NAKED_HELIUM_STAR_GIANT_BRANCH = 9’, ‘HELIUM_WHITE_DWARF = 10’, ‘CARBON_OXYGEN_WHITE_DWARF = 11’, ‘OXYGEN_NEON_WHITE_DWARF = 12’, ‘NEUTRON_STAR = 13’, ‘BLACK_HOLE = 14’, ‘STAR_WITH_NO_MASS = 15’]

max_stellar_type_1

Description:

The maximum stellar type of star 1 (internal index 0). Evolution is stopped when the star reaches this stellar type. If this is negative, massless remnants are allowed, and the maximum stellar type is the absolute value.

Parameter input type:

Integer

Default value:

16

Macros:

[‘LOW_MASS_MS = 0’, ‘MS = 1’, ‘HG = 2’, ‘GIANT_BRANCH = 3’, ‘CHeB = 4’, ‘EAGB = 5’, ‘TPAGB = 6’, ‘HeMS = 7’, ‘HeHG = 8’, ‘HeGB = 9’, ‘HeWD = 10’, ‘COWD = 11’, ‘ONeWD = 12’, ‘NS = 13’, ‘BH = 14’, ‘MASSLESS_REMNANT = 15’, ‘LOW_MASS_MAIN_SEQUENCE = 0’, ‘MAIN_SEQUENCE = 1’, ‘HERTZSPRUNG_GAP = 2’, ‘FIRST_GIANT_BRANCH = 3’, ‘CORE_HELIUM_BURNING = 4’, ‘EARLY_ASYMPTOTIC_GIANT_BRANCH = 5’, ‘THERMALLY_PULSING_ASYMPTOTIC_GIANT_BRANCH = 6’, ‘NAKED_MAIN_SEQUENCE_HELIUM_STAR = 7’, ‘NAKED_HELIUM_STAR_HERTZSPRUNG_GAP = 8’, ‘NAKED_HELIUM_STAR_GIANT_BRANCH = 9’, ‘HELIUM_WHITE_DWARF = 10’, ‘CARBON_OXYGEN_WHITE_DWARF = 11’, ‘OXYGEN_NEON_WHITE_DWARF = 12’, ‘NEUTRON_STAR = 13’, ‘BLACK_HOLE = 14’, ‘STAR_WITH_NO_MASS = 15’]

max_stellar_type_2

Description:

The maximum stellar type of star 2 (internal index 1). Evolution is stopped when the star reaches this stellar type. If this is negative, massless remnants are allowed, and the maximum stellar type is the absolute value.

Parameter input type:

Integer

Default value:

16

Macros:

[‘LOW_MASS_MS = 0’, ‘MS = 1’, ‘HG = 2’, ‘GIANT_BRANCH = 3’, ‘CHeB = 4’, ‘EAGB = 5’, ‘TPAGB = 6’, ‘HeMS = 7’, ‘HeHG = 8’, ‘HeGB = 9’, ‘HeWD = 10’, ‘COWD = 11’, ‘ONeWD = 12’, ‘NS = 13’, ‘BH = 14’, ‘MASSLESS_REMNANT = 15’, ‘LOW_MASS_MAIN_SEQUENCE = 0’, ‘MAIN_SEQUENCE = 1’, ‘HERTZSPRUNG_GAP = 2’, ‘FIRST_GIANT_BRANCH = 3’, ‘CORE_HELIUM_BURNING = 4’, ‘EARLY_ASYMPTOTIC_GIANT_BRANCH = 5’, ‘THERMALLY_PULSING_ASYMPTOTIC_GIANT_BRANCH = 6’, ‘NAKED_MAIN_SEQUENCE_HELIUM_STAR = 7’, ‘NAKED_HELIUM_STAR_HERTZSPRUNG_GAP = 8’, ‘NAKED_HELIUM_STAR_GIANT_BRANCH = 9’, ‘HELIUM_WHITE_DWARF = 10’, ‘CARBON_OXYGEN_WHITE_DWARF = 11’, ‘OXYGEN_NEON_WHITE_DWARF = 12’, ‘NEUTRON_STAR = 13’, ‘BLACK_HOLE = 14’, ‘STAR_WITH_NO_MASS = 15’]

max_stellar_type_3

Description:

The maximum stellar type of star 3 (internal index 2). Evolution is stopped when the star reaches this stellar type. If this is negative, massless remnants are allowed, and the maximum stellar type is the absolute value.

Parameter input type:

Integer

Default value:

16

Macros:

[‘LOW_MASS_MS = 0’, ‘MS = 1’, ‘HG = 2’, ‘GIANT_BRANCH = 3’, ‘CHeB = 4’, ‘EAGB = 5’, ‘TPAGB = 6’, ‘HeMS = 7’, ‘HeHG = 8’, ‘HeGB = 9’, ‘HeWD = 10’, ‘COWD = 11’, ‘ONeWD = 12’, ‘NS = 13’, ‘BH = 14’, ‘MASSLESS_REMNANT = 15’, ‘LOW_MASS_MAIN_SEQUENCE = 0’, ‘MAIN_SEQUENCE = 1’, ‘HERTZSPRUNG_GAP = 2’, ‘FIRST_GIANT_BRANCH = 3’, ‘CORE_HELIUM_BURNING = 4’, ‘EARLY_ASYMPTOTIC_GIANT_BRANCH = 5’, ‘THERMALLY_PULSING_ASYMPTOTIC_GIANT_BRANCH = 6’, ‘NAKED_MAIN_SEQUENCE_HELIUM_STAR = 7’, ‘NAKED_HELIUM_STAR_HERTZSPRUNG_GAP = 8’, ‘NAKED_HELIUM_STAR_GIANT_BRANCH = 9’, ‘HELIUM_WHITE_DWARF = 10’, ‘CARBON_OXYGEN_WHITE_DWARF = 11’, ‘OXYGEN_NEON_WHITE_DWARF = 12’, ‘NEUTRON_STAR = 13’, ‘BLACK_HOLE = 14’, ‘STAR_WITH_NO_MASS = 15’]

max_stellar_type_4

Description:

The maximum stellar type of star 4 (internal index 3). Evolution is stopped when the star reaches this stellar type. If this is negative, massless remnants are allowed, and the maximum stellar type is the absolute value.

Parameter input type:

Integer

Default value:

16

Macros:

[‘LOW_MASS_MS = 0’, ‘MS = 1’, ‘HG = 2’, ‘GIANT_BRANCH = 3’, ‘CHeB = 4’, ‘EAGB = 5’, ‘TPAGB = 6’, ‘HeMS = 7’, ‘HeHG = 8’, ‘HeGB = 9’, ‘HeWD = 10’, ‘COWD = 11’, ‘ONeWD = 12’, ‘NS = 13’, ‘BH = 14’, ‘MASSLESS_REMNANT = 15’, ‘LOW_MASS_MAIN_SEQUENCE = 0’, ‘MAIN_SEQUENCE = 1’, ‘HERTZSPRUNG_GAP = 2’, ‘FIRST_GIANT_BRANCH = 3’, ‘CORE_HELIUM_BURNING = 4’, ‘EARLY_ASYMPTOTIC_GIANT_BRANCH = 5’, ‘THERMALLY_PULSING_ASYMPTOTIC_GIANT_BRANCH = 6’, ‘NAKED_MAIN_SEQUENCE_HELIUM_STAR = 7’, ‘NAKED_HELIUM_STAR_HERTZSPRUNG_GAP = 8’, ‘NAKED_HELIUM_STAR_GIANT_BRANCH = 9’, ‘HELIUM_WHITE_DWARF = 10’, ‘CARBON_OXYGEN_WHITE_DWARF = 11’, ‘OXYGEN_NEON_WHITE_DWARF = 12’, ‘NEUTRON_STAR = 13’, ‘BLACK_HOLE = 14’, ‘STAR_WITH_NO_MASS = 15’]

probability

Description:

The probability is a weighting applied to the star based on, say, the initial mass function. When running a grid of stars to simulate all stars, the summed probability of all the stars should be 1.0.

Parameter input type:

Float

Default value:

1

phasevol

Description:

The system’s phase volume, used by binary_grid.

Parameter input type:

Float

Default value:

NULL

stellar_structure_algorithm

Description:

Set the stellar structure algorithm. 0=modified BSE (default), 1=none, 2=external function (must be defined by the calling code), 3=binary_c (not yet implemented).

Parameter input type:

Integer

Default value:

0

Macros:

[‘STELLAR_STRUCTURE_ALGORITHM_MODIFIED_BSE = 0’, ‘STELLAR_STRUCTURE_ALGORITHM_NONE = 1’, ‘STELLAR_STRUCTURE_ALGORITHM_EXTERNAL_FUNCTION = 2’, ‘STELLAR_STRUCTURE_ALGORITHM_MINT = 3’]

solver

Description:

The type of solver. Default is the Forward-Euler (0), but could be RK2 (1), RK4 (2) or a predictor-corretor (3).

Parameter input type:

Integer

Default value:

0

Macros:

[‘SOLVER_FORWARD_EULER = 0’, ‘SOLVER_RK2 = 1’, ‘SOLVER_RK4 = 2’, ‘SOLVER_PREDICTOR_CORRECTOR = 3’]

max_evolution_time

Description:

Set the maximum age for the stars (Myr).

Parameter input type:

Float

Default value:

15000

max_model_number

Description:

Set the maximum number of models, ignored if 0 (default is 0).

Parameter input type:

Integer

Default value:

0

monte_carlo_kicks

Description:

Turn on Monte-Carlo SN kicks. On (True) by default, and indeed other algorithms are probably broken.

Parameter input type:

True|False

Default value:

True

disable_debug

Description:

Disables debug output. Only has an effect when DEBUG is 1, which probably requires a rebuild. Default FALSE.

Parameter input type:

True|False

Default value:

False

debug_filter_pointers

Description:

Filter pointers out of debug output, so that “0x1234567” becomes “0xpointer”.

Parameter input type:

True|False

Default value:

False

timestep_logging

Description:

Turn on timestep logging (default is False).

Parameter input type:

True|False

Default value:

False

flush_log

Description:

Turn on flushing of the log file (default is False, which is faster).

Parameter input type:

True|False

Default value:

False

derivative_logging

Description:

Turn on derivative logging (default is False).

Parameter input type:

True|False

Default value:

False

progress_bar

Description:

Turn on progress-bar logging (default is False).

Parameter input type:

True|False

Default value:

False

log_all_reject_timestep_failures

Description:

Turn on logging of every timestep that is rejected but which cannot have its timestep shortened.

Parameter input type:

True|False

Default value:

False

cannot_shorten_timestep_policy

Description:

Turn on logging of every timestep that is rejected but which cannot have its timestep shortened.

Parameter input type:

Integer

Default value:

0

Macros:

[‘CANNOT_SHORTEN_RESTORE_AND_TRY_EVENTS = 0’, ‘CANNOT_SHORTEN_CONTINUE = 1’, ‘CANNOT_SHORTEN_FAIL = 2’]

rejects_in_log

Description:

Show timestep rejections in the main log (default is False).

Parameter input type:

True|False

Default value:

False

vandenHeuvel_logging

Description:

Turn on van den Heuvel logging (default is False).

Parameter input type:

True|False

Default value:

False

evolution_splitting

Description:

If True, turn on splitting of an evolutionary run if splitpoint (e.g. supernova) occurs.

Parameter input type:

True|False

Default value:

False

disable_events

Description:

Whether to disable the new events code (defaults to False, so we use events by default).

Parameter input type:

True|False

Default value:

False

evolution_splitting_sn_eccentricity_threshold

Description:

Threshold eccentricity above which evolution splitting happens in a system with no SN kick. (0.01).

Parameter input type:

Float

Default value:

0.01

evolution_splitting_sn_n

Description:

Number of runs to split into when a SN occurs.

Parameter input type:

Integer

Default value:

10

evolution_splitting_maxdepth

Description:

Max number of splits in an evolutionary run.

Parameter input type:

Integer

Default value:

2

equation_of_state_algorithm

Description:

Set the equation of state algorithm. 0 = Paczynski.

Parameter input type:

Integer

Default value:

NULL

Macros:

[‘EQUATION_OF_STATE_PACZYNSKI = 0’]

opacity_algorithm

Description:

Set the opacity algorithm. 0 = Paczynski, 1 = Ferguson/Opal.

Parameter input type:

Integer

Default value:

NULL

Macros:

[‘OPACITY_ALGORITHM_PACZYNSKI = 0’, ‘OPACITY_ALGORITHM_FERGUSON_OPAL = 1’, ‘OPACITY_ALGORITHM_STARS = 2’]

wind_mass_loss

Description:

Defines the algorithm used for stellar winds. 0 = none, 1 = Hurley et al. (2002), 2 = Schneider (2018).

Parameter input type:

Unsigned integer

Default value:

3

Macros:

[‘WIND_ALGORITHM_NONE = 0’, ‘WIND_ALGORITHM_HURLEY2002 = 1’, ‘WIND_ALGORITHM_SCHNEIDER2018 = 2’, ‘WIND_ALGORITHM_BINARY_C_2020 = 3’, ‘WIND_ALGORITHM_BINARY_C_2022 = 4’]

Extra:

0

gbwind

Description:

Wind prescription for first red giant branch stars. 0=Reimers (Hurley et al 2000/2002; choose gb_reimers_eta=0.5 for their mass loss rate), 1=Schroeder+Cuntz 2005 (set gb_reimers_eta=1.0 for their mass loss rate).

Parameter input type:

Integer

Default value:

0

Macros:

[‘GB_WIND_REIMERS = 0’, ‘GB_WIND_SCHROEDER_CUNTZ_2005 = 1’, ‘GB_WIND_GOLDMAN_ETAL_2017 = 2’, ‘GB_WIND_BEASOR_ETAL_2020 = 3’]

postagbwind

Description:

Apply special post-(A)GB prescription. Default is POSTAGB_WIND_USE_GIANT which means we just use whatever is prescribed on the giant branch. Other options include: POSTAGB_WIND_NONE = 1 (no wind on the post (A)GB), POSTAGB_WIND_KRTICKA2020 = 2 which uses Krticka, Kubát and Krticková (2020, A&A 635, A173).

Parameter input type:

Integer

Default value:

0

Macros:

[‘POSTAGB_WIND_GIANT = 0’, ‘POSTAGB_WIND_NONE = 1’, ‘POSTAGB_WIND_KRTICKA2020 = 2’]

Teff_postAGB_min

Description:

The minimum temperature for which we apply post-(A)GB winds. See also Teff_postAGB_max. (6000 K).

Parameter input type:

Float

Default value:

6000

Teff_postAGB_max

Description:

The maximum temperature for which we apply post-(A)GB winds. See also Teff_postAGB_min. (120000 K).

Parameter input type:

Float

Default value:

120000

mattsson_Orich_tpagbwind

Description:

Experimental: turns on Mattsson’s TPAGB wind when the star is oxygen rich. Requires MATTSSON_MASS_LOSS.

Parameter input type:

Integer

Default value:

NULL

magnetic_braking_factor

Description:

Multiplier for the magnetic braking angular momentum loss rate.

Parameter input type:

Float

Default value:

1

magnetic_braking_gamma

Description:

Gamma factor in Rappaport style magnetic braking expression.

Parameter input type:

Float

Default value:

3

magnetic_braking_algorithm

Description:

Algorithm for the magnetic braking angular momentum loss rate. 0 = Hurley et al. 2002, 1 = Andronov, Pinnsonneault and Sills 2003, 2 = Barnes and Kim 2010.

Parameter input type:

Integer

Default value:

0

Macros:

[‘MAGNETIC_BRAKING_ALGORITHM_HURLEY_2002 = 0’, ‘MAGNETIC_BRAKING_ALGORITHM_ANDRONOV_2003 = 1’, ‘MAGNETIC_BRAKING_ALGORITHM_BARNES_2010 = 2’, ‘MAGNETIC_BRAKING_ALGORITHM_RAPPAPORT_1983 = 3’]

helium_flash_mass_loss

Description:

Mass to be lost at the helium flash.

Parameter input type:

Float

Default value:

0

gb_reimers_eta

Description:

First red giant branch wind multiplication factor, cf. eta in Reimers’ mass loss formula. (This multiplies the 4e-13 in Reimers’ formula, or the 8e-14 in Schroeder and Cuntz.).

Parameter input type:

Float

Default value:

0.4

gbwindfac

Description:

Multiplier for the giant branch wind mass loss rate.

Parameter input type:

Float

Default value:

1

tpagbwindfac

Description:

Multiplier for the TPAGB wind mass loss rate.

Parameter input type:

Float

Default value:

1

eagbwindfac

Description:

Multiplier for the EAGB wind mass loss rate.

Parameter input type:

Float

Default value:

1

nieuwenhuijzen_windfac

Description:

Multiplier for the Nieuwenhuijzen & de Jager wind mass loss rate.

Parameter input type:

Float

Default value:

1

tpagbwind

Description:

Wind prescription during the TPAGB. 0=Karakas 2002 (a modified Vassiliadis and Wood 1993), 1=Hurley et al 2000/2002 (Vassiliadis and Wood 1993), 2=Reimers, 3=Bloecker, 4=Van Loon, 5=Rob’s C-wind (broken?), 6,7=Vassiliadis and Wood 1993 (Karakas,Hurley variants respectively) when C/O>1, 8=Mattsson, 9 = Goldman et al. (2017), 10 = Beasor et al. (2020).

Parameter input type:

Integer

Default value:

0

Macros:

[‘TPAGB_WIND_VW93_KARAKAS = 0’, ‘TPAGB_WIND_VW93_ORIG = 1’, ‘TPAGB_WIND_REIMERS = 2’, ‘TPAGB_WIND_BLOECKER = 3’, ‘TPAGB_WIND_VAN_LOON = 4’, ‘TPAGB_WIND_ROB_CWIND = 5’, ‘TPAGB_WIND_VW93_KARAKAS_CARBON_STARS = 6’, ‘TPAGB_WIND_VW93_ORIG_CARBON_STARS = 7’, ‘TPAGB_WIND_MATTSSON = 8’, ‘TPAGB_WIND_GOLDMAN_ETAL_2017 = 9’, ‘TPAGB_WIND_BEASOR_ETAL_2020 = 10’]

eagbwind

Description:

Wind prescription during the EAGB. 0=BSE (Hurley+2002, based on VW93), 1 = Goldman et al. (2017), 2 = Beasor et al. (2020).

Parameter input type:

Integer

Default value:

0

Macros:

[‘EAGB_WIND_BSE = 0’, ‘EAGB_WIND_GOLDMAN_ETAL_2017 = 1’, ‘EAGB_WIND_BEASOR_ETAL_2020 = 2’]

wind_gas_to_dust_ratio

Description:

Gas to dust ratio used in wind calculations (e.g. Goldman et al. 2017). Typically 200 (Milky Way)-500 (Magellanic Clouds). Default is 200, approximately as in MW stars.

Parameter input type:

Float

Default value:

200

vwind_multiplier

Description:

Multiplier for the stellar wind velocity.

Parameter input type:

Float

Default value:

1

vwind_beta

Description:

Beta for stellar wind speed calculations, where vwind=sqrt(beta) * escape velocity. Default 0.125 (from BSE, Hurley et al. 2002).

Parameter input type:

Float

Default value:

0.125

superwind_mira_switchon

Description:

In the Vassiliadis and Wood (1993) AGB wind prescription, the superwind is turned on at a given Mira period, usually 500 days. You can vary when this switch-on happens with this parameter.

Parameter input type:

Float

Default value:

500

vw93_mira_shift

Description:

In the Vassiliadis and Wood (1993) AGB wind prescription, the wind loss rate depends on the Mira period plus this offset. Requires VW93_MIRA_SHIFT.

Parameter input type:

Float

Default value:

0

vw93_multiplier

Description:

In the Vassiliadis and Wood (1993) AGB wind prescription, the wind loss rate is multiplied by this factor. Requires VW93_MULTIPLIER.

Parameter input type:

Float

Default value:

1

tpagb_reimers_eta

Description:

TPAGB Reimers wind multiplication factor, cf. eta in Reimers’ mass loss formula. (This multiplies the 4e-13 in Reimers’ formula, or the 8e-14 in Schroeder and Cuntz.) Note that Reimers is not the default TPAGB wind prescription. See also tpagbwind.

Parameter input type:

Float

Default value:

1

Tout_Pringle_1992_multiplier

Description:

Multiplier for the Tout & Pringle (1992) magnetic wind. (0.0).

Parameter input type:

Float

Default value:

0

artificial_mass_accretion_rate%d

Description:

Constant mass accretion rate for star <n>.

Parameter input type:

Float(scanf)

Default value:

NULL

artificial_mass_accretion_rate_by_stellar_type%d

Description:

Constant mass accretion rate for stellar type <n>.

Parameter input type:

Float(scanf)

Default value:

NULL

artificial_mass_accretion_rate_by_stellar_typeLOW_MASS_MS

Description:

Constant mass accretion rate for stellar type LOW_MASS_MS.

artificial_mass_accretion_rate_by_stellar_typeMS

Description:

Constant mass accretion rate for stellar type MS.

artificial_mass_accretion_rate_by_stellar_typeHG

Description:

Constant mass accretion rate for stellar type HG.

artificial_mass_accretion_rate_by_stellar_typeGIANT_BRANCH

Description:

Constant mass accretion rate for stellar type GIANT_BRANCH.

artificial_mass_accretion_rate_by_stellar_typeCHeB

Description:

Constant mass accretion rate for stellar type CHeB.

artificial_mass_accretion_rate_by_stellar_typeEAGB

Description:

Constant mass accretion rate for stellar type EAGB.

artificial_mass_accretion_rate_by_stellar_typeTPAGB

Description:

Constant mass accretion rate for stellar type TPAGB.

artificial_mass_accretion_rate_by_stellar_typeHeMS

Description:

Constant mass accretion rate for stellar type HeMS.

artificial_mass_accretion_rate_by_stellar_typeHeHG

Description:

Constant mass accretion rate for stellar type HeHG.

artificial_mass_accretion_rate_by_stellar_typeHeGB

Description:

Constant mass accretion rate for stellar type HeGB.

artificial_mass_accretion_rate_by_stellar_typeHeWD

Description:

Constant mass accretion rate for stellar type HeWD.

artificial_mass_accretion_rate_by_stellar_typeCOWD

Description:

Constant mass accretion rate for stellar type COWD.

artificial_mass_accretion_rate_by_stellar_typeONeWD

Description:

Constant mass accretion rate for stellar type ONeWD.

artificial_mass_accretion_rate_by_stellar_typeNS

Description:

Constant mass accretion rate for stellar type NS.

artificial_mass_accretion_rate_by_stellar_typeBH

Description:

Constant mass accretion rate for stellar type BH.

artificial_mass_accretion_rate_by_stellar_typeMASSLESS_REMNANT

Description:

Constant mass accretion rate for stellar type MASSLESS_REMNANT.

artificial_angular_momentum_accretion_rate%d

Description:

Constant angular momentum accretion for star <n>.

Parameter input type:

Float(scanf)

Default value:

NULL

artificial_orbital_angular_momentum_accretion_rate

Description:

Constant angular momentum accretion rate on the orbit.

Parameter input type:

Float

Default value:

0

artificial_accretion_start_time

Description:

Time at which artificial accretion stars. Ignored if <0 (default is -1).

Parameter input type:

Float

Default value:

-1

Macros:

[‘ARTIFICIAL_ACCRETION_IGNORE = -1’]

artificial_accretion_end_time

Description:

Time at which artificial accretion ends. Ignored if <0 (default is -1).

Parameter input type:

Float

Default value:

-1

Macros:

[‘ARTIFICIAL_ACCRETION_IGNORE = -1’]

wr_wind

Description:

Massive-star (WR) wind prescription. 0 = Hurley et al 2000/2002, 1=Maeder and Meynet, 2=Nugis and Lamers, 3=John Eldridge’s version of Vink’s early-2000s wind (See Lynnette Dray’s thesis, or John Eldridge’s thesis).

Parameter input type:

Integer

Default value:

0

Macros:

[‘WR_WIND_BSE = 0’, ‘WR_WIND_MAEDER_MEYNET = 1’, ‘WR_WIND_NUGIS_LAMERS = 2’, ‘WR_WIND_ELDRIDGE = 3’]

wr_wind_fac

Description:

Massive-star (WR) wind multiplication factor.

Parameter input type:

Float

Default value:

1

wrwindfac

Description:

Massive-star (WR) wind multiplication factor. Synonymous with wr_wind_fac (which you should use instead).

Parameter input type:

Float

Default value:

1

BH_prescription

Description:

Black hole mass prescrition: relates the mass of a newly formed black hole to its progenitor’s (CO) core mass. BH_HURLEY2002 = 0 = Hurley et al 2000/2002, BH_BELCZYNSKI = 1 = Belczynski (early 2000s), BH_SPERA2015 = Spera+ 2015, BH_FRYER12_DELAYED = 3 = Fryer et al. (2012) delayed prescription, BH_FRYER12_RAPID = 4 = Fryer et al. (2012) rapid prescription, BH_FRYER12_STARTRACK = 5 = Fryer et al. (2012) startrack prescription, BH_DARK = 6 = set BH mass to the CO core mass so no mass is ejected, BH_MARASSI2019 = 7 use Marassi 2019 prescription.

Parameter input type:

Integer

Default value:

0

Macros:

[‘BH_HURLEY2002 = 0’, ‘BH_BELCZYNSKI = 1’, ‘BH_SPERA2015 = 2’, ‘BH_FRYER12_DELAYED = 3’, ‘BH_FRYER12_RAPID = 4’, ‘BH_FRYER12_STARTRACK = 5’, ‘BH_DARK = 6’, ‘BH_MARASSI2019 = 7’]

sn_kick_distribution_II

Description:

Set the distribution of speeds applied to kick type II core collapse supernova systems. 0=fixed, 1=maxwellian (hurley/BSE), 2=custom function (see monte_carlo_kicks.c).

Parameter input type:

Integer

Default value:

1

Macros:

[‘KICK_VELOCITY_FIXED = 0’, ‘KICK_VELOCITY_MAXWELLIAN = 1’, ‘KICK_VELOCITY_CUSTOM = 2’, ‘KICK_VELOCITY_GIACOBBO_MAPELLI_2020 = 3’]

sn_kick_distribution_ECAP

Description:

Set the distribution of speeds applied to the remnants of electron-capture supernovae. 0=fixed, 1=maxwellian (hurley/BSE), 2=custom function (see monte_carlo_kicks.c).

Parameter input type:

Integer

Default value:

0

Macros:

[‘KICK_VELOCITY_FIXED = 0’, ‘KICK_VELOCITY_MAXWELLIAN = 1’, ‘KICK_VELOCITY_CUSTOM = 2’, ‘KICK_VELOCITY_GIACOBBO_MAPELLI_2020 = 3’]

sn_kick_distribution_NS_NS

Description:

Set the distribution of speeds applied to kick neutron stars and black holes that survive a NS-NS merger. 0=fixed, 1=maxwellian (hurley/BSE), 2=custom function (see monte_carlo_kicks.c).

Parameter input type:

Integer

Default value:

0

Macros:

[‘KICK_VELOCITY_FIXED = 0’, ‘KICK_VELOCITY_MAXWELLIAN = 1’, ‘KICK_VELOCITY_CUSTOM = 2’, ‘KICK_VELOCITY_GIACOBBO_MAPELLI_2020 = 3’]

sn_kick_distribution_IBC

Description:

Set the distribution of speeds applied to kick newly-born neutron stars and black holes after a type Ib/c core-collapse supernova. 0=fixed, 1=maxwellian (hurley/BSE), 2=custom function (see monte_carlo_kicks.c).

Parameter input type:

Integer

Default value:

1

Macros:

[‘KICK_VELOCITY_FIXED = 0’, ‘KICK_VELOCITY_MAXWELLIAN = 1’, ‘KICK_VELOCITY_CUSTOM = 2’, ‘KICK_VELOCITY_GIACOBBO_MAPELLI_2020 = 3’]

sn_kick_distribution_GRB_COLLAPSAR

Description:

Set the distribution of speeds applied to kick newly-born neutron stars and black holes after a type Ib/c core-collapse supernova which is also a collapsar. 0=fixed, 1=maxwellian (hurley/BSE), 2=custom function (see monte_carlo_kicks.c).

Parameter input type:

Integer

Default value:

1

Macros:

[‘KICK_VELOCITY_FIXED = 0’, ‘KICK_VELOCITY_MAXWELLIAN = 1’, ‘KICK_VELOCITY_CUSTOM = 2’, ‘KICK_VELOCITY_GIACOBBO_MAPELLI_2020 = 3’]

sn_kick_distribution_TZ

Description:

Set the distribution of speeds applied to kick newly-born neutron stars and black holes at the death of a Thorne-Zytkow object. 0=fixed, 1=maxwellian (hurley/BSE), 2=custom function (see monte_carlo_kicks.c).

Parameter input type:

Integer

Default value:

0

Macros:

[‘KICK_VELOCITY_FIXED = 0’, ‘KICK_VELOCITY_MAXWELLIAN = 1’, ‘KICK_VELOCITY_CUSTOM = 2’, ‘KICK_VELOCITY_GIACOBBO_MAPELLI_2020 = 3’]

sn_kick_distribution_AIC_BH

Description:

Set the distribution of speeds applied to kick newly-born neutron stars black holes after accretion induced collapse of a neutron star. 0=fixed, 1=maxwellian (hurley/BSE), 2=custom function (see monte_carlo_kicks.c).

Parameter input type:

Integer

Default value:

0

Macros:

[‘KICK_VELOCITY_FIXED = 0’, ‘KICK_VELOCITY_MAXWELLIAN = 1’, ‘KICK_VELOCITY_CUSTOM = 2’, ‘KICK_VELOCITY_GIACOBBO_MAPELLI_2020 = 3’]

sn_kick_distribution_BH_BH

Description:

Set the distribution of speeds applied to black holes formed by the merger of two black holes. 0=fixed, 1=maxwellian (hurley/BSE), 2=custom function (see monte_carlo_kicks.c).

Parameter input type:

Integer

Default value:

0

Macros:

[‘KICK_VELOCITY_FIXED = 0’, ‘KICK_VELOCITY_MAXWELLIAN = 1’, ‘KICK_VELOCITY_CUSTOM = 2’, ‘KICK_VELOCITY_GIACOBBO_MAPELLI_2020 = 3’]

sn_kick_distribution_BH_NS

Description:

Set the distribution of speeds applied to black holes formed by the merger of a neutron star and a black holes. 0=fixed, 1=maxwellian (hurley/BSE), 2=custom function (see monte_carlo_kicks.c).

Parameter input type:

Integer

Default value:

0

Macros:

[‘KICK_VELOCITY_FIXED = 0’, ‘KICK_VELOCITY_MAXWELLIAN = 1’, ‘KICK_VELOCITY_CUSTOM = 2’, ‘KICK_VELOCITY_GIACOBBO_MAPELLI_2020 = 3’]

sn_kick_distribution_IA_Hybrid_HeCOWD

Description:

Set the distribution of speeds applied to any survivor of a hybrid He-COWD SNIa explosion. 0=fixed, 1=maxwellian (hurley/BSE), 2=custom function (see monte_carlo_kicks.c).

Parameter input type:

Integer

Default value:

0

Macros:

[‘KICK_VELOCITY_FIXED = 0’, ‘KICK_VELOCITY_MAXWELLIAN = 1’, ‘KICK_VELOCITY_CUSTOM = 2’, ‘KICK_VELOCITY_GIACOBBO_MAPELLI_2020 = 3’]

sn_kick_distribution_IA_Hybrid_HeCOWD_subluminous

Description:

Set the distribution of speeds applied to any survivor of a subluminous hybrid He-COWD SNIa explosion. 0=fixed, 1=maxwellian (hurley/BSE), 2=custom function (see monte_carlo_kicks.c).

Parameter input type:

Integer

Default value:

0

Macros:

[‘KICK_VELOCITY_FIXED = 0’, ‘KICK_VELOCITY_MAXWELLIAN = 1’, ‘KICK_VELOCITY_CUSTOM = 2’, ‘KICK_VELOCITY_GIACOBBO_MAPELLI_2020 = 3’]

sn_kick_dispersion_II

Description:

Set the dispersion of speeds applied to kick type II core collapse supernova systems. 0=fixed, 1=maxwellian (hurley/BSE), 2=custom function (see monte_carlo_kicks.c).

Parameter input type:

Float

Default value:

190

sn_kick_dispersion_ECAP

Description:

Set the dispersion of speeds applied to the remnants of electron-capture supernovae. 0=fixed, 1=maxwellian (hurley/BSE), 2=custom function (see monte_carlo_kicks.c).

Parameter input type:

Float

Default value:

0

sn_kick_dispersion_NS_NS

Description:

Set the dispersion of speeds applied to kick neutron stars and black holes that survive a NS-NS merger. 0=fixed, 1=maxwellian (hurley/BSE), 2=custom function (see monte_carlo_kicks.c).

Parameter input type:

Float

Default value:

0

sn_kick_dispersion_IBC

Description:

Set the dispersion of speeds applied to kick newly-born neutron stars and black holes after a type Ib/c core-collapse supernova. 0=fixed, 1=maxwellian (hurley/BSE), 2=custom function (see monte_carlo_kicks.c).

Parameter input type:

Float

Default value:

190

sn_kick_dispersion_GRB_COLLAPSAR

Description:

Set the dispersion of speeds applied to kick newly-born neutron stars and black holes after a type Ib/c core-collapse supernova which is also a collapsar. 0=fixed, 1=maxwellian (hurley/BSE), 2=custom function (see monte_carlo_kicks.c).

Parameter input type:

Float

Default value:

190

sn_kick_dispersion_TZ

Description:

Set the dispersion of speeds applied to kick newly-born neutron stars and black holes at the death of a Thorne-Zytkow object. 0=fixed, 1=maxwellian (hurley/BSE), 2=custom function (see monte_carlo_kicks.c).

Parameter input type:

Float

Default value:

0

sn_kick_dispersion_AIC_BH

Description:

Set the dispersion of speeds applied to kick newly-born neutron stars black holes after accretion induced collapse of a neutron star. 0=fixed, 1=maxwellian (hurley/BSE), 2=custom function (see monte_carlo_kicks.c).

Parameter input type:

Float

Default value:

0

sn_kick_dispersion_BH_BH

Description:

Set the dispersion of speeds applied to black holes formed by the merger of two black holes. 0=fixed, 1=maxwellian (hurley/BSE), 2=custom function (see monte_carlo_kicks.c).

Parameter input type:

Float

Default value:

0

sn_kick_dispersion_BH_NS

Description:

Set the dispersion of speeds applied to black holes formed by the merger of a neutron star and a black holes. 0=fixed, 1=maxwellian (hurley/BSE), 2=custom function (see monte_carlo_kicks.c).

Parameter input type:

Float

Default value:

0

sn_kick_dispersion_IA_Hybrid_HeCOWD

Description:

Set the dispersion of speeds applied to the survivor of a SNIa explosion of a hybrid He-COWD. 0=fixed, 1=maxwellian (hurley/BSE), 2=custom function (see monte_carlo_kicks.c).

Parameter input type:

Float

Default value:

0

sn_kick_dispersion_IA_Hybrid_HeCOWD_subluminous

Description:

Set the dispersion of speeds applied to the survivor of a subluminous SNIa explosion of a hybrid He-COWD. 0=fixed, 1=maxwellian (hurley/BSE), 2=custom function (see monte_carlo_kicks.c).

Parameter input type:

Float

Default value:

0

sn_kick_companion_IA_He

Description:

Set the speed (if >=0) of, or the algothim (if <0) used to calculate the, kick on the companion when a Ia He supernova occurs. 0 = none, 1 = Liu+2015, 2 = Wheeler+ 1975.

Parameter input type:

Float

Default value:

0

Macros:

[‘SN_IMPULSE_NONE = 0’, ‘SN_IMPULSE_LIU2015 = 1’, ‘SN_IMPULSE_WHEELER1975 = 2’]

sn_kick_companion_IA_ELD

Description:

Set the speed (if >=0) of, or the algothim (if <0) used to calculate the, kick on the companion when a Ia ELD (sub-Mch) supernova occurs. 0 = none, 1 = Liu+2015, 2 = Wheeler+ 1975.

Parameter input type:

Float

Default value:

0

Macros:

[‘SN_IMPULSE_NONE = 0’, ‘SN_IMPULSE_LIU2015 = 1’, ‘SN_IMPULSE_WHEELER1975 = 2’]

sn_kick_companion_IA_CHAND

Description:

Set the speed (if >=0) of, or the algothim (if <0) used to calculate the, kick on the companion when a Ia Mch supernova occurs. 0 = none, 1 = Liu+2015, 2 = Wheeler+ 1975.

Parameter input type:

Float

Default value:

0

Macros:

[‘SN_IMPULSE_NONE = 0’, ‘SN_IMPULSE_LIU2015 = 1’, ‘SN_IMPULSE_WHEELER1975 = 2’]

sn_kick_companion_AIC

Description:

Set the speed (if >=0) of, or the algothim (if <0) used to calculate the, kick on the companion when an accretion induced collapse (supernova) occurs. 0 = none, 1 = Liu+2015, 2 = Wheeler+ 1975.

Parameter input type:

Float

Default value:

0

Macros:

[‘SN_IMPULSE_NONE = 0’, ‘SN_IMPULSE_LIU2015 = 1’, ‘SN_IMPULSE_WHEELER1975 = 2’]

sn_kick_companion_ECAP

Description:

Set the speed (if >=0) of, or the algothim (if <0) used to calculate the, kick on the companion when an electron capture supernova occurs. 0 = none, 1 = Liu+2015, 2 = Wheeler+ 1975.

Parameter input type:

Float

Default value:

0

Macros:

[‘SN_IMPULSE_NONE = 0’, ‘SN_IMPULSE_LIU2015 = 1’, ‘SN_IMPULSE_WHEELER1975 = 2’]

sn_kick_companion_IA_He_Coal

Description:

Set the speed (if >=0) of, or the algothim (if <0) used to calculate the, kick on the companion when a Ia helium merger supernova occurs. 0 = none, 1 = Liu+2015, 2 = Wheeler+ 1975.

Parameter input type:

Float

Default value:

0

Macros:

[‘SN_IMPULSE_NONE = 0’, ‘SN_IMPULSE_LIU2015 = 1’, ‘SN_IMPULSE_WHEELER1975 = 2’]

sn_kick_companion_IA_CHAND_Coal

Description:

Set the speed (if >=0) of, or the algothim (if <0) used to calculate the, kick on the companion when a Ia Mch merger supernova occurs. 0 = none, 1 = Liu+2015, 2 = Wheeler+ 1975.

Parameter input type:

Float

Default value:

0

Macros:

[‘SN_IMPULSE_NONE = 0’, ‘SN_IMPULSE_LIU2015 = 1’, ‘SN_IMPULSE_WHEELER1975 = 2’]

sn_kick_companion_NS_NS

Description:

Set the speed (if >=0) of, or the algothim (if <0) used to calculate the, kick on the companion when a neutron-star neutron-star merger. 0 = none, 1 = Liu+2015, 2 = Wheeler+ 1975.

Parameter input type:

Float

Default value:

0

Macros:

[‘SN_IMPULSE_NONE = 0’, ‘SN_IMPULSE_LIU2015 = 1’, ‘SN_IMPULSE_WHEELER1975 = 2’]

sn_kick_companion_GRB_COLLAPSAR

Description:

Set the speed (if >=0) of, or the algothim (if <0) used to calculate the, kick on the companion when a GRB Collapsar (rapidly rotating SN Ibc) supernova occurs. 0 = none, 1 = Liu+2015, 2 = Wheeler+ 1975.

Parameter input type:

Float

Default value:

0

Macros:

[‘SN_IMPULSE_NONE = 0’, ‘SN_IMPULSE_LIU2015 = 1’, ‘SN_IMPULSE_WHEELER1975 = 2’]

sn_kick_companion_HeStarIa

Description:

Set the speed (if >=0) of, or the algothim (if <0) used to calculate the, kick on the companion when a He-star Ia supernova occurs. 0 = none, 1 = Liu+2015, 2 = Wheeler+ 1975.

Parameter input type:

Float

Default value:

0

Macros:

[‘SN_IMPULSE_NONE = 0’, ‘SN_IMPULSE_LIU2015 = 1’, ‘SN_IMPULSE_WHEELER1975 = 2’]

sn_kick_companion_IBC

Description:

Set the speed (if >=0) of, or the algothim (if <0) used to calculate the, kick on the companion when a type Ib/c supernova occurs. 0 = none, 1 = Liu+2015, 2 = Wheeler+ 1975.

Parameter input type:

Float

Default value:

0

Macros:

[‘SN_IMPULSE_NONE = 0’, ‘SN_IMPULSE_LIU2015 = 1’, ‘SN_IMPULSE_WHEELER1975 = 2’]

sn_kick_companion_II

Description:

Set the speed (if >=0) of, or the algothim (if <0) used to calculate the, kick on the companion when a type II supernova occurs. 0 = none, 1 = Liu+2015, 2 = Wheeler+ 1975.

Parameter input type:

Float

Default value:

0

Macros:

[‘SN_IMPULSE_NONE = 0’, ‘SN_IMPULSE_LIU2015 = 1’, ‘SN_IMPULSE_WHEELER1975 = 2’]

sn_kick_companion_IIa

Description:

Set the speed (if >=0) of, or the algothim (if <0) used to calculate the, kick on the companion when a type IIa supernova occurs. 0 = none, 1 = Liu+2015, 2 = Wheeler+ 1975.

Parameter input type:

Float

Default value:

0

Macros:

[‘SN_IMPULSE_NONE = 0’, ‘SN_IMPULSE_LIU2015 = 1’, ‘SN_IMPULSE_WHEELER1975 = 2’]

sn_kick_companion_WDKICK

Description:

Set the speed (if >=0) of, or the algothim (if <0) used to calculate the, kick on the companion when a WD is kicked. 0 = none, 1 = Liu+2015, 2 = Wheeler+ 1975.

Parameter input type:

Float

Default value:

0

Macros:

[‘SN_IMPULSE_NONE = 0’, ‘SN_IMPULSE_LIU2015 = 1’, ‘SN_IMPULSE_WHEELER1975 = 2’]

sn_kick_companion_TZ

Description:

Set the speed (if >=0) of, or the algothim (if <0) used to calculate the, kick on the companion when a Thorne-Zytkow object is formed. 0 = none, 1 = Liu+2015, 2 = Wheeler+ 1975.

Parameter input type:

Float

Default value:

0

Macros:

[‘SN_IMPULSE_NONE = 0’, ‘SN_IMPULSE_LIU2015 = 1’, ‘SN_IMPULSE_WHEELER1975 = 2’]

sn_kick_companion_AIC_BH

Description:

Set the speed (if >=0) of, or the algothim (if <0) used to calculate the, kick on the companion when a neutron star collapses to a black hole. 0 = none, 1 = Liu+2015, 2 = Wheeler+ 1975.

Parameter input type:

Float

Default value:

0

Macros:

[‘SN_IMPULSE_NONE = 0’, ‘SN_IMPULSE_LIU2015 = 1’, ‘SN_IMPULSE_WHEELER1975 = 2’]

sn_kick_companion_BH_BH

Description:

Set the speed (if >=0) of, or the algothim (if <0) used to calculate the, kick on the companion when two black holes merge. 0 = none, 1 = Liu+2015, 2 = Wheeler+ 1975.

Parameter input type:

Float

Default value:

0

Macros:

[‘SN_IMPULSE_NONE = 0’, ‘SN_IMPULSE_LIU2015 = 1’, ‘SN_IMPULSE_WHEELER1975 = 2’]

sn_kick_companion_BH_NS

Description:

Set the speed (if >=0) of, or the algothim (if <0) used to calculate the, kick on the companion when a black hole merges with a neutron star. 0 = none, 1 = Liu+2015, 2 = Wheeler+ 1975.

Parameter input type:

Float

Default value:

0

Macros:

[‘SN_IMPULSE_NONE = 0’, ‘SN_IMPULSE_LIU2015 = 1’, ‘SN_IMPULSE_WHEELER1975 = 2’]

sn_kick_companion_IA_Hybrid_HeCOWD

Description:

Set the speed (if >=0) of, or the algothim (if <0) used to calculate the kick on the companion, if it survives, in a hybrid He-COWD type Ia explosion. 0 = none, 1 = Liu+2015, 2 = Wheeler+ 1975.

Parameter input type:

Float

Default value:

0

Macros:

[‘SN_IMPULSE_NONE = 0’, ‘SN_IMPULSE_LIU2015 = 1’, ‘SN_IMPULSE_WHEELER1975 = 2’]

sn_kick_companion_IA_Hybrid_HeCOWD_subluminous

Description:

Set the speed (if >=0) of, or the algothim (if <0) used to calculate the kick on the companion, if it survives, in a subluminous hybrid He-COWD type Ia explosion. 0 = none, 1 = Liu+2015, 2 = Wheeler+ 1975.

Parameter input type:

Float

Default value:

0

Macros:

[‘SN_IMPULSE_NONE = 0’, ‘SN_IMPULSE_LIU2015 = 1’, ‘SN_IMPULSE_WHEELER1975 = 2’]

wd_sigma

Description:

Set the speed at which white dwarfs are kicked when they form, in km/s. Default is zero (i.e. no kick). Requires WD_KICKS.

Parameter input type:

Float

Default value:

0

wd_kick_direction

Description:

Set the direction of white dwarf kicks. 0 = random, 1 = up, 2 = forward, 3 = backward, 4 = inward, 5 = outward. Requires WD_KICKS.

Parameter input type:

Integer

Default value:

0

Macros:

[‘KICK_RANDOM = 0’, ‘KICK_STRAIGHT_UP = 1’, ‘KICK_FORWARD = 2’, ‘KICK_BACKWARD = 3’, ‘KICK_INWARD = 4’, ‘KICK_OUTWARD = 5’]

wd_kick_when

Description:

Decide when to kick a white dwarf. 0=at birth, 1=at first RLOF, 2=at given pulse number (see wd_kick_pulse_number), 3 at every pulse Requires WD_KICKS.

Parameter input type:

Integer

Default value:

0

Macros:

[‘WD_KICK_END_AGB = 0’, ‘WD_KICK_FIRST_RLOF = 1’, ‘WD_KICK_AT_GIVEN_PULSE = 2’, ‘WD_KICK_AT_EVERY_PULSE = 3’]

wd_kick_pulse_number

Description:

Apply a kick to a star at a desired pulse number on the TPAGB (i.e. pre-WD). Requires WD_KICKS.

Parameter input type:

Integer

Default value:

0

minimum_helium_ignition_core_mass

Description:

Minimum helium core mass required to ignite helium in the case that the hydrogen envelope is stripped on the giant branch, e.g. to make an sdB or sdO star. Typically 0.4, if 0.0 then the BSE algorithm (based on the total mass) is used.

Parameter input type:

Float

Default value:

0

minimum_CO_core_mass_for_carbon_ignition

Description:

Minimum CO core mass for carbon ignition, assuming Mc,bagb>1.6Msun. Typically around 1.08Msun (Pols+1998).

Parameter input type:

Float

Default value:

1.08

minimum_CO_core_mass_for_neon_ignition

Description:

Minimum CO core mass for neon ignition. Typically around 1.42Msun. Stars that have cores that ignite carbon, but not neon explode in electon-capture supernovae.

Parameter input type:

Float

Default value:

1.42

minimum_mcbagb_for_nondegenerate_carbon_ignition

Description:

Minimum Mc,bagb (core mass at the base of the AGB) for non-degenerate carbon ignition. Typically around 2.25Msun (Pols+1998).

Parameter input type:

Float

Default value:

2.25

maximum_mcbagb_for_degenerate_carbon_ignition

Description:

Maximum Mc,bagb (core mass at the base of the AGB) for degenerate carbon ignition. Typically around 1.6Msun (Pols+1998).

Parameter input type:

Float

Default value:

1.6

max_neutron_star_mass

Description:

Maximum mass of a neutron star before it collapses to a black hole. Typically around 2Msun.

Parameter input type:

Float

Default value:

2.2

chandrasekhar_mass

Description:

The Chandrasekhar mass, usually 1.44Msun.

Parameter input type:

Float

Default value:

1.38

delta_mcmin

Description:

A parameter to reduce the minimum core mass for third dredge up to occur on the TPAGB. As used by Izzard and Tout (2004) to increase the amount of dredge up, hence carbon, in Magellanic cloud stars.

Parameter input type:

Float

Default value:

NULL

lambda_min

Description:

A parameter to increase the efficiency of third dredge up on the TPAGB. The efficiency is lambda * lambda_mult, and setting lambda_min>0 implies that, once Mc>Mcmin (see delta_mcmin) lambda=Max(lambda(fit to Karakas), lambda_min). As used by Izzard and Tout (2004) to increase the amount of dredge up, hence carbon, in Magellanic cloud stars. You can set this to THIRD_DREDGE_UP_LAMBDA_MIN_AUTO to match the results of Izzard and Tout (2004). See also lambda_multiplier.

Parameter input type:

Float

Default value:

NULL

Macros:

[‘THIRD_DREDGE_UP_LAMBDA_MIN_AUTO = -1’]

lambda_multiplier

Description:

A parameter to increase the efficiency of third dredge up on the TPAGB. The efficiency is lambda * lambda_mult, and setting lambda_min>0 implies that, once Mc>Mcmin (see delta_mcmin) lambda=Max(lambda(fit to Karakas), lambda_min). As used by Izzard and Tout (2004) to increase the amount of dredge up, hence carbon, in Magellanic cloud stars.

Parameter input type:

Float

Default value:

NULL

minimum_envelope_mass_for_third_dredgeup

Description:

The minimum envelope mass for third dredge up on the TPAGB. Early, solar metallicity models by Straniero et al suggested 0.5Msun is typical. However, circumstantial evidence (Izzard et al 2009) as well as newer models by Stancliffe and Karakas suggest that at low metallicity a value nearer zero is more appropriate.

Parameter input type:

Float

Default value:

NULL

mass_of_pmz

Description:

The mass in the partial mixing zone of a TPAGB star, using the Karakas 2012 tables. Ask Carlo Abate for more details, or see the series of papers Abate et al 2012, 2013, 2014. Requires NUCSYN and USE_TABULAR_INTERSHELL_ABUNDANCES_KARAKAS_2012.

Parameter input type:

Float

Default value:

NULL

c13_eff

Description:

The “efficiency” of partial mixing in a TPAGB star intershell region, when using the s-process tables of Gallino, Busso, Lugaro et al. as provided by Maria Lugaro for the Izzard et al. 2009 paper. Requires NUCSYN and NUCSYN_S_PROCESS.

Parameter input type:

Float

Default value:

NULL

mc13_pocket_multiplier

Description:

Multiplies the mass in the partial mixing zone of a TPAGB star, when using the s-process tables of Gallino, Busso, Lugaro et al. as provided by Maria Lugaro for the Izzard et al. 2009 paper. Requires NUCSYN and NUCSYN_S_PROCESS.

Parameter input type:

Float

Default value:

NULL

tides_convective_damping

Description:

Tidal convective damping algorithm. 0=TIDES_HURLEY2002 Zahn 197x timescales + Hut, as in Hurley et al (2002), 1 = TIDES_ZAHN1989: Zahn 1989 lambdas + Hut.

Parameter input type:

Integer

Default value:

0

Macros:

[‘TIDES_HURLEY2002 = 0’, ‘TIDES_ZAHN1989 = 1’, ‘TIDES_HURLEY2002 = 0’, ‘TIDES_ZAHN1989 = 1’]

E2_prescription

Description:

Choose how to calculate the E2 structural parameter (used in tidal timescale calculations). 0=Hurley 1=Izzard (see Siess et al 2013).

Parameter input type:

Integer

Default value:

0

Macros:

[‘E2_HURLEY_2002 = 0’, ‘E2_IZZARD = 1’, ‘E2_MINT = 2’, ‘Random variation : integer between 0 and 1 ‘]

dtfac

Description:

A parameter to decrease the timestep ONLY during the TPAGB phase.

Parameter input type:

Float

Default value:

1

hbbtfac

Description:

A parameter to modulate the temperature at the base of the hot-bottom burning zone in TPAGB stars. (Works only if NUCSYN is defined).

Parameter input type:

Float

Default value:

NULL

wind_multiplier_%d

Description:

Wind multiplier for the stellar type specified by the intger %d. By default these are all 1.0.

Parameter input type:

Float(scanf)

Default value:

NULL

wind_multiplier_LOW_MASS_MS

Description:

Wind multiplier for the stellar type specified by the intger -2104370523. By default these are all 1.0.

wind_multiplier_MS

Description:

Wind multiplier for the stellar type specified by the intger -2104376897. By default these are all 1.0.

wind_multiplier_HG

Description:

Wind multiplier for the stellar type specified by the intger -2104376864. By default these are all 1.0.

wind_multiplier_GIANT_BRANCH

Description:

Wind multiplier for the stellar type specified by the intger -2104370511. By default these are all 1.0.

wind_multiplier_CHeB

Description:

Wind multiplier for the stellar type specified by the intger -2104376857. By default these are all 1.0.

wind_multiplier_EAGB

Description:

Wind multiplier for the stellar type specified by the intger -2104376852. By default these are all 1.0.

wind_multiplier_TPAGB

Description:

Wind multiplier for the stellar type specified by the intger -2104376823. By default these are all 1.0.

wind_multiplier_HeMS

Description:

Wind multiplier for the stellar type specified by the intger -2104376722. By default these are all 1.0.

wind_multiplier_HeHG

Description:

Wind multiplier for the stellar type specified by the intger -2104370498. By default these are all 1.0.

wind_multiplier_HeGB

Description:

Wind multiplier for the stellar type specified by the intger -2104370493. By default these are all 1.0.

wind_multiplier_HeWD

Description:

Wind multiplier for the stellar type specified by the intger -2104370488. By default these are all 1.0.

wind_multiplier_COWD

Description:

Wind multiplier for the stellar type specified by the intger -2104370483. By default these are all 1.0.

wind_multiplier_ONeWD

Description:

Wind multiplier for the stellar type specified by the intger -2104370478. By default these are all 1.0.

wind_multiplier_NS

Description:

Wind multiplier for the stellar type specified by the intger -2104377658. By default these are all 1.0.

wind_multiplier_BH

Description:

Wind multiplier for the stellar type specified by the intger -2104370472. By default these are all 1.0.

wind_multiplier_MASSLESS_REMNANT

Description:

Wind multiplier for the stellar type specified by the intger -2104376290. By default these are all 1.0.

wind_type_multiplier_%d

Description:

Wind multiplier for different types of wind (MS, GB, AGB, WR, LBV, OTHER), given by the integer %d. By default these are all 1.0.

Parameter input type:

Float(scanf)

Default value:

NULL

wind_type_multiplier_LOW_MASS_MS

Description:

Wind multiplier for different types of wind (MS, GB, AGB, WR, LBV, OTHER), given by the integer -2104370523. By default these are all 1.0.

wind_type_multiplier_MS

Description:

Wind multiplier for different types of wind (MS, GB, AGB, WR, LBV, OTHER), given by the integer -2104376897. By default these are all 1.0.

wind_type_multiplier_HG

Description:

Wind multiplier for different types of wind (MS, GB, AGB, WR, LBV, OTHER), given by the integer -2104376864. By default these are all 1.0.

wind_type_multiplier_GIANT_BRANCH

Description:

Wind multiplier for different types of wind (MS, GB, AGB, WR, LBV, OTHER), given by the integer -2104370511. By default these are all 1.0.

wind_type_multiplier_CHeB

Description:

Wind multiplier for different types of wind (MS, GB, AGB, WR, LBV, OTHER), given by the integer -2104376857. By default these are all 1.0.

wind_type_multiplier_EAGB

Description:

Wind multiplier for different types of wind (MS, GB, AGB, WR, LBV, OTHER), given by the integer -2104376852. By default these are all 1.0.

wind_type_multiplier_TPAGB

Description:

Wind multiplier for different types of wind (MS, GB, AGB, WR, LBV, OTHER), given by the integer -2104376823. By default these are all 1.0.

wind_type_multiplier_HeMS

Description:

Wind multiplier for different types of wind (MS, GB, AGB, WR, LBV, OTHER), given by the integer -2104376722. By default these are all 1.0.

wind_type_multiplier_HeHG

Description:

Wind multiplier for different types of wind (MS, GB, AGB, WR, LBV, OTHER), given by the integer -2104370498. By default these are all 1.0.

wind_type_multiplier_HeGB

Description:

Wind multiplier for different types of wind (MS, GB, AGB, WR, LBV, OTHER), given by the integer -2104370493. By default these are all 1.0.

wind_type_multiplier_HeWD

Description:

Wind multiplier for different types of wind (MS, GB, AGB, WR, LBV, OTHER), given by the integer -2104370488. By default these are all 1.0.

wind_type_multiplier_COWD

Description:

Wind multiplier for different types of wind (MS, GB, AGB, WR, LBV, OTHER), given by the integer -2104370483. By default these are all 1.0.

wind_type_multiplier_ONeWD

Description:

Wind multiplier for different types of wind (MS, GB, AGB, WR, LBV, OTHER), given by the integer -2104370478. By default these are all 1.0.

wind_type_multiplier_NS

Description:

Wind multiplier for different types of wind (MS, GB, AGB, WR, LBV, OTHER), given by the integer -2104377658. By default these are all 1.0.

wind_type_multiplier_BH

Description:

Wind multiplier for different types of wind (MS, GB, AGB, WR, LBV, OTHER), given by the integer -2104370472. By default these are all 1.0.

wind_type_multiplier_MASSLESS_REMNANT

Description:

Wind multiplier for different types of wind (MS, GB, AGB, WR, LBV, OTHER), given by the integer -2104376290. By default these are all 1.0.

pre_main_sequence

Description:

Set to True to turn on pre-main sequence evolution. Currently this is not a special stellar type, rather the first (small) fraction of the main sequence has increased radii to match the Railton et al 2014 fits to Tout’s pre-main sequence stars. Requires PRE_MAIN_SEQUENCE. See also pre_main_sequence_fit_lobes.

Parameter input type:

True|False

Default value:

False

pre_main_sequence_fit_lobes

Description:

Set to True force a pre-main sequence star into its Roche lobe. This is done by artificially aging it. Requires PRE_MAIN_SEQUENCE.

Parameter input type:

True|False

Default value:

False

small_envelope_method

Description:

Choose the method used to determine the stellar radius when the envelope mass is very thin. 0 = Hurley et al. (2002), 1 = Miller Bertolami et al. (2016+) for GB and AGB stars only.

Parameter input type:

Integer

Default value:

0

Macros:

[‘SMALL_ENVELOPE_METHOD_BSE = 0’, ‘SMALL_ENVELOPE_METHOD_MILLER_BERTOLAMI = 1’]

timestep_modulator

Description:

Multiplier applied to the global timestep. Requires TIMESTEP_MODULATION.

Parameter input type:

Float

Default value:

1

timestep_multiplier%d

Description:

Multiplier applied to timestep limit <n>.

Parameter input type:

Float(scanf)

Default value:

NULL

timestep_multiplierNONE

Description:

Multiplier applied to timestep limit NONE.

timestep_multiplierMS

Description:

Multiplier applied to timestep limit MS.

timestep_multiplierPREMS

Description:

Multiplier applied to timestep limit PREMS.

timestep_multiplierPREROCHE_MS

Description:

Multiplier applied to timestep limit PREROCHE_MS.

timestep_multiplierPREROCHE_HG

Description:

Multiplier applied to timestep limit PREROCHE_HG.

timestep_multiplierHG

Description:

Multiplier applied to timestep limit HG.

timestep_multiplierFGB

Description:

Multiplier applied to timestep limit FGB.

timestep_multiplierCHeB

Description:

Multiplier applied to timestep limit CHeB.

timestep_multiplierEAGB

Description:

Multiplier applied to timestep limit EAGB.

timestep_multiplierEAGB_AXEL

Description:

Multiplier applied to timestep limit EAGB_AXEL.

timestep_multiplierEAGB_PREROCHE

Description:

Multiplier applied to timestep limit EAGB_PREROCHE.

timestep_multiplierTPAGB

Description:

Multiplier applied to timestep limit TPAGB.

timestep_multiplierTPAGB_NUCSYN_INTERPULSE

Description:

Multiplier applied to timestep limit TPAGB_NUCSYN_INTERPULSE.

timestep_multiplierTPAGB_NUCSYN_SPEEDUP

Description:

Multiplier applied to timestep limit TPAGB_NUCSYN_SPEEDUP.

timestep_multiplierTPAGB_NUCSYN_PREROCHE

Description:

Multiplier applied to timestep limit TPAGB_NUCSYN_PREROCHE.

timestep_multiplierTPAGB_NUCSYN_KARAKAS_SMOOTH

Description:

Multiplier applied to timestep limit TPAGB_NUCSYN_KARAKAS_SMOOTH.

timestep_multiplierHeMS

Description:

Multiplier applied to timestep limit HeMS.

timestep_multiplierHeHG_GB

Description:

Multiplier applied to timestep limit HeHG_GB.

timestep_multiplierOTHER_STELLAR_TYPES

Description:

Multiplier applied to timestep limit OTHER_STELLAR_TYPES.

timestep_multiplierSTELLAR_ANGMOM

Description:

Multiplier applied to timestep limit STELLAR_ANGMOM.

timestep_multiplierSTELLAR_MASS_LOSS

Description:

Multiplier applied to timestep limit STELLAR_MASS_LOSS.

timestep_multiplierSTELLAR_MAGNETIC_BRAKING

Description:

Multiplier applied to timestep limit STELLAR_MAGNETIC_BRAKING.

timestep_multiplierCIRCUMBINARY_DISC

Description:

Multiplier applied to timestep limit CIRCUMBINARY_DISC.

timestep_multiplier

Description:

Multiplier applied to timestep limit .

timestep_multiplierDISC

Description:

Multiplier applied to timestep limit DISC.

timestep_multiplierFASTWIND

Description:

Multiplier applied to timestep limit FASTWIND.

timestep_multiplierSELMA

Description:

Multiplier applied to timestep limit SELMA.

timestep_multiplierCEMP_POSTMS

Description:

Multiplier applied to timestep limit CEMP_POSTMS.

timestep_multiplierCEMP_NOTEMP

Description:

Multiplier applied to timestep limit CEMP_NOTEMP.

timestep_multiplierCEMP_EMP

Description:

Multiplier applied to timestep limit CEMP_EMP.

timestep_multiplierCEMP_NEARLY

Description:

Multiplier applied to timestep limit CEMP_NEARLY.

timestep_multiplierCEMP_FLOOR

Description:

Multiplier applied to timestep limit CEMP_FLOOR.

timestep_multiplierFABIAN_IMF_LOG

Description:

Multiplier applied to timestep limit FABIAN_IMF_LOG.

timestep_multiplierHRD1

Description:

Multiplier applied to timestep limit HRD1.

timestep_multiplierHRD2

Description:

Multiplier applied to timestep limit HRD2.

timestep_multiplier16

Description:

Multiplier applied to timestep limit 16.

timestep_multiplierYVT

Description:

Multiplier applied to timestep limit YVT.

timestep_multiplierMINIMUM_TIMESTEP

Description:

Multiplier applied to timestep limit MINIMUM_TIMESTEP.

timestep_multiplierMAXIMUM_TIMESTEP

Description:

Multiplier applied to timestep limit MAXIMUM_TIMESTEP.

timestep_multiplierMAXIMUM_TIMESTEP_BY_STELLAR_TYPE

Description:

Multiplier applied to timestep limit MAXIMUM_TIMESTEP_BY_STELLAR_TYPE.

timestep_multiplierNOVAE

Description:

Multiplier applied to timestep limit NOVAE.

timestep_multiplierARTIFICIAL_ACCRETION

Description:

Multiplier applied to timestep limit ARTIFICIAL_ACCRETION.

timestep_multiplierSN

Description:

Multiplier applied to timestep limit SN.

timestep_multiplierMASS_GAIN

Description:

Multiplier applied to timestep limit MASS_GAIN.

timestep_multiplierMASS_LOSS

Description:

Multiplier applied to timestep limit MASS_LOSS.

timestep_multiplierTIDES

Description:

Multiplier applied to timestep limit TIDES.

timestep_multiplierNUCSYN_ANGELOU_LITHIUM

Description:

Multiplier applied to timestep limit NUCSYN_ANGELOU_LITHIUM.

timestep_multiplierCARBON_BURNING

Description:

Multiplier applied to timestep limit CARBON_BURNING.

timestep_multiplierBURN_IN

Description:

Multiplier applied to timestep limit BURN_IN.

timestep_multiplierRADIUS_CHANGES

Description:

Multiplier applied to timestep limit RADIUS_CHANGES.

timestep_multiplierMASSLESS_REMNANT

Description:

Multiplier applied to timestep limit MASSLESS_REMNANT.

timestep_multiplierORBITAL_ANGMOM

Description:

Multiplier applied to timestep limit ORBITAL_ANGMOM.

timestep_multiplierSTELLAR_ANGMOM2

Description:

Multiplier applied to timestep limit STELLAR_ANGMOM2.

timestep_multiplierTIDES2

Description:

Multiplier applied to timestep limit TIDES2.

timestep_multiplierMASS_GAIN2

Description:

Multiplier applied to timestep limit MASS_GAIN2.

timestep_multiplierMASS_LOSS2

Description:

Multiplier applied to timestep limit MASS_LOSS2.

timestep_multiplierGRAVITATIONAL_WAVE_RADIATION

Description:

Multiplier applied to timestep limit GRAVITATIONAL_WAVE_RADIATION.

timestep_multiplierRLOF

Description:

Multiplier applied to timestep limit RLOF.

timestep_multiplierFIXED_TIMESTEP

Description:

Multiplier applied to timestep limit FIXED_TIMESTEP.

timestep_multiplierTIME_REMAINING

Description:

Multiplier applied to timestep limit TIME_REMAINING.

timestep_multiplierZOOMFAC

Description:

Multiplier applied to timestep limit ZOOMFAC.

timestep_multiplierTHERMAL

Description:

Multiplier applied to timestep limit THERMAL.

timestep_multiplierNUMBER

Description:

Multiplier applied to timestep limit NUMBER.

maximum_timestep

Description:

The maximum timestep (MYr). Must not be zero.

Parameter input type:

Float

Default value:

1e+20

maximum_timestep_by_stellar_type_%d

Description:

The maximum timestep by stellar type (MYr). Ignored if zero (default 0).

Parameter input type:

Float(scanf)

Default value:

NULL

maximum_timestep_by_stellar_type_LOW_MASS_MS

Description:

The maximum timestep by stellar type (MYr). Ignored if zero (default 0).

maximum_timestep_by_stellar_type_MS

Description:

The maximum timestep by stellar type (MYr). Ignored if zero (default 0).

maximum_timestep_by_stellar_type_HG

Description:

The maximum timestep by stellar type (MYr). Ignored if zero (default 0).

maximum_timestep_by_stellar_type_GIANT_BRANCH

Description:

The maximum timestep by stellar type (MYr). Ignored if zero (default 0).

maximum_timestep_by_stellar_type_CHeB

Description:

The maximum timestep by stellar type (MYr). Ignored if zero (default 0).

maximum_timestep_by_stellar_type_EAGB

Description:

The maximum timestep by stellar type (MYr). Ignored if zero (default 0).

maximum_timestep_by_stellar_type_TPAGB

Description:

The maximum timestep by stellar type (MYr). Ignored if zero (default 0).

maximum_timestep_by_stellar_type_HeMS

Description:

The maximum timestep by stellar type (MYr). Ignored if zero (default 0).

maximum_timestep_by_stellar_type_HeHG

Description:

The maximum timestep by stellar type (MYr). Ignored if zero (default 0).

maximum_timestep_by_stellar_type_HeGB

Description:

The maximum timestep by stellar type (MYr). Ignored if zero (default 0).

maximum_timestep_by_stellar_type_HeWD

Description:

The maximum timestep by stellar type (MYr). Ignored if zero (default 0).

maximum_timestep_by_stellar_type_COWD

Description:

The maximum timestep by stellar type (MYr). Ignored if zero (default 0).

maximum_timestep_by_stellar_type_ONeWD

Description:

The maximum timestep by stellar type (MYr). Ignored if zero (default 0).

maximum_timestep_by_stellar_type_NS

Description:

The maximum timestep by stellar type (MYr). Ignored if zero (default 0).

maximum_timestep_by_stellar_type_BH

Description:

The maximum timestep by stellar type (MYr). Ignored if zero (default 0).

maximum_timestep_by_stellar_type_MASSLESS_REMNANT

Description:

The maximum timestep by stellar type (MYr). Ignored if zero (default 0).

zoomfac_multiplier_decrease

Description:

When a timestep is rejected, decrease the timestep by this factor (0.5).

Parameter input type:

Float

Default value:

0.5

zoomfac_multiplier_increase

Description:

When a timestep is rejected, zooms, then succeeds, increase the timestep by this factor (1.2).

Parameter input type:

Float

Default value:

1.2

maximum_timestep_factor

Description:

The maximum factor between two subsequent timesteps (1.2).

Parameter input type:

Float

Default value:

0

maximum_nuclear_burning_timestep

Description:

The maximum timestep (MYr) in any nuclear burning phase.

Parameter input type:

Float

Default value:

1e+20

nova_retention_method

Description:

Algorithm used to calculate the amount of mass retained during a nova explosion. 0=use nova_retention_fraction. (other methods pending).

Parameter input type:

Integer

Default value:

0

Macros:

[‘NOVA_RETENTION_ALGORITHM_CONSTANT = 0’, ‘NOVA_RETENTION_ALGORITHM_CLAEYS2014 = 1’, ‘NOVA_RETENTION_ALGORITHM_HILLMAN2015 = 2’, ‘NOVA_RETENTION_ALGORITHM_WANGWU = 3’]

nova_yield_CO_algorithm

Description:

Algorithm used to calculate the yields of COWD novae. Current options: NOVA_YIELD_CO_ALGOTHIM_JOSE_HERNANZ_1998 (0).

Parameter input type:

Integer

Default value:

NULL

nova_yield_ONe_algorithm

Description:

Algorithm used to calculate the yields of ONeWD novae. Current options: NOVA_YIELD_ONe_ALGOTHIM_JOSE_HERNANZ_1998 (0), NOVA_YIELD_ONe_ALGORITHM_JOSE2022 (1).

Parameter input type:

Integer

Default value:

NULL

MINT_metallicity

Description:

This sets the metallicity for MINT. It is ignored if set to -1.0, the default, in which case the normal metallicity parameter is used.

Parameter input type:

Float

Default value:

-1

gaia_Teff_binwidth

Description:

Log10(Effective temperature) bin width used to make Gaia-like HRDs.

Parameter input type:

Float

Default value:

NULL

gaia_L_binwidth

Description:

Log10(luminosity) bin width used to make Gaia-like HRDs.

Parameter input type:

Float

Default value:

NULL

gaia_colour_transform_method

Description:

Use this to select the method used to transform to Gaia colours from other colour schemes. GAIA_CONVERSION_UBVRI_UNIVARIATE_JORDI2010 = 0 Jordi et al.’s univariate UBVRI fits, GAIA_CONVERSION_UBVRI_BIVARIATE_JORDI2010 = 1 Jordi et al.’s bivariate UBVRI fits, GAIA_CONVERSION_ugriz_UNIVARIATE_JORDI2010 = 2 Jordi et al.’s univariate UBVRI fits, GAIA_CONVERSION_ugriz_BIVARIATE_JORDI2010 = 3 Jordi et al.’s univariate ugriv fits, GAIA_CONVERSION_UBVRI_UNIVARIATE_EVANS2018 = 4 Evans et al. (2018, DR2) fits, GAIA_CONVERSION_ugriz_UNIVARIATE_EVANS2018 = 5 Evans et al. (2018, DR2) fits, GAIA_CONVERSION_UBVRI_RIELLO2020 = 6 Riello et al. (2020, DR3) fits, GAIA_CONVERSION_ugriz_RIELLO2020 = 7 Riello et al. (2020, DR3) fits.

Parameter input type:

Integer

Default value:

4

Macros:

[‘GAIA_CONVERSION_UBVRI_UNIVARIATE_JORDI2010 = 0’, ‘GAIA_CONVERSION_UBVRI_BIVARIATE_JORDI2010 = 1’, ‘GAIA_CONVERSION_ugriz_UNIVARIATE_JORDI2010 = 2’, ‘GAIA_CONVERSION_ugriz_BIVARIATE_JORDI2010 = 3’, ‘GAIA_CONVERSION_UBVRI_UNIVARIATE_EVANS2018 = 4’, ‘GAIA_CONVERSION_ugriz_UNIVARIATE_EVANS2018 = 5’, ‘GAIA_CONVERSION_UBVRI_RIELLO2020 = 6’, ‘GAIA_CONVERSION_ugriz_RIELLO2020 = 7’]

rotationally_enhanced_mass_loss

Description:

Set to 1 to enable rotationally enhanced mass loss rate algorithms: 0= none, 1=formula cf. Langer models (=ROTATIONALLY_ENHANCED_MASSLOSS_LANGER_FORMULA), 2=limit accretion rate before wind loss is applied, 3 = both 1 and 2. See also rotationally_enhanced_exponent.

Parameter input type:

Integer

Default value:

0

Macros:

[‘ROTATIONALLY_ENHANCED_MASSLOSS_NONE = 0’, ‘ROTATIONALLY_ENHANCED_MASSLOSS_LANGER_FORMULA = 1’, ‘ROTATIONALLY_ENHANCED_MASSLOSS_ANGMOM = 2’, ‘ROTATIONALLY_ENHANCED_MASSLOSS_LANGER_FORMULA_AND_ANGMOM = 3’]

AGB_core_algorithm

Description:

Algorithm to use for calculating AGB core masses. 0=Hurley et al. 2002 if no NUCSYN, Karakas 2002 if NUCSYN is defined; 1=Hurley et al. 2002 (overshooting models); 1=Karakas 2002 (non-overshooting models).

Parameter input type:

Integer

Default value:

1

Macros:

[‘AGB_CORE_ALGORITHM_DEFAULT = 0’, ‘AGB_CORE_ALGORITHM_HURLEY = 1’, ‘AGB_CORE_ALGORITHM_KARAKAS = 2’]

AGB_radius_algorithm

Description:

Algorithm to use for calculating radii on the TPAGB.

Parameter input type:

Integer

Default value:

1

Macros:

[‘AGB_RADIUS_ALGORITHM_DEFAULT = 0’, ‘AGB_RADIUS_ALGORITHM_HURLEY = 1’, ‘AGB_RADIUS_ALGORITHM_KARAKAS = 2’]

AGB_luminosity_algorithm

Description:

Algorithm to use for calculating luminosities on the TPAGB.

Parameter input type:

Integer

Default value:

1

Macros:

[‘AGB_LUMINOSITY_ALGORITHM_DEFAULT = 0’, ‘AGB_LUMINOSITY_ALGORITHM_HURLEY = 1’, ‘AGB_LUMINOSITY_ALGORITHM_KARAKAS = 2’]

AGB_3dup_algorithm

Description:

Algorithm to use for calculating third dredge up efficiency on the TPAGB.

Parameter input type:

Integer

Default value:

1

Macros:

[‘AGB_THIRD_DREDGE_UP_ALGORITHM_DEFAULT = 0’, ‘AGB_THIRD_DREDGE_UP_ALGORITHM_HURLEY = 1’, ‘AGB_THIRD_DREDGE_UP_ALGORITHM_KARAKAS = 2’, ‘AGB_THIRD_DREDGE_UP_ALGORITHM_STANCLIFFE = 3’]

overspin_algorithm

Description:

Algorithm to determine what we do when a star is rotating at its breakup velocity. OVERSPIN_BSE (0) conservatively transfers the angular momentum back to the orbit. OVERSPIN_MASSLOSS uses the rotationally_enhanced_massloss parameter to lose mass which carries away the angular momentum.

Parameter input type:

Integer

Default value:

0

Macros:

[‘OVERSPIN_BSE = 0’, ‘OVERSPIN_MASSLOSS = 1’]

rotationally_enhanced_exponent

Description:

The exponent (power) by which rotationally enhanced mass loss is raised. Requires ROTATIONALLY_ENHANCED_MASS_LOSS. See also rotationally_enhanced_mass_loss.

Parameter input type:

Float

Default value:

1

batchmode

Description:

Set the batchmode control variable. Use only if you know what you are doing!.

Parameter input type:

Integer

Default value:

3

speedtests

Description:

If True, turns on speedtests during version information (off by default).

Parameter input type:

True|False

Default value:

False

use_fixed_timestep_%d

Description:

Set to True to use fixed timestep <n>, False to turn off. Fixed timesteps are on (this is True) by default.

Parameter input type:

Boolean(scanf)

Default value:

NULL

use_fixed_timestep_ENSEMBLE

Description:

Set to True to use fixed timestep ENSEMBLE, False to turn off. Fixed timesteps are on (this is True) by default.

use_fixed_timestep_TEST

Description:

Set to True to use fixed timestep TEST, False to turn off. Fixed timesteps are on (this is True) by default.

use_fixed_timestep_NUMBER

Description:

Set to True to use fixed timestep NUMBER, False to turn off. Fixed timesteps are on (this is True) by default.

task%d

Description:

Control tasks to be performed by binary_c. By default, these are all TRUE. For more information see binary_c_macros.h, particularly the BINARY_C_TASK_* macros.

Parameter input type:

Boolean(scanf)

Default value:

NULL

taskTIME_EVOLUTION

Description:

Control tasks to be performed by binary_c. By default, these are all TRUE. For more information see binary_c_macros.h, particularly the BINARY_C_TASK_* macros.

taskCALCULATE_DERIVATIVES

Description:

Control tasks to be performed by binary_c. By default, these are all TRUE. For more information see binary_c_macros.h, particularly the BINARY_C_TASK_* macros.

taskAPPLY_DERIVATIVES

Description:

Control tasks to be performed by binary_c. By default, these are all TRUE. For more information see binary_c_macros.h, particularly the BINARY_C_TASK_* macros.

taskCALCULATE_STELLAR_EVOLUTION

Description:

Control tasks to be performed by binary_c. By default, these are all TRUE. For more information see binary_c_macros.h, particularly the BINARY_C_TASK_* macros.

taskCALCULATE_BINARY_EVOLUTION

Description:

Control tasks to be performed by binary_c. By default, these are all TRUE. For more information see binary_c_macros.h, particularly the BINARY_C_TASK_* macros.

taskNUMBER

Description:

Control tasks to be performed by binary_c. By default, these are all TRUE. For more information see binary_c_macros.h, particularly the BINARY_C_TASK_* macros.

orbiting_object

Description:

Usage: –orbiting_object mass,spinrate,central_object,period.

Parameter input type:
Default value:

NULL

Extra:

1.0

orbiting_objects_log

Description:

If True, turn on orbiting-object log.

Parameter input type:

True|False

Default value:

False

orbiting_objects_wind_accretion_multiplier

Description:

Multiplier for wind accretion on orbiting objects. Hurley et al 2002 use 1.5, which is the default.

Parameter input type:

Float

Default value:

1.5

orbiting_objects_close_pc_threshold

Description:

How close are orbiting objects allowed to be? Set this to be the absolute percentage difference minimum.

Parameter input type:

Float

Default value:

1

orbiting_objects_tides_multiplier

Description:

Multiplier for tidal torques on orbiting objects.

Parameter input type:

Float

Default value:

1

evaporate_escaped_orbiting_objects

Description:

If True, evaporate orbiting objects that have escaped the system.

Parameter input type:

True|False

Default value:

False

RLOF_transition_objects_escape

Description:

If True, objects that escape their Roche lobe are ejected from the system, otherwise they are placed just outside the minimum stable orbit.

Parameter input type:

True|False

Default value:

False

PN_resolve

Description:

If True, the timestep will be shortened to resolve better the PN phase (FALSE).

Parameter input type:

True|False

Default value:

False

PN_resolve_minimum_luminosity

Description:

The luminosity above which extra time resolution for PNe is applied (see PN_resolve).

Parameter input type:

Float

Default value:

31.62

PN_resolve_maximum_envelope_mass

Description:

The envelope mass below which extra time resolution for PNe is applied (see PN_resolve).

Parameter input type:

Float

Default value:

0.1

PN_resolve_minimum_effective_temperature

Description:

The minimum effective temperature above which extra time resolution for PNe is applied (see PN_resolve).

Parameter input type:

Float

Default value:

12500

PN_fast_wind

Description:

If True, thin-envelope PNe will have fast winds (FALSE).

Parameter input type:

True|False

Default value:

False

PN_fast_wind_dm_GB

Description:

The envelope mass below which fast wind used during the GB if PN_fast_wind is TRUE. (See also PN_fast_wind, PN_fast_wind_mdot_GB).

Parameter input type:

Float

Default value:

0.01

PN_fast_wind_mdot_GB

Description:

The envelope mass below which fast wind used during the GB if PN_fast_wind is TRUE. (See also PN_fast_wind, PN_fast_wind_mdot_GB).

Parameter input type:

Float

Default value:

1e-06

PN_fast_wind_dm_AGB

Description:

The envelope mass below which fast wind used during the AGB if PN_fast_wind is TRUE. (See also PN_fast_wind, PN_fast_wind_mdot_AGB).

Parameter input type:

Float

Default value:

0.001

PN_fast_wind_mdot_AGB

Description:

The envelope mass below which fast wind used during the GB if PN_fast_wind is TRUE. (See also PN_fast_wind, PN_fast_wind_mdot_AGB).

Parameter input type:

Float

Default value:

1e-06

HeWD_HeWD_ignition_mass

Description:

HeWD-HeWD mergers above this mass reignite helium. (0.3).

Parameter input type:

Float

Default value:

0.3

wind_Nieuwenhuijzen_luminosity_lower_limit

Description:

Above this luminosity we activate the Nieuwenhuijzen and de Jager wind (4e3 Lsun).

Parameter input type:

Float

Default value:

4000

wind_LBV_luminosity_lower_limit

Description:

Above this luminosity we activate the LBV wind (6e5 Lsun).

Parameter input type:

Float

Default value:

600000

colour%d

Description:

Sets colour %d (0 to NUM_ANSI_COLOURS-1) to the extended ANSI set colour you choose (1-255, 0 means ignore). The colour numbers are defined in src/logging/ansi_colours.h.

Parameter input type:

Int(scanf)

Default value:

NULL

apply_Darwin_Radau_correction

Description:

Apply Darwin-Radau correction to the moment of inertia to take rotation into account?.

Parameter input type:

True|False

Default value:

False

adjust_structure_from_mass_changes

Description:

If True, adjust L and R based on changes in M on the previous thermal timescale. (Currently experimental hence defaults to FALSE).

Parameter input type:

True|False

Default value:

NULL

save_mass_history_n_thermal

Description:

Number of thermal timesteps of mass change history to be saved (4).

Parameter input type:

Float

Default value:

NULL

degenerate_core_merger_nucsyn

Description:

If TRUE, assume that in a degnerate core merger, energy is generated from nucleosynthesis of the whole core, and that this can disrupt the core. The BSE algorithm (Hurley et al. 2002) assumes this to be TRUE, but binary_c assumes FALSE by default. (FALSE).

Parameter input type:

True|False

Default value:

False

degenerate_core_helium_merger_ignition

Description:

If TRUE, assume that when there is a degenerate helium core merger, the star reignites helium. This is required to make R-type carbon stars. (TRUE).

Parameter input type:

True|False

Default value:

True

degenerate_core_merger_dredgeup_fraction

Description:

If non-zero, mix this fraction of the degenerate core during a merger.(0.0).

Parameter input type:

Float

Default value:

NULL

PPISN_prescription

Description:

(Pulsational) Pair-Instability Supernova prescription: handles remnant mass calculation for stars that undergo PPISN or PISN. Requires PPISN to be defined.

Parameter input type:

Integer

Default value:

0

Macros:

[‘PPISN_NONE = 0’, ‘PPISN_BELCZYNSKI16 = 1’, ‘PPISN_SPERA17 = 2’, ‘PPISN_STEVENSON19 = 3’, ‘PPISN_FARMER19 = 4’, ‘PPISN_MAPELLI20 = 5’, ‘PPISN_HENDRIKS23 = 6’]

Extra:

Ignore

PPISN_additional_massloss

Description:

(Pulsational) Pair-Instability Supernova additional mass loss on top of that given by the chosen prescription. Requires PPISN to be defined.

Parameter input type:

Float

Default value:

0

Extra:

Ignore

PPISN_core_mass_range_shift

Description:

(Pulsational) Pair-Instability Supernova CO core mass range shift: Mass by which we shift the range of CO core mass that undergoes PPISN. This is taken into account in the fit that calculates mass removal. Negative shifts the range to lower masses, positive to higher masses. Requires PPISN to be defined.

Parameter input type:

Float

Default value:

0

Extra:

Ignore

PPISN_massloss_multiplier

Description:

(Pulsational) Pair-Instability Supernova CO mass loss multiplier. Currently only works with the new prescription. Does not change additional mass loss. Requires PPISN to be defined.

Parameter input type:

Float

Default value:

1

Extra:

Ignore

fixed_beta_mass_transfer_efficiency

Description:

Parameter to fix mass transfer efficiency. Choosing a value between 0 and 1 enables this functionality and disables the rate limiters accretion_limit_thermal_multiplier and accretion_limit_dynamical_multiplier, and sets the accretion efficiency to the value provided or that of the Eddington accretion limit (whichever has the most excess). Default value is FIXED_BETA_MASS_TRANSFER_EFFICIENCY_DISABLED == -1.0, which disables this functionality.

Parameter input type:

Float

Default value:

-1

Macros:

[‘FIXED_BETA_MASS_TRANSFER_EFFICIENCY_DISABLED = -1’]

Extra:

Ignore

Section: binary

binary

Option

Description

separation

Description:

Set the orbital separation (actually the semi-major axis) of the binary (internal index 0, stellar indices 0 and 1) in solar radii. Note that if the orbital period is given, it is used to calculate the separation and the value you set with this argument is ignored. So if you want to set the separation instead, either do not set the orbital period or set the orbital period to zero (0.0).

Parameter input type:

Float

Default value:

0

separation_triple

Description:

Set the orbital separation (actually the semi-major axis) of the triple (internal index 1) in solar radii. Note that if the orbital period is given, it is used to calculate the separation. So if you want to set the separation instead, either do not set the orbital period or set the orbital period to zero (0.0).

Parameter input type:

Float

Default value:

0

separation_quadruple

Description:

Set the orbital separation (actually the semi-major axis) of the quadruple (internal index 2) in solar radii. Note that if the orbital period is given, it is used to calculate the separation. So if you want to set the separation instead, either do not set the orbital period or set the orbital period to zero (0.0).

Parameter input type:

Float

Default value:

0

orbital_period

Description:

Set the initial orbital period of the binary, stars 1 and 2 (internal indices 0 and 1) in days. See also separation.

Parameter input type:

Float

Default value:

0

orbital_period_triple

Description:

Set the initial orbital period of the triple in days. See also separation.

Parameter input type:

Float

Default value:

0

orbital_period_quadruple

Description:

Set the orbital period of the outer binary in a quadrulple (stars 3 and 4, internal indices 2 and 3) in days. See also separation.

Parameter input type:

Float

Default value:

0

eccentricity

Description:

Set the initial eccentricity of the binary orbit (stars 1 and 2, internal indices 0 and 1).

Parameter input type:

Float

Default value:

0

eccentricity_triple

Description:

Set the initial eccentricity of the triple orbit.

Parameter input type:

Float

Default value:

0

eccentricity_quadruple

Description:

Set the initial eccentricity of the outer binary of a quadruple (stars 3 and 4, internal indices 2 and 3).

Parameter input type:

Float

Default value:

0

incliniation

Description:

Set the initial orbital_inclination of the binary relative to zero.

Parameter input type:

Float

Default value:

0

incliniation_triple

Description:

Set the initial orbital_inclination of the triple orbit relative to zero.

Parameter input type:

Float

Default value:

0

incliniation_quadruple

Description:

Set the initial orbital_inclinationy of the quadruple orbit relative to zero.

Parameter input type:

Float

Default value:

0

orbital_phase

Description:

Set the initial orbital phase of the binary orbit.

Parameter input type:

Float

Default value:

0

orbital_phase_triple

Description:

Set the initial orbital phase of the triple orbit.

Parameter input type:

Float

Default value:

0

orbital_phase_quadruple

Description:

Set the initial orbital phase of the quadruple orbit.

Parameter input type:

Float

Default value:

0

argument_of_periastron

Description:

Set the initial argument of periastron of the binary orbit.

Parameter input type:

Float

Default value:

0

argument_of_periastron_triple

Description:

Set the initial argument of periastron of the triple orbit.

Parameter input type:

Float

Default value:

0

argument_of_periastron_quadruple

Description:

Set the initial argument of periastron of the quadruple orbit.

Parameter input type:

Float

Default value:

0

disc_timestep_factor

Description:

Factor that multiplies the natural timestep of a disc.

Parameter input type:

Float

Default value:

0.01

white_dwarf_cooling_model

Description:

White dwarf cooling model, relates age to luminosity. WHITE_DWARF_COOLING_MESTEL = 0 is Mestel’s model, WHITE_DWARF_COOLING_MESTEL_MODIFIED = 1 is Hurley’s modified Mestel model, WHITE_DWARF_COOLING_CARRASCO2014 = 2 is based on Carrasco (2014) tables.

Parameter input type:

Integer

Default value:

0

Macros:

[‘WHITE_DWARF_COOLING_MESTEL = 0’, ‘WHITE_DWARF_COOLING_MESTEL_MODIFIED = 1’, ‘WHITE_DWARF_COOLING_CARRASCO2014 = 2’]

white_dwarf_radius_model

Description:

White dwarf radius model, radius to mass (and perhaps age). WHITE_DWARF_RADIUS_NAUENBERG1972 = 0 Nauenberg (1972), WHITE_DWARF_RADIUS_MU = 1 mu-dependent variant, WHITE_DWARF_RADIUS_CARRASCO2014 = 2 is based on Carrasco (2014) tables.

Parameter input type:

Integer

Default value:

0

Macros:

[‘WHITE_DWARF_RADIUS_NAUENBERG1972 = 0’, ‘WHITE_DWARF_RADIUS_MU = 1’, ‘WHITE_DWARF_RADIUS_CARRASCO2014 = 2’, ‘Random variation : integer between 0 and 2 ‘]

cbdisc_mass_loss_inner_viscous_accretion_method

Description:

Chooses where the mass that is accreted from the inner edge of a circumbinary disc goes, i.e. to which star. 0 = Young and Clarke 2015, 1 = Gerosa et al 2015, 2 = 50:50 (i.e. not dependence on mass).

Parameter input type:

Integer

Default value:

0

Macros:

[‘CBDISC_MASS_LOSS_INNER_VISCOUS_ACCRETION_METHOD_YOUNG_CLARKE_2015 = 0’, ‘CBDISC_MASS_LOSS_INNER_VISCOUS_ACCRETION_METHOD_GEROSA_2015 = 1’, ‘CBDISC_MASS_LOSS_INNER_VISCOUS_ACCRETION_METHOD_EQUAL = 2’, ‘CBDISC_MASS_LOSS_INNER_VISCOUS_ACCRETION_METHOD_NONE = 3’]

cbdisc_inner_edge_stripping

Description:

If True, allow inner edge mass stripping.

Parameter input type:

True|False

Default value:

True

cbdisc_end_evolution_after_disc

Description:

If True, stop evolution when a disc evaporates.

Parameter input type:

True|False

Default value:

False

cbdisc_no_wind_if_cbdisc

Description:

If True, disable stellar winds when there is a circumbinary disc.

Parameter input type:

True|False

Default value:

False

cbdisc_outer_edge_stripping

Description:

If True, allow outer edge mass stripping.

Parameter input type:

True|False

Default value:

True

disc_n_monte_carlo_guesses

Description:

Number of monte carlo guesses to try in the disc solver if the normal list of guesses fails (0).

Parameter input type:

Integer

Default value:

0

disc_log

Description:

If 1, turn on the disc log. Requires DISC_LOG to be defined on build.

Parameter input type:

Integer

Default value:

0

Macros:

[‘DISC_LOG_LEVEL_NONE = 0’, ‘DISC_LOG_LEVEL_NORMAL = 1’, ‘DISC_LOG_LEVEL_SUBTIMESTEP = 2’, ‘DISC_LOG_LEVEL_NORMAL_FIRST_DISC_ONLY = -1’, ‘DISC_LOG_LEVEL_SUBTIMESTEP_FIRST_DISC_ONLY = -2’]

disc_log2d

Description:

If 1, turn on the 2d disc log. Requires DISC_LOG to be defined on build.

Parameter input type:

Integer

Default value:

0

Macros:

[‘DISC_LOG_LEVEL_NONE = 0’, ‘DISC_LOG_LEVEL_NORMAL = 1’, ‘DISC_LOG_LEVEL_SUBTIMESTEP = 2’, ‘DISC_LOG_LEVEL_NORMAL_FIRST_DISC_ONLY = -1’, ‘DISC_LOG_LEVEL_SUBTIMESTEP_FIRST_DISC_ONLY = -2’]

disc_log_dt

Description:

If non-zero, only allows disc log output every disc_log_dt Myr.

Parameter input type:

Float

Default value:

0

disc_log_directory

Description:

Directory into which disc logging is sent (must exist!).

Parameter input type:

String

Default value:

/tmp/

Extra:

/tmp/

post_ce_adaptive_menv

Description:

If TRUE, and if post_ce_objects_have_envelopes is TRUE, then the envelope mass of a post-CE star is such that it sits just inside its Roche lobe. If FALSE then a fixed (thin) envelope mass is applied that depends on the stellar type (see macros POST_CE_ENVELOPE_DM_GB, POST_CE_ENVELOPE_DM_EAGB and POST_CE_ENVELOPE_DM_TPAGB).

Parameter input type:

True|False

Default value:

False

post_ce_objects_have_envelopes

Description:

If TRUE then post-common-envelope objects have thin envelopes. You need this if you are to have post-CE post-AGB stars. Note that this may be unstable, i.e. you may end up having many CEEs. The mass in the envelope is controlled by post_ce_adaptive_menv. TRUE by default.

Parameter input type:

True|False

Default value:

False

PN_comenv_transition_time

Description:

Post-common envelope transition time in years (1e2). This is the time taken to move from CEE ejection to Teff > 30e4 K. Hall et al. (2013) suggest ~100 years.

Parameter input type:

Float

Default value:

100

minimum_time_between_PNe

Description:

The minimum time (Myr) between planetary nebula detections. This prevents multiple, fast common envelopes triggering two PNe (0.1).

Parameter input type:

Float

Default value:

0.1

PN_Hall_fading_time_algorithm

Description:

In stars with low mass (<0.45Msun) cores, you can choose to set the PN fading time to either the minimum (PN_HALL_FADING_TIME_ALGORITHM_MINIMUM) or maximum (PN_HALL_FADING_TIME_ALGORITHM_MAXIMUM) as shown in Fig. 6 of Hall et al. (2013).

Parameter input type:

Integer

Default value:

0

Macros:

[‘PN_HALL_FADING_TIME_ALGORITHM_MINIMUM = 0’, ‘PN_HALL_FADING_TIME_ALGORITHM_MAXIMUM = 1’]

PPN_envelope_mass

Description:

Desired pre-planetary nebula (post-AGB) envelope mass.

Parameter input type:

Float

Default value:

0.01

cbdisc_eccentricity_pumping_method

Description:

Select from various eccentricity-pumping methods when there is a circumbinary disc. Requires DISCS. 0 = off.

Parameter input type:

Integer

Default value:

1

Macros:

[‘CBDISC_ECCENTRICITY_PUMPING_NONE = 0’, ‘CBDISC_ECCENTRICITY_PUMPING_DERMINE = 1’]

cbdisc_viscous_photoevaporative_coupling

Description:

Set to 1 to turn on viscous-photoevaporative coupling in circumbinary discs. Requires DISCS. 0 = CBDISC_VISCOUS_PHOTOEVAPORATIVE_COUPLING_NONE = off, 1 = CBDISC_VISCOUS_PHOTOEVAPORATIVE_COUPLING_INSTANT instant, 2 = CBDISC_VISCOUS_PHOTOEVAPORATIVE_COUPLING_VISCOUS slow, viscous wind.

Parameter input type:

Integer

Default value:

1

Macros:

[‘CBDISC_VISCOUS_PHOTOEVAPORATIVE_COUPLING_NONE = 0’, ‘CBDISC_VISCOUS_PHOTOEVAPORATIVE_COUPLING_INSTANT = 1’, ‘CBDISC_VISCOUS_PHOTOEVAPORATIVE_COUPLING_VISCOUS = 2’]

cbdisc_inner_edge_stripping_timescale

Description:

Defines the timescale for mass loss from by inner edge stripping. 0 = instant, 1 = very long, 2 = viscous at Revap_in, 3 = orbital at Revap_in.

Parameter input type:

Integer

Default value:

1

Macros:

[‘DISC_STRIPPING_TIMESCALE_INSTANT = 1’, ‘DISC_STRIPPING_TIMESCALE_INFINITE = 2’, ‘DISC_STRIPPING_TIMESCALE_VISCOUS = 3’, ‘DISC_STRIPPING_TIMESCALE_ORBIT = 4’]

cbdisc_outer_edge_stripping_timescale

Description:

Defines the timescale for mass loss from by outer edge stripping. 0 = instant, 1 = very long, 2 = viscous at Revap_in, 3 = orbital at Revap_out.

Parameter input type:

Integer

Default value:

1

Macros:

[‘DISC_STRIPPING_TIMESCALE_INSTANT = 1’, ‘DISC_STRIPPING_TIMESCALE_INFINITE = 2’, ‘DISC_STRIPPING_TIMESCALE_VISCOUS = 3’, ‘DISC_STRIPPING_TIMESCALE_ORBIT = 4’]

cbdisc_viscous_L2_coupling

Description:

Set to 1 to turn on viscous-L2-loss coupling in circumbinary discs. Requires DISCS. 0 = off.

Parameter input type:

Integer

Default value:

1

gravitational_radiation_model

Description:

Model for gravitational radiation from the system. 0=Hurley et al 2002 (Landau and Lifshitz 1951). 1 = as 0 but only when there is no RLOF. 2 = none.

Parameter input type:

Integer

Default value:

0

Macros:

[‘GRAVITATIONAL_RADIATION_BSE = 0’, ‘GRAVITATIONAL_RADIATION_BSE_WHEN_NO_RLOF = 1’, ‘GRAVITATIONAL_RADIATION_NONE = 2’, ‘GRAVITATIONAL_RADIATION_LANDAU_LIFSHITZ = 3’, ‘GRAVITATIONAL_RADIATION_LANDAU_LIFSHITZ_WHEN_NO_RLOF = 4’, ‘Random variation : integer between 0 and 4 ‘]

nova_irradiation_multiplier

Description:

Multiplier for nova-radiative induced mass loss. (Shara+1986).

Parameter input type:

Float

Default value:

0

gravitational_radiation_modulator_J

Description:

Modulator for gravitational wave radiation angular momentum loss rate (1.0).

Parameter input type:

Float

Default value:

1

gravitational_radiation_modulator_e

Description:

Modulator for gravitational wave radiation eccentricity pumping rate (1.0).

Parameter input type:

Float

Default value:

1

nova_faml_multiplier

Description:

Nova friction-induced angular momentum loss multiplier. (Shara+1986).

Parameter input type:

Float

Default value:

0

RLOF_angular_momentum_transfer_model

Description:

Choose angular momentum transfer model in RLOF. 0=BSE (with discs), 1=conservative.

Parameter input type:

Integer

Default value:

0

Macros:

[‘RLOF_ANGULAR_MOMENTUM_TRANSFER_MODEL_BSE = 0’, ‘RLOF_ANGULAR_MOMENTUM_TRANSFER_MODEL_CONSERVATIVE = 1’, ‘Random variation : integer between 0 and 1 ‘]

post_SN_orbit_method

Description:

Method by which the post-SN orbit is calculated. 0=BSE, 1=Tauris&Taken 1998.

Parameter input type:

Integer

Default value:

1

Macros:

[‘POST_SN_ORBIT_BSE = 0’, ‘POST_SN_ORBIT_TT98 = 1’]

multiplicity

Description:

Multiplicity: 1=single star, 2=binary, 3=triple, 4=quadruple.

Parameter input type:

Integer

Default value:

0

accretion_limit_eddington_steady_multiplier

Description:

Steady accretion is limited by the Eddington instability, with limiting rate given by the accretion_limit_eddington_steady_multiplier * the normal (spherically symmetric) Eddington rate. This is known in the trade as the Eddington factor, and anything greater than 1.0 potentially gives you super-Eddington accretion.

Parameter input type:

Float

Default value:

1

accretion_limit_eddington_LMMS_multiplier

Description:

Accretion from a low-mass, convective, main_sequence star is limited by the Eddington instability, with limiting rate given by the accretion_limit_eddington_LMMS_multiplier * the normal (spherically symmetric) Eddington rate. This is known in the trade as the Eddington factor, and anything greater than 1.0 potentially gives you super-Eddington accretion.

Parameter input type:

Float

Default value:

1

accretion_limit_eddington_WD_to_remnant_multiplier

Description:

Accretion from a WD onto a remnant star (e.g. another white dwarf, neutron star or black hole) is limited by the Eddington instability, with limiting rate given by the accretion_limit_eddington_WD_to_remnant_multiplier * the normal (spherically symmetric) Eddington rate. This is known in the trade as the Eddington factor, and anything greater than 1.0 potentially gives you super-Eddington accretion.

Parameter input type:

Float

Default value:

-1

accretion_limit_thermal_multiplier

Description:

Mass transfer onto a MS, HG or CHeB star is limited by the accretor’s thermal rate times this multiplier.

Parameter input type:

Float

Default value:

1

accretion_limit_dynamical_multiplier

Description:

Mass transfer is limited by the accretor’s dynamical rate times this multiplier.

Parameter input type:

Float

Default value:

1

donor_limit_envelope_multiplier

Description:

Mass transfer by RLOF is limited by this fraction of the donor’s envelope mass per timestep.

Parameter input type:

Float

Default value:

0

donor_limit_thermal_multiplier

Description:

Mass transfer by RLOF is limited by the accretor’s thermal rate times this multiplier.

Parameter input type:

Float

Default value:

1

donor_limit_dynamical_multiplier

Description:

Mass transfer by RLOF is limited by the donor’s dynamical rate times this multiplier.

Parameter input type:

Float

Default value:

1

Bondi_Hoyle_accretion_factor

Description:

Wind accretion rate, as calculated by the Bondi-Hoyle-Littleton formula, multiplcation factor. Hurley et al 2002 use 1.5, which is the default.

Parameter input type:

Float

Default value:

1.5

tidal_strength_factor

Description:

A modulator for the tidal strength. If this factor > 1 then tides are stronger, i.e. tidal timescales are reduced.

Parameter input type:

Float

Default value:

1

Macros:

[‘TIDAL_STRENGTH_VARIABLE = -1’, ‘Random variation : log-spaced double between 0.0001 and 100 ‘]

tides_variable_filename

Description:

Filename pointing to a table of data to be interpolated to compute the tidal strength multiplier.

Parameter input type:

String

Default value:

hachisu_qcrit

Description:

Critical q=Maccretor/Mdonor above which Hachisu’s disk wind turns on.

Parameter input type:

Float

Default value:

1.15

Macros:

[‘HACHISU_IGNORE_QCRIT = -1’, ‘Random variation : double array or double scalar from [-1] or in range 0.5 to 10’]

hachisu_disk_wind

Description:

Set to True to turn on Hachisu’s disk wind when material accretes too fast onto a white dwarf. This helps to make more SNeIa. See also hachisu_qcrit.

Parameter input type:

True|False

Default value:

False

mass_accretion_for_eld

Description:

The mass that must be accreted onto a COWD for it to ignite as an edge-lit detonation SNIa.

Parameter input type:

Float

Default value:

0.15

WDWD_merger_algorithm

Description:

Algorithm to be used when merging two white dwarfs. 0 = Hurley et al. (2002), 1 = Perets+ (2019), 2 = Chen+ (2016, todo).

Parameter input type:

Integer

Default value:

0

Macros:

[‘WDWD_MERGER_ALGORITHM_BSE = 0’, ‘WDWD_MERGER_ALGORITHM_PERETS2019 = 1’, ‘WDWD_MERGER_ALGORITHM_CHEN2016 = 2’, ‘WDWD_MERGER_ALGORITHM_SATO2016 = 3’, ‘WDWD_MERGER_ALGORITHM_HYBRID_PERETS2019_SATO2016 = 4’, ‘WDWD_MERGER_ALGORITHM_RUITER2013 = 5’]

eta_violent_WDWD_merger

Description:

Energy factor used to calculate q_crit for WDWD mergers that are “violent”. Default = 0.75 (Ruiter 2013).

Parameter input type:

Float

Default value:

0.75

COWD_COWD_explode_above_mass

Description:

COWD-COWD mergers above this mass trigger a SNIa.

Parameter input type:

Float

Default value:

0

HeWD_COWD_explode_above_mass

Description:

HeWD-COWD mergers above this mass trigger a SNIa.

Parameter input type:

Float

Default value:

0

COWD_to_ONeWD_accretion_rate

Description:

Accretion rate above which COWDs are converted to ONeWDs, in Msun/yr. A typical value is 2.05e-6 (Wang et al. 2017). Ignored if 0 (default 0).

Parameter input type:

Float

Default value:

2.05e-06

type_Ia_MCh_supernova_algorithm

Description:

Algorithm to be used when calculating type Ia yields from Chandrasekhar-mass exploders. 0 = DD7 (Iwamoto 1999), 1 = Seitenzahl 2013 3D hydro yields (you must also set Seitenzahl2013_model).

Parameter input type:

Integer

Default value:

NULL

core_collapse_supernova_algorithm

Description:

Algorithm to be used to determine core-collapse supernova yields.

Parameter input type:

Integer

Default value:

NULL

electron_capture_supernova_algorithm

Description:

Algorithm to be used to determine electron-capture supernova yields.

Parameter input type:

Integer

Default value:

NULL

core_collapse_rprocess_algorithm

Description:

Algorithm to be used to determine core-collapse r-process yields. (Default 0 = NUCSYN_CCSNE_RPROCESS_SIMMERER2004.

Parameter input type:

Integer

Default value:

NULL

core_collapse_rprocess_mass

Description:

Mass of r-process material to be associated with a core-collapse supernova.

Parameter input type:

Float

Default value:

NULL

Seitenzahl2013_model

Description:

Which of Seitenzahl et al. 2013’s models to use? One of N1,N3,N5,N10,N20,N40,N100L,N100,N100H,N150,N200,N300C,N1600,N1600C,N100_Z0.5,N100_Z0.1,N100_Z0.01 (defaults to N100).

Parameter input type:

String

Default value:

NULL

Extra:

N1

type_Ia_sub_MCh_supernova_algorithm

Description:

Algorithm to be used when calculating type Ia yields from sub-Chandrasekhar-mass exploders. (Currently unused.).

Parameter input type:

Integer

Default value:

NULL

max_HeWD_mass

Description:

The maximum mass a HeWD can have before it ignites helium (0.7).

Parameter input type:

Float

Default value:

0.7

merger_mass_loss_fraction

Description:

Fraction of the total mass which is lost when stars merge. Note: if merger_mass_loss_fraction_by_stellar_type_<n> is set, where <n> is the stellar type, it is used in place of the value of merger_mass_loss_fraction. Can be a positive number, zero, or an algorithm from: MERGER_MASS_LOSS_FRACTION_GLEBBEEK_2013.

Parameter input type:

Float

Default value:

0.1

Macros:

[‘MERGER_MASS_LOSS_FRACTION_NONE = 0’, ‘MERGER_MASS_LOSS_FRACTION_GLEBBEEK2013 = -1’, ‘MERGER_MASS_LOSS_FRACTION_UNUSED = 666666’]

merger_mass_loss_fraction_nondegenerate

Description:

Equivalent to setting the merger_mass_loss_fraction_by_stellar_type for all nuclear-burning, non-degenerate stellar types (including low-mass main-sequence stars). Note: you should not combine this with attempts to set the mass loss fractions individually.

Parameter input type:
Default value:

NULL

merger_mass_loss_fraction_degenerate

Description:

Equivalent to setting the merger_mass_loss_fraction_by_stellar_type for all (compact) degenerate (i.e. WD, NS) stellar types. Note: you should not combine this with attempts to set the mass loss fractions individually.

Parameter input type:
Default value:

NULL

merger_mass_loss_fraction_by_stellar_type_%d

Description:

Fraction of the total mass which is lost when stars merge to stellar type <n>.

Parameter input type:

Float(scanf)

Default value:

NULL

merger_mass_loss_fraction_by_stellar_type_LOW_MASS_MS

Description:

Fraction of the total mass which is lost when stars merge to stellar type LOW_MASS_MS.

merger_mass_loss_fraction_by_stellar_type_MS

Description:

Fraction of the total mass which is lost when stars merge to stellar type MS.

merger_mass_loss_fraction_by_stellar_type_HG

Description:

Fraction of the total mass which is lost when stars merge to stellar type HG.

merger_mass_loss_fraction_by_stellar_type_GIANT_BRANCH

Description:

Fraction of the total mass which is lost when stars merge to stellar type GIANT_BRANCH.

merger_mass_loss_fraction_by_stellar_type_CHeB

Description:

Fraction of the total mass which is lost when stars merge to stellar type CHeB.

merger_mass_loss_fraction_by_stellar_type_EAGB

Description:

Fraction of the total mass which is lost when stars merge to stellar type EAGB.

merger_mass_loss_fraction_by_stellar_type_TPAGB

Description:

Fraction of the total mass which is lost when stars merge to stellar type TPAGB.

merger_mass_loss_fraction_by_stellar_type_HeMS

Description:

Fraction of the total mass which is lost when stars merge to stellar type HeMS.

merger_mass_loss_fraction_by_stellar_type_HeHG

Description:

Fraction of the total mass which is lost when stars merge to stellar type HeHG.

merger_mass_loss_fraction_by_stellar_type_HeGB

Description:

Fraction of the total mass which is lost when stars merge to stellar type HeGB.

merger_mass_loss_fraction_by_stellar_type_HeWD

Description:

Fraction of the total mass which is lost when stars merge to stellar type HeWD.

merger_mass_loss_fraction_by_stellar_type_COWD

Description:

Fraction of the total mass which is lost when stars merge to stellar type COWD.

merger_mass_loss_fraction_by_stellar_type_ONeWD

Description:

Fraction of the total mass which is lost when stars merge to stellar type ONeWD.

merger_mass_loss_fraction_by_stellar_type_NS

Description:

Fraction of the total mass which is lost when stars merge to stellar type NS.

merger_mass_loss_fraction_by_stellar_type_BH

Description:

Fraction of the total mass which is lost when stars merge to stellar type BH.

merger_mass_loss_fraction_by_stellar_type_MASSLESS_REMNANT

Description:

Fraction of the total mass which is lost when stars merge to stellar type MASSLESS_REMNANT.

transient_method

Description:

Choose the method used to model transients. Options: TRANSIENT_METHOD_NONE, TRANSIENT_METHOD_IVANOVA2013.

Parameter input type:

Integer

Default value:

0

Macros:

[‘TRANSIENT_METHOD_NONE = 0’, ‘TRANSIENT_METHOD_IVANOVA2013 = 1’, ‘TRANSIENT_METHOD_VILLAR2017 = 2’]

merger_angular_momentum_factor

Description:

When two stars merge the resulting single star retains a fraction of the total system angular momentum (or the critical spin angular momentum, if it is smaller) multiplied by this factor.

Parameter input type:

Float

Default value:

1

wind_angular_momentum_loss

Description:

Prescription for losing angular momentum in a stellar wind. 0=Hurley et al 2002 (‘Tout’) prescription, 1=lw i.e. a factor multiplying the specific orbital angular momentum, 2=lw hybrid for fast winds. Set wind_djorb_fac to the desired factor..

Parameter input type:

Integer

Default value:

0

Macros:

[‘WIND_ANGMOM_LOSS_BSE = 0’, ‘WIND_ANGMOM_LOSS_LW = 1’, ‘WIND_ANGMOM_LOSS_LW_HYBRID = 2’, ‘WIND_ANGMOM_LOSS_SPHERICALLY_SYMMETRIC = 3’, ‘Random variation : integer between 0 and 3 ‘]

wind_djorb_fac

Description:

Factor multiplying angular momentum loss in a stellar wind when wind_angular_momentum_loss=0 (the Tout/Hurley et al 2002 prescription). See wind_angular_momentum_loss.

Parameter input type:

Float

Default value:

1

lw

Description:

Factor multiplying angular momentum loss in a stellar wind when wind_angular_momentum_loss=1,2 (the ‘lw’ prescription). See wind_angular_momentum_loss.

Parameter input type:

Float

Default value:

1

VW93_EAGB_wind_speed

Description:

Activate this to use Vassiliadis and Wood (1993) wind speed during the EAGB.

Parameter input type:

True|False

Default value:

False

VW93_TPAGB_wind_speed

Description:

Activate this to use Vassiliadis and Wood (1993) wind speed during the EAGB.

Parameter input type:

True|False

Default value:

False

use_periastron_Roche_radius

Description:

Set this to True to use the Roche lobe radius at periastron, rather than (the default to) assume a circular orbit. This will be useful one day when we treat RLOF in eccentric orbits properly, hopefully.

Parameter input type:

True|False

Default value:

False

qcrit_nuclear_burning

Description:

Set qcrit, the critical mass ratio for stable mass transfer, for all nuclear-burning donors irrespective of the accretor stellar type.

Parameter input type:
Default value:

NULL

Macros:

[‘QCRIT_BSE = -1’, ‘QCRIT_HJELLMING_WEBBINK = -2’, ‘QCRIT_Q_NO_COMENV = -3’, ‘QCRIT_CHEN_HAN_TABLE = -4’, ‘QCRIT_CHEN_HAN_FORMULA = -5’, ‘QCRIT_GE2015 = -6’, ‘QCRIT_VOS2018 = -7’, ‘QCRIT_TEMMINK2022 = -8’, ‘QCRIT_DEFAULT = -9’, ‘QCRIT_GB_BSE = -1’, ‘QCRIT_GB_HJELLMING_WEBBINK = -2’, ‘QCRIT_GB_Q_NO_COMENV = -3’, ‘QCRIT_GB_CHEN_HAN_TABLE = -4’, ‘QCRIT_GB_CHEN_HAN_FORMULA = -5’, ‘QCRIT_GB_GE2015 = -6’, ‘QCRIT_GB_VOS2018 = -7’, ‘QCRIT_GB_TEMMINK2022 = -8’, ‘QCRIT_GB_DEFAULT = -1’]

qcrit_WD

Description:

Set qcrit, the critical mass ratio for stable mass transfer, for all WD donors irrespective of the accretor stellar type.

Parameter input type:
Default value:

NULL

qcrit_all

Description:

Set qcrit, the critical mass ratio for stable mass transfer, for all donors irrespective of donor or accretor stellar type.

Parameter input type:
Default value:

NULL

qcrit_LMMS

Description:

Apply critical q=Mdonor/Maccretor value for low-mass main sequence stars to determine the stability of Roche-lobe overflow for non-degenerate accretors. See also qcrits_*, qcrits_degenerate_*.

Parameter input type:

Float

Default value:

0.6944

Macros:

[‘QCRIT_BSE = -1’, ‘QCRIT_HJELLMING_WEBBINK = -2’, ‘QCRIT_Q_NO_COMENV = -3’, ‘QCRIT_CHEN_HAN_TABLE = -4’, ‘QCRIT_CHEN_HAN_FORMULA = -5’, ‘QCRIT_GE2015 = -6’, ‘QCRIT_VOS2018 = -7’, ‘QCRIT_TEMMINK2022 = -8’, ‘QCRIT_DEFAULT = -9’, ‘QCRIT_GB_BSE = -1’, ‘QCRIT_GB_HJELLMING_WEBBINK = -2’, ‘QCRIT_GB_Q_NO_COMENV = -3’, ‘QCRIT_GB_CHEN_HAN_TABLE = -4’, ‘QCRIT_GB_CHEN_HAN_FORMULA = -5’, ‘QCRIT_GB_GE2015 = -6’, ‘QCRIT_GB_VOS2018 = -7’, ‘QCRIT_GB_TEMMINK2022 = -8’, ‘QCRIT_GB_DEFAULT = -1’]

qcrit_MS

Description:

Apply critical q=Mdonor/Maccretor value for (non-low mass) main sequence stars to determine the stability of Roche-lobe overflow for non-degenerate accretors. See also qcrits_*, qcrits_degenerate_*.

Parameter input type:

Float

Default value:

1.6

Macros:

[‘QCRIT_BSE = -1’, ‘QCRIT_HJELLMING_WEBBINK = -2’, ‘QCRIT_Q_NO_COMENV = -3’, ‘QCRIT_CHEN_HAN_TABLE = -4’, ‘QCRIT_CHEN_HAN_FORMULA = -5’, ‘QCRIT_GE2015 = -6’, ‘QCRIT_VOS2018 = -7’, ‘QCRIT_TEMMINK2022 = -8’, ‘QCRIT_DEFAULT = -9’, ‘QCRIT_GB_BSE = -1’, ‘QCRIT_GB_HJELLMING_WEBBINK = -2’, ‘QCRIT_GB_Q_NO_COMENV = -3’, ‘QCRIT_GB_CHEN_HAN_TABLE = -4’, ‘QCRIT_GB_CHEN_HAN_FORMULA = -5’, ‘QCRIT_GB_GE2015 = -6’, ‘QCRIT_GB_VOS2018 = -7’, ‘QCRIT_GB_TEMMINK2022 = -8’, ‘QCRIT_GB_DEFAULT = -1’]

qcrit_HG

Description:

Apply critical q=Mdonor/Maccretor value for Hertzsprung gap stars to determine the stability of Roche-lobe overflow for non-degenerate accretors. See also qcrits_*, qcrits_degenerate_*.

Parameter input type:

Float

Default value:

4

Macros:

[‘QCRIT_BSE = -1’, ‘QCRIT_HJELLMING_WEBBINK = -2’, ‘QCRIT_Q_NO_COMENV = -3’, ‘QCRIT_CHEN_HAN_TABLE = -4’, ‘QCRIT_CHEN_HAN_FORMULA = -5’, ‘QCRIT_GE2015 = -6’, ‘QCRIT_VOS2018 = -7’, ‘QCRIT_TEMMINK2022 = -8’, ‘QCRIT_DEFAULT = -9’, ‘QCRIT_GB_BSE = -1’, ‘QCRIT_GB_HJELLMING_WEBBINK = -2’, ‘QCRIT_GB_Q_NO_COMENV = -3’, ‘QCRIT_GB_CHEN_HAN_TABLE = -4’, ‘QCRIT_GB_CHEN_HAN_FORMULA = -5’, ‘QCRIT_GB_GE2015 = -6’, ‘QCRIT_GB_VOS2018 = -7’, ‘QCRIT_GB_TEMMINK2022 = -8’, ‘QCRIT_GB_DEFAULT = -1’]

qcrit_GB

Description:

Apply critical q=Mdonor/Maccretor value for first red giant branch stars to determine the stability of Roche-lobe overflow for non-degenerate accretors. See also qcrits_*, qcrits_degenerate_*.

Parameter input type:

Float

Default value:

-1

Macros:

[‘QCRIT_BSE = -1’, ‘QCRIT_HJELLMING_WEBBINK = -2’, ‘QCRIT_Q_NO_COMENV = -3’, ‘QCRIT_CHEN_HAN_TABLE = -4’, ‘QCRIT_CHEN_HAN_FORMULA = -5’, ‘QCRIT_GE2015 = -6’, ‘QCRIT_VOS2018 = -7’, ‘QCRIT_TEMMINK2022 = -8’, ‘QCRIT_DEFAULT = -9’, ‘QCRIT_GB_BSE = -1’, ‘QCRIT_GB_HJELLMING_WEBBINK = -2’, ‘QCRIT_GB_Q_NO_COMENV = -3’, ‘QCRIT_GB_CHEN_HAN_TABLE = -4’, ‘QCRIT_GB_CHEN_HAN_FORMULA = -5’, ‘QCRIT_GB_GE2015 = -6’, ‘QCRIT_GB_VOS2018 = -7’, ‘QCRIT_GB_TEMMINK2022 = -8’, ‘QCRIT_GB_DEFAULT = -1’]

qcrit_CHeB

Description:

Apply critical q=Mdonor/Maccretor value for core helium burning stars to determine the stability of Roche-lobe overflow for non-degenerate accretors. See also qcrits_*, qcrits_degenerate_*.

Parameter input type:

Float

Default value:

3

Macros:

[‘QCRIT_BSE = -1’, ‘QCRIT_HJELLMING_WEBBINK = -2’, ‘QCRIT_Q_NO_COMENV = -3’, ‘QCRIT_CHEN_HAN_TABLE = -4’, ‘QCRIT_CHEN_HAN_FORMULA = -5’, ‘QCRIT_GE2015 = -6’, ‘QCRIT_VOS2018 = -7’, ‘QCRIT_TEMMINK2022 = -8’, ‘QCRIT_DEFAULT = -9’, ‘QCRIT_GB_BSE = -1’, ‘QCRIT_GB_HJELLMING_WEBBINK = -2’, ‘QCRIT_GB_Q_NO_COMENV = -3’, ‘QCRIT_GB_CHEN_HAN_TABLE = -4’, ‘QCRIT_GB_CHEN_HAN_FORMULA = -5’, ‘QCRIT_GB_GE2015 = -6’, ‘QCRIT_GB_VOS2018 = -7’, ‘QCRIT_GB_TEMMINK2022 = -8’, ‘QCRIT_GB_DEFAULT = -1’]

qcrit_EAGB

Description:

Apply critical q=Mdonor/Maccretor value for early-AGB stars to determine the stability of Roche-lobe overflow for non-degenerate accretors. See also qcrits_*, qcrits_degenerate_*.

Parameter input type:

Float

Default value:

-1

Macros:

[‘QCRIT_BSE = -1’, ‘QCRIT_HJELLMING_WEBBINK = -2’, ‘QCRIT_Q_NO_COMENV = -3’, ‘QCRIT_CHEN_HAN_TABLE = -4’, ‘QCRIT_CHEN_HAN_FORMULA = -5’, ‘QCRIT_GE2015 = -6’, ‘QCRIT_VOS2018 = -7’, ‘QCRIT_TEMMINK2022 = -8’, ‘QCRIT_DEFAULT = -9’, ‘QCRIT_GB_BSE = -1’, ‘QCRIT_GB_HJELLMING_WEBBINK = -2’, ‘QCRIT_GB_Q_NO_COMENV = -3’, ‘QCRIT_GB_CHEN_HAN_TABLE = -4’, ‘QCRIT_GB_CHEN_HAN_FORMULA = -5’, ‘QCRIT_GB_GE2015 = -6’, ‘QCRIT_GB_VOS2018 = -7’, ‘QCRIT_GB_TEMMINK2022 = -8’, ‘QCRIT_GB_DEFAULT = -1’]

qcrit_TPAGB

Description:

Apply critical q=Mdonor/Maccretor value for TP-AGB stars to determine the stability of Roche-lobe overflow for non-degenerate accretors. See also qcrits_*, qcrits_degenerate_*.

Parameter input type:

Float

Default value:

-1

Macros:

[‘QCRIT_BSE = -1’, ‘QCRIT_HJELLMING_WEBBINK = -2’, ‘QCRIT_Q_NO_COMENV = -3’, ‘QCRIT_CHEN_HAN_TABLE = -4’, ‘QCRIT_CHEN_HAN_FORMULA = -5’, ‘QCRIT_GE2015 = -6’, ‘QCRIT_VOS2018 = -7’, ‘QCRIT_TEMMINK2022 = -8’, ‘QCRIT_DEFAULT = -9’, ‘QCRIT_GB_BSE = -1’, ‘QCRIT_GB_HJELLMING_WEBBINK = -2’, ‘QCRIT_GB_Q_NO_COMENV = -3’, ‘QCRIT_GB_CHEN_HAN_TABLE = -4’, ‘QCRIT_GB_CHEN_HAN_FORMULA = -5’, ‘QCRIT_GB_GE2015 = -6’, ‘QCRIT_GB_VOS2018 = -7’, ‘QCRIT_GB_TEMMINK2022 = -8’, ‘QCRIT_GB_DEFAULT = -1’]

qcrit_HeMS

Description:

Apply critical q=Mdonor/Maccretor value for helium main sequence stars to determine the stability of Roche-lobe overflow for non-degenerate accretors. See also qcrits_*, qcrits_degenerate_*.

Parameter input type:

Float

Default value:

3

Macros:

[‘QCRIT_BSE = -1’, ‘QCRIT_HJELLMING_WEBBINK = -2’, ‘QCRIT_Q_NO_COMENV = -3’, ‘QCRIT_CHEN_HAN_TABLE = -4’, ‘QCRIT_CHEN_HAN_FORMULA = -5’, ‘QCRIT_GE2015 = -6’, ‘QCRIT_VOS2018 = -7’, ‘QCRIT_TEMMINK2022 = -8’, ‘QCRIT_DEFAULT = -9’, ‘QCRIT_GB_BSE = -1’, ‘QCRIT_GB_HJELLMING_WEBBINK = -2’, ‘QCRIT_GB_Q_NO_COMENV = -3’, ‘QCRIT_GB_CHEN_HAN_TABLE = -4’, ‘QCRIT_GB_CHEN_HAN_FORMULA = -5’, ‘QCRIT_GB_GE2015 = -6’, ‘QCRIT_GB_VOS2018 = -7’, ‘QCRIT_GB_TEMMINK2022 = -8’, ‘QCRIT_GB_DEFAULT = -1’]

qcrit_HeHG

Description:

Apply critical q=Mdonor/Maccretor value for helium Hertzsprung gap stars to determine the stability of Roche-lobe overflow for non-degenerate accretors. See also qcrits_*, qcrits_degenerate_*.

Parameter input type:

Float

Default value:

4

Macros:

[‘QCRIT_BSE = -1’, ‘QCRIT_HJELLMING_WEBBINK = -2’, ‘QCRIT_Q_NO_COMENV = -3’, ‘QCRIT_CHEN_HAN_TABLE = -4’, ‘QCRIT_CHEN_HAN_FORMULA = -5’, ‘QCRIT_GE2015 = -6’, ‘QCRIT_VOS2018 = -7’, ‘QCRIT_TEMMINK2022 = -8’, ‘QCRIT_DEFAULT = -9’, ‘QCRIT_GB_BSE = -1’, ‘QCRIT_GB_HJELLMING_WEBBINK = -2’, ‘QCRIT_GB_Q_NO_COMENV = -3’, ‘QCRIT_GB_CHEN_HAN_TABLE = -4’, ‘QCRIT_GB_CHEN_HAN_FORMULA = -5’, ‘QCRIT_GB_GE2015 = -6’, ‘QCRIT_GB_VOS2018 = -7’, ‘QCRIT_GB_TEMMINK2022 = -8’, ‘QCRIT_GB_DEFAULT = -1’]

qcrit_HeGB

Description:

Apply critical q=Mdonor/Maccretor value for helium red giant stars to determine the stability of Roche-lobe overflow for non-degenerate accretors. See also qcrits_*, qcrits_degenerate_*.

Parameter input type:

Float

Default value:

0.78125

Macros:

[‘QCRIT_BSE = -1’, ‘QCRIT_HJELLMING_WEBBINK = -2’, ‘QCRIT_Q_NO_COMENV = -3’, ‘QCRIT_CHEN_HAN_TABLE = -4’, ‘QCRIT_CHEN_HAN_FORMULA = -5’, ‘QCRIT_GE2015 = -6’, ‘QCRIT_VOS2018 = -7’, ‘QCRIT_TEMMINK2022 = -8’, ‘QCRIT_DEFAULT = -9’, ‘QCRIT_GB_BSE = -1’, ‘QCRIT_GB_HJELLMING_WEBBINK = -2’, ‘QCRIT_GB_Q_NO_COMENV = -3’, ‘QCRIT_GB_CHEN_HAN_TABLE = -4’, ‘QCRIT_GB_CHEN_HAN_FORMULA = -5’, ‘QCRIT_GB_GE2015 = -6’, ‘QCRIT_GB_VOS2018 = -7’, ‘QCRIT_GB_TEMMINK2022 = -8’, ‘QCRIT_GB_DEFAULT = -1’]

qcrit_HeWD

Description:

Apply critical q=Mdonor/Maccretor value for helium white dwarf stars to determine the stability of Roche-lobe overflow for non-degenerate accretors. See also qcrits_*, qcrits_degenerate_*.

Parameter input type:

Float

Default value:

3

Macros:

[‘QCRIT_BSE = -1’, ‘QCRIT_HJELLMING_WEBBINK = -2’, ‘QCRIT_Q_NO_COMENV = -3’, ‘QCRIT_CHEN_HAN_TABLE = -4’, ‘QCRIT_CHEN_HAN_FORMULA = -5’, ‘QCRIT_GE2015 = -6’, ‘QCRIT_VOS2018 = -7’, ‘QCRIT_TEMMINK2022 = -8’, ‘QCRIT_DEFAULT = -9’, ‘QCRIT_GB_BSE = -1’, ‘QCRIT_GB_HJELLMING_WEBBINK = -2’, ‘QCRIT_GB_Q_NO_COMENV = -3’, ‘QCRIT_GB_CHEN_HAN_TABLE = -4’, ‘QCRIT_GB_CHEN_HAN_FORMULA = -5’, ‘QCRIT_GB_GE2015 = -6’, ‘QCRIT_GB_VOS2018 = -7’, ‘QCRIT_GB_TEMMINK2022 = -8’, ‘QCRIT_GB_DEFAULT = -1’]

qcrit_COWD

Description:

Apply critical q=Mdonor/Maccretor value for carbon-oxygen white dwarf stars to determine the stability of Roche-lobe overflow for non-degenerate accretors. See also qcrits_*, qcrits_degenerate_*.

Parameter input type:

Float

Default value:

3

Macros:

[‘QCRIT_BSE = -1’, ‘QCRIT_HJELLMING_WEBBINK = -2’, ‘QCRIT_Q_NO_COMENV = -3’, ‘QCRIT_CHEN_HAN_TABLE = -4’, ‘QCRIT_CHEN_HAN_FORMULA = -5’, ‘QCRIT_GE2015 = -6’, ‘QCRIT_VOS2018 = -7’, ‘QCRIT_TEMMINK2022 = -8’, ‘QCRIT_DEFAULT = -9’, ‘QCRIT_GB_BSE = -1’, ‘QCRIT_GB_HJELLMING_WEBBINK = -2’, ‘QCRIT_GB_Q_NO_COMENV = -3’, ‘QCRIT_GB_CHEN_HAN_TABLE = -4’, ‘QCRIT_GB_CHEN_HAN_FORMULA = -5’, ‘QCRIT_GB_GE2015 = -6’, ‘QCRIT_GB_VOS2018 = -7’, ‘QCRIT_GB_TEMMINK2022 = -8’, ‘QCRIT_GB_DEFAULT = -1’]

qcrit_ONeWD

Description:

Apply critical q=Mdonor/Maccretor value for oxygen-neon white dwarf stars to determine the stability of Roche-lobe overflow for non-degenerate accretors. See also qcrits_*, qcrits_degenerate_*.

Parameter input type:

Float

Default value:

3

Macros:

[‘QCRIT_BSE = -1’, ‘QCRIT_HJELLMING_WEBBINK = -2’, ‘QCRIT_Q_NO_COMENV = -3’, ‘QCRIT_CHEN_HAN_TABLE = -4’, ‘QCRIT_CHEN_HAN_FORMULA = -5’, ‘QCRIT_GE2015 = -6’, ‘QCRIT_VOS2018 = -7’, ‘QCRIT_TEMMINK2022 = -8’, ‘QCRIT_DEFAULT = -9’, ‘QCRIT_GB_BSE = -1’, ‘QCRIT_GB_HJELLMING_WEBBINK = -2’, ‘QCRIT_GB_Q_NO_COMENV = -3’, ‘QCRIT_GB_CHEN_HAN_TABLE = -4’, ‘QCRIT_GB_CHEN_HAN_FORMULA = -5’, ‘QCRIT_GB_GE2015 = -6’, ‘QCRIT_GB_VOS2018 = -7’, ‘QCRIT_GB_TEMMINK2022 = -8’, ‘QCRIT_GB_DEFAULT = -1’]

qcrit_NS

Description:

Apply critical q=Mdonor/Maccretor value for neutron stars to determine the stability of Roche-lobe overflow for non-degenerate accretors. See also qcrits_*, qcrits_degenerate_*.

Parameter input type:

Float

Default value:

3

Macros:

[‘QCRIT_BSE = -1’, ‘QCRIT_HJELLMING_WEBBINK = -2’, ‘QCRIT_Q_NO_COMENV = -3’, ‘QCRIT_CHEN_HAN_TABLE = -4’, ‘QCRIT_CHEN_HAN_FORMULA = -5’, ‘QCRIT_GE2015 = -6’, ‘QCRIT_VOS2018 = -7’, ‘QCRIT_TEMMINK2022 = -8’, ‘QCRIT_DEFAULT = -9’, ‘QCRIT_GB_BSE = -1’, ‘QCRIT_GB_HJELLMING_WEBBINK = -2’, ‘QCRIT_GB_Q_NO_COMENV = -3’, ‘QCRIT_GB_CHEN_HAN_TABLE = -4’, ‘QCRIT_GB_CHEN_HAN_FORMULA = -5’, ‘QCRIT_GB_GE2015 = -6’, ‘QCRIT_GB_VOS2018 = -7’, ‘QCRIT_GB_TEMMINK2022 = -8’, ‘QCRIT_GB_DEFAULT = -1’]

qcrit_BH

Description:

Apply critical q=Mdonor/Maccretor value for black holes to determine the stability of Roche-lobe overflow for non-degenerate accretors. See also qcrits_*, qcrits_degenerate_*.

Parameter input type:

Float

Default value:

3

Macros:

[‘QCRIT_BSE = -1’, ‘QCRIT_HJELLMING_WEBBINK = -2’, ‘QCRIT_Q_NO_COMENV = -3’, ‘QCRIT_CHEN_HAN_TABLE = -4’, ‘QCRIT_CHEN_HAN_FORMULA = -5’, ‘QCRIT_GE2015 = -6’, ‘QCRIT_VOS2018 = -7’, ‘QCRIT_TEMMINK2022 = -8’, ‘QCRIT_DEFAULT = -9’, ‘QCRIT_GB_BSE = -1’, ‘QCRIT_GB_HJELLMING_WEBBINK = -2’, ‘QCRIT_GB_Q_NO_COMENV = -3’, ‘QCRIT_GB_CHEN_HAN_TABLE = -4’, ‘QCRIT_GB_CHEN_HAN_FORMULA = -5’, ‘QCRIT_GB_GE2015 = -6’, ‘QCRIT_GB_VOS2018 = -7’, ‘QCRIT_GB_TEMMINK2022 = -8’, ‘QCRIT_GB_DEFAULT = -1’]

qcrit_degenerate_LMMS

Description:

Apply critical q=Mdonor/Maccretor value for (low mass) main sequence stars to determine the stability of Roche-lobe overflow for degenerate accretors. See also qcrits_*, qcrits_degenerate_*.

Parameter input type:

Float

Default value:

1

Macros:

[‘QCRIT_BSE = -1’, ‘QCRIT_HJELLMING_WEBBINK = -2’, ‘QCRIT_Q_NO_COMENV = -3’, ‘QCRIT_CHEN_HAN_TABLE = -4’, ‘QCRIT_CHEN_HAN_FORMULA = -5’, ‘QCRIT_GE2015 = -6’, ‘QCRIT_VOS2018 = -7’, ‘QCRIT_TEMMINK2022 = -8’, ‘QCRIT_DEFAULT = -9’, ‘QCRIT_GB_BSE = -1’, ‘QCRIT_GB_HJELLMING_WEBBINK = -2’, ‘QCRIT_GB_Q_NO_COMENV = -3’, ‘QCRIT_GB_CHEN_HAN_TABLE = -4’, ‘QCRIT_GB_CHEN_HAN_FORMULA = -5’, ‘QCRIT_GB_GE2015 = -6’, ‘QCRIT_GB_VOS2018 = -7’, ‘QCRIT_GB_TEMMINK2022 = -8’, ‘QCRIT_GB_DEFAULT = -1’]

qcrit_degenerate_MS

Description:

Apply critical q=Mdonor/Maccretor value for (non-low mass) main sequence stars to determine the stability of Roche-lobe overflow for degenerate accretors. See also qcrits_*, qcrits_degenerate_*.

Parameter input type:

Float

Default value:

1

Macros:

[‘QCRIT_BSE = -1’, ‘QCRIT_HJELLMING_WEBBINK = -2’, ‘QCRIT_Q_NO_COMENV = -3’, ‘QCRIT_CHEN_HAN_TABLE = -4’, ‘QCRIT_CHEN_HAN_FORMULA = -5’, ‘QCRIT_GE2015 = -6’, ‘QCRIT_VOS2018 = -7’, ‘QCRIT_TEMMINK2022 = -8’, ‘QCRIT_DEFAULT = -9’, ‘QCRIT_GB_BSE = -1’, ‘QCRIT_GB_HJELLMING_WEBBINK = -2’, ‘QCRIT_GB_Q_NO_COMENV = -3’, ‘QCRIT_GB_CHEN_HAN_TABLE = -4’, ‘QCRIT_GB_CHEN_HAN_FORMULA = -5’, ‘QCRIT_GB_GE2015 = -6’, ‘QCRIT_GB_VOS2018 = -7’, ‘QCRIT_GB_TEMMINK2022 = -8’, ‘QCRIT_GB_DEFAULT = -1’]

qcrit_degenerate_HG

Description:

Apply critical q=Mdonor/Maccretor value for Hertzsprung gap stars to determine the stability of Roche-lobe overflow for degenerate accretors. See also qcrits_*, qcrits_degenerate_*.

Parameter input type:

Float

Default value:

4.7619

Macros:

[‘QCRIT_BSE = -1’, ‘QCRIT_HJELLMING_WEBBINK = -2’, ‘QCRIT_Q_NO_COMENV = -3’, ‘QCRIT_CHEN_HAN_TABLE = -4’, ‘QCRIT_CHEN_HAN_FORMULA = -5’, ‘QCRIT_GE2015 = -6’, ‘QCRIT_VOS2018 = -7’, ‘QCRIT_TEMMINK2022 = -8’, ‘QCRIT_DEFAULT = -9’, ‘QCRIT_GB_BSE = -1’, ‘QCRIT_GB_HJELLMING_WEBBINK = -2’, ‘QCRIT_GB_Q_NO_COMENV = -3’, ‘QCRIT_GB_CHEN_HAN_TABLE = -4’, ‘QCRIT_GB_CHEN_HAN_FORMULA = -5’, ‘QCRIT_GB_GE2015 = -6’, ‘QCRIT_GB_VOS2018 = -7’, ‘QCRIT_GB_TEMMINK2022 = -8’, ‘QCRIT_GB_DEFAULT = -1’]

qcrit_degenerate_GB

Description:

Apply critical q=Mdonor/Maccretor value for first red giant branch stars to determine the stability of Roche-lobe overflow for degenerate accretors. See also qcrits_*, qcrits_degenerate_*.

Parameter input type:

Float

Default value:

1.15

Macros:

[‘QCRIT_BSE = -1’, ‘QCRIT_HJELLMING_WEBBINK = -2’, ‘QCRIT_Q_NO_COMENV = -3’, ‘QCRIT_CHEN_HAN_TABLE = -4’, ‘QCRIT_CHEN_HAN_FORMULA = -5’, ‘QCRIT_GE2015 = -6’, ‘QCRIT_VOS2018 = -7’, ‘QCRIT_TEMMINK2022 = -8’, ‘QCRIT_DEFAULT = -9’, ‘QCRIT_GB_BSE = -1’, ‘QCRIT_GB_HJELLMING_WEBBINK = -2’, ‘QCRIT_GB_Q_NO_COMENV = -3’, ‘QCRIT_GB_CHEN_HAN_TABLE = -4’, ‘QCRIT_GB_CHEN_HAN_FORMULA = -5’, ‘QCRIT_GB_GE2015 = -6’, ‘QCRIT_GB_VOS2018 = -7’, ‘QCRIT_GB_TEMMINK2022 = -8’, ‘QCRIT_GB_DEFAULT = -1’]

qcrit_degenerate_CHeB

Description:

Apply critical q=Mdonor/Maccretor value for core helium burning stars to determine the stability of Roche-lobe overflow for degenerate accretors. See also qcrits_*, qcrits_degenerate_*.

Parameter input type:

Float

Default value:

3

Macros:

[‘QCRIT_BSE = -1’, ‘QCRIT_HJELLMING_WEBBINK = -2’, ‘QCRIT_Q_NO_COMENV = -3’, ‘QCRIT_CHEN_HAN_TABLE = -4’, ‘QCRIT_CHEN_HAN_FORMULA = -5’, ‘QCRIT_GE2015 = -6’, ‘QCRIT_VOS2018 = -7’, ‘QCRIT_TEMMINK2022 = -8’, ‘QCRIT_DEFAULT = -9’, ‘QCRIT_GB_BSE = -1’, ‘QCRIT_GB_HJELLMING_WEBBINK = -2’, ‘QCRIT_GB_Q_NO_COMENV = -3’, ‘QCRIT_GB_CHEN_HAN_TABLE = -4’, ‘QCRIT_GB_CHEN_HAN_FORMULA = -5’, ‘QCRIT_GB_GE2015 = -6’, ‘QCRIT_GB_VOS2018 = -7’, ‘QCRIT_GB_TEMMINK2022 = -8’, ‘QCRIT_GB_DEFAULT = -1’]

qcrit_degenerate_EAGB

Description:

Apply critical q=Mdonor/Maccretor value for early-AGB stars to determine the stability of Roche-lobe overflow for degenerate accretors. See also qcrits_*, qcrits_degenerate_*.

Parameter input type:

Float

Default value:

1.15

Macros:

[‘QCRIT_BSE = -1’, ‘QCRIT_HJELLMING_WEBBINK = -2’, ‘QCRIT_Q_NO_COMENV = -3’, ‘QCRIT_CHEN_HAN_TABLE = -4’, ‘QCRIT_CHEN_HAN_FORMULA = -5’, ‘QCRIT_GE2015 = -6’, ‘QCRIT_VOS2018 = -7’, ‘QCRIT_TEMMINK2022 = -8’, ‘QCRIT_DEFAULT = -9’, ‘QCRIT_GB_BSE = -1’, ‘QCRIT_GB_HJELLMING_WEBBINK = -2’, ‘QCRIT_GB_Q_NO_COMENV = -3’, ‘QCRIT_GB_CHEN_HAN_TABLE = -4’, ‘QCRIT_GB_CHEN_HAN_FORMULA = -5’, ‘QCRIT_GB_GE2015 = -6’, ‘QCRIT_GB_VOS2018 = -7’, ‘QCRIT_GB_TEMMINK2022 = -8’, ‘QCRIT_GB_DEFAULT = -1’]

qcrit_degenerate_TPAGB

Description:

Apply critical q=Mdonor/Maccretor value for TP-AGB stars to determine the stability of Roche-lobe overflow for degenerate accretors. See also qcrits_*, qcrits_degenerate_*.

Parameter input type:

Float

Default value:

1.15

Macros:

[‘QCRIT_BSE = -1’, ‘QCRIT_HJELLMING_WEBBINK = -2’, ‘QCRIT_Q_NO_COMENV = -3’, ‘QCRIT_CHEN_HAN_TABLE = -4’, ‘QCRIT_CHEN_HAN_FORMULA = -5’, ‘QCRIT_GE2015 = -6’, ‘QCRIT_VOS2018 = -7’, ‘QCRIT_TEMMINK2022 = -8’, ‘QCRIT_DEFAULT = -9’, ‘QCRIT_GB_BSE = -1’, ‘QCRIT_GB_HJELLMING_WEBBINK = -2’, ‘QCRIT_GB_Q_NO_COMENV = -3’, ‘QCRIT_GB_CHEN_HAN_TABLE = -4’, ‘QCRIT_GB_CHEN_HAN_FORMULA = -5’, ‘QCRIT_GB_GE2015 = -6’, ‘QCRIT_GB_VOS2018 = -7’, ‘QCRIT_GB_TEMMINK2022 = -8’, ‘QCRIT_GB_DEFAULT = -1’]

qcrit_degenerate_HeMS

Description:

Apply critical q=Mdonor/Maccretor value for helium main sequence stars to determine the stability of Roche-lobe overflow for degenerate accretors. See also qcrits_*, qcrits_degenerate_*.

Parameter input type:

Float

Default value:

3

Macros:

[‘QCRIT_BSE = -1’, ‘QCRIT_HJELLMING_WEBBINK = -2’, ‘QCRIT_Q_NO_COMENV = -3’, ‘QCRIT_CHEN_HAN_TABLE = -4’, ‘QCRIT_CHEN_HAN_FORMULA = -5’, ‘QCRIT_GE2015 = -6’, ‘QCRIT_VOS2018 = -7’, ‘QCRIT_TEMMINK2022 = -8’, ‘QCRIT_DEFAULT = -9’, ‘QCRIT_GB_BSE = -1’, ‘QCRIT_GB_HJELLMING_WEBBINK = -2’, ‘QCRIT_GB_Q_NO_COMENV = -3’, ‘QCRIT_GB_CHEN_HAN_TABLE = -4’, ‘QCRIT_GB_CHEN_HAN_FORMULA = -5’, ‘QCRIT_GB_GE2015 = -6’, ‘QCRIT_GB_VOS2018 = -7’, ‘QCRIT_GB_TEMMINK2022 = -8’, ‘QCRIT_GB_DEFAULT = -1’]

qcrit_degenerate_HeHG

Description:

Apply critical q=Mdonor/Maccretor value for helium Hertzsprung gap stars to determine the stability of Roche-lobe overflow for degenerate accretors. See also qcrits_*, qcrits_degenerate_*.

Parameter input type:

Float

Default value:

4.7619

Macros:

[‘QCRIT_BSE = -1’, ‘QCRIT_HJELLMING_WEBBINK = -2’, ‘QCRIT_Q_NO_COMENV = -3’, ‘QCRIT_CHEN_HAN_TABLE = -4’, ‘QCRIT_CHEN_HAN_FORMULA = -5’, ‘QCRIT_GE2015 = -6’, ‘QCRIT_VOS2018 = -7’, ‘QCRIT_TEMMINK2022 = -8’, ‘QCRIT_DEFAULT = -9’, ‘QCRIT_GB_BSE = -1’, ‘QCRIT_GB_HJELLMING_WEBBINK = -2’, ‘QCRIT_GB_Q_NO_COMENV = -3’, ‘QCRIT_GB_CHEN_HAN_TABLE = -4’, ‘QCRIT_GB_CHEN_HAN_FORMULA = -5’, ‘QCRIT_GB_GE2015 = -6’, ‘QCRIT_GB_VOS2018 = -7’, ‘QCRIT_GB_TEMMINK2022 = -8’, ‘QCRIT_GB_DEFAULT = -1’]

qcrit_degenerate_HeGB

Description:

Apply critical q=Mdonor/Maccretor value for helium red giant stars to determine the stability of Roche-lobe overflow for degenerate accretors. See also qcrits_*, qcrits_degenerate_*.

Parameter input type:

Float

Default value:

1.15

Macros:

[‘QCRIT_BSE = -1’, ‘QCRIT_HJELLMING_WEBBINK = -2’, ‘QCRIT_Q_NO_COMENV = -3’, ‘QCRIT_CHEN_HAN_TABLE = -4’, ‘QCRIT_CHEN_HAN_FORMULA = -5’, ‘QCRIT_GE2015 = -6’, ‘QCRIT_VOS2018 = -7’, ‘QCRIT_TEMMINK2022 = -8’, ‘QCRIT_DEFAULT = -9’, ‘QCRIT_GB_BSE = -1’, ‘QCRIT_GB_HJELLMING_WEBBINK = -2’, ‘QCRIT_GB_Q_NO_COMENV = -3’, ‘QCRIT_GB_CHEN_HAN_TABLE = -4’, ‘QCRIT_GB_CHEN_HAN_FORMULA = -5’, ‘QCRIT_GB_GE2015 = -6’, ‘QCRIT_GB_VOS2018 = -7’, ‘QCRIT_GB_TEMMINK2022 = -8’, ‘QCRIT_GB_DEFAULT = -1’]

qcrit_degenerate_HeWD

Description:

Apply critical q=Mdonor/Maccretor value for helium white dwarf stars to determine the stability of Roche-lobe overflow for degenerate accretors. See also qcrits_*, qcrits_degenerate_*.

Parameter input type:

Float

Default value:

0.625

Macros:

[‘QCRIT_BSE = -1’, ‘QCRIT_HJELLMING_WEBBINK = -2’, ‘QCRIT_Q_NO_COMENV = -3’, ‘QCRIT_CHEN_HAN_TABLE = -4’, ‘QCRIT_CHEN_HAN_FORMULA = -5’, ‘QCRIT_GE2015 = -6’, ‘QCRIT_VOS2018 = -7’, ‘QCRIT_TEMMINK2022 = -8’, ‘QCRIT_DEFAULT = -9’, ‘QCRIT_GB_BSE = -1’, ‘QCRIT_GB_HJELLMING_WEBBINK = -2’, ‘QCRIT_GB_Q_NO_COMENV = -3’, ‘QCRIT_GB_CHEN_HAN_TABLE = -4’, ‘QCRIT_GB_CHEN_HAN_FORMULA = -5’, ‘QCRIT_GB_GE2015 = -6’, ‘QCRIT_GB_VOS2018 = -7’, ‘QCRIT_GB_TEMMINK2022 = -8’, ‘QCRIT_GB_DEFAULT = -1’]

qcrit_degenerate_COWD

Description:

Apply critical q=Mdonor/Maccretor value for carbon-oxygen white dwarf stars to determine the stability of Roche-lobe overflow for degenerate accretors. See also qcrits_*, qcrits_degenerate_*.

Parameter input type:

Float

Default value:

0.625

Macros:

[‘QCRIT_BSE = -1’, ‘QCRIT_HJELLMING_WEBBINK = -2’, ‘QCRIT_Q_NO_COMENV = -3’, ‘QCRIT_CHEN_HAN_TABLE = -4’, ‘QCRIT_CHEN_HAN_FORMULA = -5’, ‘QCRIT_GE2015 = -6’, ‘QCRIT_VOS2018 = -7’, ‘QCRIT_TEMMINK2022 = -8’, ‘QCRIT_DEFAULT = -9’, ‘QCRIT_GB_BSE = -1’, ‘QCRIT_GB_HJELLMING_WEBBINK = -2’, ‘QCRIT_GB_Q_NO_COMENV = -3’, ‘QCRIT_GB_CHEN_HAN_TABLE = -4’, ‘QCRIT_GB_CHEN_HAN_FORMULA = -5’, ‘QCRIT_GB_GE2015 = -6’, ‘QCRIT_GB_VOS2018 = -7’, ‘QCRIT_GB_TEMMINK2022 = -8’, ‘QCRIT_GB_DEFAULT = -1’]

qcrit_degenerate_ONeWD

Description:

Apply critical q=Mdonor/Maccretor value for oxygen-neon white dwarf stars to determine the stability of Roche-lobe overflow for degenerate accretors. See also qcrits_*, qcrits_degenerate_*.

Parameter input type:

Float

Default value:

0.625

Macros:

[‘QCRIT_BSE = -1’, ‘QCRIT_HJELLMING_WEBBINK = -2’, ‘QCRIT_Q_NO_COMENV = -3’, ‘QCRIT_CHEN_HAN_TABLE = -4’, ‘QCRIT_CHEN_HAN_FORMULA = -5’, ‘QCRIT_GE2015 = -6’, ‘QCRIT_VOS2018 = -7’, ‘QCRIT_TEMMINK2022 = -8’, ‘QCRIT_DEFAULT = -9’, ‘QCRIT_GB_BSE = -1’, ‘QCRIT_GB_HJELLMING_WEBBINK = -2’, ‘QCRIT_GB_Q_NO_COMENV = -3’, ‘QCRIT_GB_CHEN_HAN_TABLE = -4’, ‘QCRIT_GB_CHEN_HAN_FORMULA = -5’, ‘QCRIT_GB_GE2015 = -6’, ‘QCRIT_GB_VOS2018 = -7’, ‘QCRIT_GB_TEMMINK2022 = -8’, ‘QCRIT_GB_DEFAULT = -1’]

qcrit_degenerate_NS

Description:

Apply critical q=Mdonor/Maccretor value for neutron stars to determine the stability of Roche-lobe overflow for degenerate accretors. See also qcrits_*, qcrits_degenerate_*.

Parameter input type:

Float

Default value:

0.625

Macros:

[‘QCRIT_BSE = -1’, ‘QCRIT_HJELLMING_WEBBINK = -2’, ‘QCRIT_Q_NO_COMENV = -3’, ‘QCRIT_CHEN_HAN_TABLE = -4’, ‘QCRIT_CHEN_HAN_FORMULA = -5’, ‘QCRIT_GE2015 = -6’, ‘QCRIT_VOS2018 = -7’, ‘QCRIT_TEMMINK2022 = -8’, ‘QCRIT_DEFAULT = -9’, ‘QCRIT_GB_BSE = -1’, ‘QCRIT_GB_HJELLMING_WEBBINK = -2’, ‘QCRIT_GB_Q_NO_COMENV = -3’, ‘QCRIT_GB_CHEN_HAN_TABLE = -4’, ‘QCRIT_GB_CHEN_HAN_FORMULA = -5’, ‘QCRIT_GB_GE2015 = -6’, ‘QCRIT_GB_VOS2018 = -7’, ‘QCRIT_GB_TEMMINK2022 = -8’, ‘QCRIT_GB_DEFAULT = -1’]

qcrit_degenerate_BH

Description:

Apply critical q=Mdonor/Maccretor value for black holes to determine the stability of Roche-lobe overflow for degenerate accretors. See also qcrits_*, qcrits_degenerate_*.

Parameter input type:

Float

Default value:

0.625

Macros:

[‘QCRIT_BSE = -1’, ‘QCRIT_HJELLMING_WEBBINK = -2’, ‘QCRIT_Q_NO_COMENV = -3’, ‘QCRIT_CHEN_HAN_TABLE = -4’, ‘QCRIT_CHEN_HAN_FORMULA = -5’, ‘QCRIT_GE2015 = -6’, ‘QCRIT_VOS2018 = -7’, ‘QCRIT_TEMMINK2022 = -8’, ‘QCRIT_DEFAULT = -9’, ‘QCRIT_GB_BSE = -1’, ‘QCRIT_GB_HJELLMING_WEBBINK = -2’, ‘QCRIT_GB_Q_NO_COMENV = -3’, ‘QCRIT_GB_CHEN_HAN_TABLE = -4’, ‘QCRIT_GB_CHEN_HAN_FORMULA = -5’, ‘QCRIT_GB_GE2015 = -6’, ‘QCRIT_GB_VOS2018 = -7’, ‘QCRIT_GB_TEMMINK2022 = -8’, ‘QCRIT_GB_DEFAULT = -1’]

mass_for_Hestar_Ia_upper

Description:

Only helium stars below this mass can explode as SNIa. Default is zero, i.e. it never happens. See also mass_for_Hestar_Ia_lower.

Parameter input type:

Float

Default value:

0

mass_for_Hestar_Ia_lower

Description:

Only helium stars above this mass can explode as SNIa. Default is zero, i.e. it never happens. See also mass_for_Hestar_Ia_upper.

Parameter input type:

Float

Default value:

0

alphaCB

Description:

Circumbinary disk viscosity parameter, alpha.

Parameter input type:

Float

Default value:

NULL

minimum_donor_menv_for_comenv

Description:

Minimum donor envelope mass for common envelope evolution to be triggered (Msun). Default 0.

Parameter input type:

Float

Default value:

0

comenv_prescription

Description:

Use this to choose which common envelope prescription you should use. 0=Hurley et al 2002 (based on the Paczyński energy model) or 1=Nelemans and Tout (angular momentum model). See also alpha_ce, comenv_ms_accretion_mass, comenv_ms_accretion_fraction, comenv_ns_accretion_fraction, comenv_ns_accretion_mass, nelemans_gamma, nelemans_minq, nelemans_max_frac_j_change, nelemans_n_comenvs, lambda_ce, lambda_ionisation.

Parameter input type:

Integer

Default value:

0

Macros:

[‘COMENV_UNDEF = -1’, ‘COMENV_BSE = 0’, ‘COMENV_NELEMANS_TOUT = 1’, ‘COMENV_NANDEZ2016 = 2’, ‘COMENV_GE2022 = 3’, ‘Random variation : integer between 0 and 1 ‘]

comenv_prescription%d

Description:

Use this to choose which common envelope prescription you should use. 0=Hurley et al 2002 (based on the Paczyński energy model) or 1=Nelemans and Tout (angular momentum model). See also alpha_ce, comenv_ms_accretion_mass, comenv_ms_accretion_fraction, comenv_ns_accretion_fraction, comenv_ns_accretion_mass, nelemans_gamma, nelemans_minq, nelemans_max_frac_j_change, nelemans_n_comenvs, lambda_ce, lambda_ionisation.

Parameter input type:

Int(scanf)

Default value:

NULL

Macros:

[‘COMENV_UNDEF = -1’, ‘COMENV_BSE = 0’, ‘COMENV_NELEMANS_TOUT = 1’, ‘COMENV_NANDEZ2016 = 2’, ‘COMENV_GE2022 = 3’]

comenv_ejection_spin_method

Description:

When a common envelope is ejected, we need to decide how fast the stars are left spinning. COMENV_EJECTION_SPIN_METHOD_DO_NOTHING (0) is the default, this just leaves the stars/stellar cores spinning with the same spin rate (omega = angular velocity) with which they entered the common envelope phase. COMENV_EJECTION_SPIN_METHOD_SYCHRONIZE instead tidally synchronizes the stars with their new orbital angular velocity.

Parameter input type:

Integer

Default value:

0

Macros:

[‘COMENV_EJECTION_SPIN_METHOD_DO_NOTHING = 0’, ‘COMENV_EJECTION_SPIN_METHOD_SYNCHRONIZE = 1’, ‘Random variation : integer between 0 and 1 ‘]

comenv_merger_spin_method

Description:

When a common envelope binary merges, we need to decide how fast the resulting single star is left spinning. COMENV_MERGER_SPIN_METHOD_SPECIFIC (0) is the default, this preserves angular momentum but limits the specific angular momentum of the merged star to the specific angular momentum of the system at the onset of common envelope evolution. COMENV_MERGER_SPIN_METHOD_CONSERVE_ANGMOM (1) sets the merger’s angular momentum to be that of the system at the onset of common envelope evolution (which means the star may be rotating supercritically). COMENV_MERGER_SPIN_METHOD_CONSERVE_OMEGA (2) sets the spin rate (angular frequency = omega) of the merged star to be that of the orbit just at the onset of common envelope evolution.

Parameter input type:

Integer

Default value:

0

Macros:

[‘COMENV_MERGER_SPIN_METHOD_SPECIFIC = 0’, ‘COMENV_MERGER_SPIN_METHOD_CONSERVE_ANGMOM = 1’, ‘COMENV_MERGER_SPIN_METHOD_CONSERVE_OMEGA = 2’, ‘COMENV_MERGER_SPIN_METHOD_BREAKUP = 3’, ‘Random variation : integer between 0 and 2 ‘]

comenv_ms_accretion_mass

Description:

Experimental. During common envelope evolution, a main sequence star may accrete some of the envelope’s mass. Requires COMENV_MS_ACCRETION. See also comenv_ms_accretion_fraction.

Parameter input type:

Float

Default value:

NULL

comenv_ms_accretion_fraction

Description:

Experimental. During common envelope evolution, a main sequence may accrete a fraction of the envelope’s mass. Requires COMENV_MS_ACCRETION. See also comenv_ms_accretion_mass.

Parameter input type:

Float

Default value:

NULL

comenv_ns_accretion_mass

Description:

Experimental. During common envelope evolution, a neutron star may accrete some of the envelope’s mass. Requires COMENV_NS_ACCRETION. See also comenv_ns_accretion_fraction.

Parameter input type:

Float

Default value:

NULL

comenv_ns_accretion_fraction

Description:

Experimental. During common envelope evolution, a neutron star may accrete a fraction of the envelope’s mass. Requires COMENV_NS_ACCRETION. See also comenv_ns_accretion_mass.

Parameter input type:

Float

Default value:

NULL

alpha_ce

Description:

Common envelope energy formalism parameter. A fraction alpha of the orbital energy is used to eject the envelope. See Hurley et al 2002 for details.

Parameter input type:

Float

Default value:

1

alpha_ce%d

Description:

Common envelope energy formalism parameter. A fraction alpha of the orbital energy is used to eject the envelope. See Hurley et al. 2002 for details.

Parameter input type:

Float(scanf)

Default value:

NULL

lambda_ce

Description:

Common envelope parameter. The binding energy of the common envelope is G*M*Menv/(lambda*R). Typically this is taken to be 0.5, but if set to LAMBDA_CE_DEWI_TAURIS == -1 binary_c uses the Dewi and Tauris fits instead, LAMBDA_CE_WANG_2016 == -2 uses the formalism of Wang, Jia and Li (2016), if LAMBDA_CE_POLYTROPE == -3 then a polytropic formalism is used (see also comenv_splitmass) and if LAMBDA_CE_KLENCKI_2020 == -4 use Klencki et al. (2020).

Parameter input type:

Float

Default value:

0.5

Macros:

[‘LAMBDA_CE_DEWI_TAURIS = -1’, ‘LAMBDA_CE_WANG_2016 = -2’, ‘LAMBDA_CE_POLYTROPE = -3’, ‘LAMBDA_CE_KLENCKI_2020 = -4’]

lambda_ce%d

Description:

Common envelope parameter. The binding energy of the common envelope is G*M*Menv/(lambda*R). Typically this is taken to be 0.5, but if set to -1 binary_c uses the Dewi and Tauris fits instead, -2 uses the formalism of Wang, Jia and Li (2016) and if -3 then a polytropic formalism is used (see also comenv_splitmass).

Parameter input type:

Float(scanf)

Default value:

NULL

comenv_splitmass

Description:

When lambda_ce=-2, the envelope binding energy, lambda, is calculated using a polytropic formalism. The comenv_splitmass defines the point, in the units of the core mass, above which material is ejected.

Parameter input type:

Float

Default value:

NULL

nelemans_recalc_eccentricity

Description:

If True, recalculate the eccentricity after angular momentum is removed.

Parameter input type:

True|False

Default value:

False

comenv_post_eccentricity

Description:

Eccentricity remaining after common envelope ejection.

Parameter input type:

Float

Default value:

1e-05

dedmRLOF

Description:

De/dM during RLOF (defualt 0.0).

Parameter input type:

Float

Default value:

0

nelemans_gamma

Description:

Set the fraction of the orbital specific angular momentum that is used to eject the common envelope according to the Nelemans and Tout prescription. See also nelemans_minq, nelemans_max_frac_j_change, nelemans_n_comenvs.

Parameter input type:

Float

Default value:

1.75

nelemans_minq

Description:

Only activate the Nelemans and Tout common envelope prescription for q>nelemans_minq. See also nelemans_gamma, nelemans_max_frac_j_change, nelemans_n_comenvs.

Parameter input type:

Float

Default value:

0.2

nelemans_max_frac_j_change

Description:

Maximum fractional angular momentum change in the Nelemans and Tout common envelope prescription. See also nelemans_gamma, nelemans_minq, nelemans_n_comenvs.

Parameter input type:

Float

Default value:

1

nelemans_n_comenvs

Description:

Set the maximum number of common envelope ejections allowed to follow the Nelemans and Tout prescription, after which the standard alpha prescription is used.

Parameter input type:

Integer

Default value:

1

lambda_ionisation

Description:

A fraction lambda_ionisation of the recombination energy in the common envelope goes into ejecting the envelope. This is usually 0.0, but a positive value can make a big difference to the outcome of common envelope evolution.

Parameter input type:

Float

Default value:

0.5

lambda_ionisation%d

Description:

A fraction lambda_ionisation of the recombination energy in the common envelope goes into ejecting the envelope. This is usually 0.0, but a positive value can make a big difference to the outcome of common envelope evolution.

Parameter input type:

Float(scanf)

Default value:

NULL

lambda_enthalpy

Description:

A fraction of the enthalpy to be included in the common envelope evolution binding energy. Only used for the Wang 2016 prescription (so far).

Parameter input type:

Float

Default value:

0

lambda_enthalpy%d

Description:

A fraction of the enthalpy to be included in the common envelope evolution binding energy. Only used for the Wang 2016 prescription (so far).

Parameter input type:

Float(scanf)

Default value:

NULL

cbdisc_albedo

Description:

Circumbinary-disc albedo. Requires DISCS.

Parameter input type:

Float

Default value:

0

cbdisc_gamma

Description:

Circumbinary disc gamma (equation of state) parameter. Requires DISCS.

Parameter input type:

Float

Default value:

1.4

cbdisc_alpha

Description:

Circumbinary disc alpha (viscosity) parameter. Requires DISCS.

Parameter input type:

Float

Default value:

0.001

cbdisc_kappa

Description:

Circumbinary disc kappa (opacity) parameter. Requires DISCS.

Parameter input type:

Float

Default value:

0.01

cbdisc_minimum_evaporation_timescale

Description:

Circumbinary disc minimum evaporation timescale (years). If (slow, not edge stripped) mass loss would evaporate the disc on a timescale less than this, simply evaporate the disc immediated. Usually set to 1y, ignore if zero. Requires DISCS.

Parameter input type:

Float

Default value:

1

cbdisc_torquef

Description:

Circumbinary disc binary torque multiplier. Requires DISCS.

Parameter input type:

Float

Default value:

0.001

cbdisc_max_lifetime

Description:

Circumbinary disc maximum lifetime (years, ignored if 0). Requires DISCS.

Parameter input type:

Float

Default value:

1e+06

cbdisc_init_dM

Description:

On cbdisc start, reduce mass by a fraction dM if it won’t converge. Requires DISCS.

Parameter input type:

Float

Default value:

0.1

cbdisc_init_dJdM

Description:

On cbdisc start, reduce angular momentum by a fraction dJ/dM*dM if it won’t converge. Requires DISCS.

Parameter input type:

Float

Default value:

0.5

cbdisc_mass_loss_constant_rate

Description:

Circumbinary disc constant mass loss rate (Msun/year). Requires DISCS.

Parameter input type:

Float

Default value:

0

cbdisc_mass_loss_FUV_multiplier

Description:

Circumbinary disc FUV mass loss rate multiplier (no units). Requires DISCS.

Parameter input type:

Float

Default value:

1

cbdisc_mass_loss_Xray_multiplier

Description:

Circumbinary disc X-ray mass loss rate multiplier (no units). Requires DISCS.

Parameter input type:

Float

Default value:

1

cbdisc_mass_loss_ISM_ram_pressure_multiplier

Description:

Circumbinary disc interstellar medium ram pressure stripping mass loss rate multiplier (no units). Requires DISCS.

Parameter input type:

Float

Default value:

1

cbdisc_mass_loss_ISM_pressure

Description:

Circumbinary disc interstellar medium ram pressure in units of Boltzmann constant per Kelvin (I think…). Requires DISCS. Typically 3000.0. See e.g. http://www.astronomy.ohio-state.edu/~pogge/Ast871/Notes/Intro.pdf page 15 or https://arxiv.org/pdf/0902.0820.pdf Fig. 1 (left panel).

Parameter input type:

Float

Default value:

3000

cbdisc_mass_loss_inner_viscous_multiplier

Description:

Circumbinary disc inner edge viscous mass loss rate multiplier (no units). Requires DISCS.

Parameter input type:

Float

Default value:

1

cbdisc_mass_loss_inner_viscous_angular_momentum_multiplier

Description:

Circumbinary disc inner edge viscous angular momentum multiplier (no units). The inner edge angular momentum Requires DISCS.

Parameter input type:

Float

Default value:

1

cbdisc_resonance_multiplier

Description:

Circumbinary disc resonant interaction multiplier, affects eccentricity pumping and angular momentum rates. Requires DISCS.

Parameter input type:

Float

Default value:

1

cbdisc_resonance_damping

Description:

Circumbinary disc resonant interaction damping: should be on (True) to damp the l=1, m=2 resonance when the disc inner edge lies outside the resonance location. Requires DISCS.

Parameter input type:

True|False

Default value:

True

cbdisc_fail_ring_inside_separation

Description:

If True, while converging on a structure, circumbinary discs with Rring < the binary separation are immediately failed.

Parameter input type:

True|False

Default value:

False

cbdisc_mass_loss_inner_L2_cross_multiplier

Description:

Circumbinary disc inner edge L2-crossing mass loss rate multiplier (no units). Requires DISCS.

Parameter input type:

Float

Default value:

1

cbdisc_minimum_luminosity

Description:

Circumbinary disc minimum luminosity. If the disc becomes dimmer than this, the disc is evaporated instantly. Requires DISCS.

Parameter input type:

Float

Default value:

0

cbdisc_minimum_mass

Description:

Circumbinary disc minimum mass. If the disc becomes less massive than this, the disc is evaporated instantly. Requires DISCS.

Parameter input type:

Float

Default value:

1e-06

cbdisc_minimum_fRing

Description:

Circumbinary disc minimum fRing. If the disc becomes a ring, and fRing = abs(Rout/Rin-1) < this value (and this value is non-zero), the disc is evaporated instantly. Requires DISCS.

Parameter input type:

Float

Default value:

0.2

comenv_disc_angmom_fraction

Description:

If >0 Fraction of the common envelope’s angular momentum that goes into the circumbinary disc. If -1 then uses the moments of inertia to calculate (deprecated), if -2 use the common envelope’s specific angular momentum, if -3 uses the L2 point at the end of the common envelope to set the angular momentum. Requires DISCS and DISCS_CIRCUMBINARY_FROM_COMENV.

Parameter input type:

Float

Default value:

0

comenv_disc_mass_fraction

Description:

Fraction of the common envelope’s mass that goes into the circumbinary disc. Requires DISCS and DISCS_CIRCUMBINARY_FROM_COMENV.

Parameter input type:

Float

Default value:

0

wind_disc_angmom_fraction

Description:

If >0 Fraction of the wind envelope’s angular momentum that goes into the circumbinary disc. If -1 then uses the L2 point’s specific angular momentum. Requires DISCS and DISCS_CIRCUMBINARY_FROM_WIND.

Parameter input type:

Float

Default value:

0

wind_disc_mass_fraction

Description:

Fraction of the stellar wind’s mass that goes into the circumbinary disc. Requires DISCS and DISCS_CIRCUMBINARY_FROM_WIND.

Parameter input type:

Float

Default value:

0

WRLOF_method

Description:

Choose whether and how to apply wind-Roche-lobe-overflow. 0=none, 1=q-dependent, 2=quadratic See Abate et al 2013/14 for details. Requires WRLOF_MASS_TRANSFER.

Parameter input type:

Integer

Default value:

0

Macros:

[‘WRLOF_NONE = 0’, ‘WRLOF_Q_DEPENDENT = 1’, ‘WRLOF_QUADRATIC = 2’]

minimum_timestep

Description:

The minimum timestep (Myr).

Parameter input type:

Float

Default value:

1e-06

timestep_solver_factor

Description:

Factor applied in timestep_limits, e.g. to prevent X changing too fast by comparing to X/dX/dt, which is usually 1 but can be higher to lengthen timesteps when using an alternative solver.

Parameter input type:

Float

Default value:

1

RLOF_mdot_factor

Description:

Multiplier applied to the mass transfer rate during Roche-lobe overflow. Requires RLOF_MDOT_MODULATION.

Parameter input type:

Float

Default value:

1

RLOF_f

Description:

Factor to enlarge a Roche lobe, nominally because of radiation pressure (see Dermine et al paper). Requires RLOF_RADIATION_CORRECTION.

Parameter input type:

Float

Default value:

NULL

unstable_RLOF_can_trigger_SNIa

Description:

If true, allow unstable (dynamical) RLOF to trigger a sub-MCh SNIa. (False).

Parameter input type:

True|False

Default value:

False

triggered_SNIa_algorithm

Description:

Algorithm to use for triggered SNeIa.

Parameter input type:

Integer

Default value:

NULL

minimum_separation_for_instant_RLOF

Description:

If True, instead of evolving the system just report the minimum separation (on the zero-age main sequence) that would lead to instant RLOF. Used by binary_grid. See also minimum_orbital_period_for_instant_RLOF and maximum_mass_ratio_for_instant_RLOF.

Parameter input type:

True|False

Default value:

False

minimum_orbital_period_for_instant_RLOF

Description:

If True, instead of evolving the system just report the minimum orbital period (on the zero-age main sequence) that would lead to instant RLOF. Used by binary_grid. See also minimum_separation_for_instant_RLOF and maximum_mass_ratio_for_instant_RLOF.

Parameter input type:

True|False

Default value:

False

maximum_mass_ratio_for_instant_RLOF

Description:

If True, instead of evolving the system just report the maximum mass ratio (on the zero-age main sequence) that would lead to instant RLOF, given M1 and orbital period. Used by binary_grid. See also minimum_separation_for_instant_RLOF and minimum_orbital_period_for_instant_RLOF.

Parameter input type:

True|False

Default value:

False

RLOF_method

Description:

Use RLOF_method to choose the algorithm you use for Roche-lobe overflow mass loss rate calculations. 0=Hurley et al 2002, 1=Adaptive (for radiative stars) R=RL method, 2=Ritter (probably broken), 3=Claeys etal 2014 variant on Hurley et al 2002.

Parameter input type:

Integer

Default value:

0

Macros:

[‘RLOF_METHOD_BSE = 0’, ‘RLOF_METHOD_ADAPTIVE = 1’, ‘RLOF_METHOD_RITTER = 2’, ‘RLOF_METHOD_CLAEYS = 3’, ‘RLOF_METHOD_ADAPTIVE2 = 4’, ‘Random variation : int array from 0,3’]

RLOF_interpolation_method

Description:

When a star overflows its Roche lobe, it always has R>RL because of the limited time resolution of the simulation. Binary_c then uses an algorithm to get back to when R~RL (within a desired tolerance, set in RLOF_ENTRY_THRESHOLD which is usually 1.02, i.e. overflow of 2%). You can choose algorithm 0, the Hurley et al 2002 method which reverses time (i.e. uses a Newton-like scheme), or 1 to use the binary_c method which rejects a timestep (and hence does no logging on that timestep) and repeats with half the timestep until R~RL. The latter is now the default, because this means there are no negative timesteps which break various other algorithms (e.g. nucleosynthesis).

Parameter input type:

Integer

Default value:

0

Macros:

[‘RLOF_INTERPOLATION_BINARY_C = 0’, ‘RLOF_INTERPOLATION_BSE = 1’]

He_nova_ELDs

Description:

Enable edge-lit detonations from WDs that have He novae. (Default: False, i.e. off.).

Parameter input type:

True|False

Default value:

False

nova_retention_fraction_H

Description:

The mass accreted during a hydrogen nova as fraction of mass transferred.

Parameter input type:

Float

Default value:

0.001

nova_retention_fraction_He

Description:

The mass accreted during a helium nova as fraction of mass transferred.

Parameter input type:

Float

Default value:

0.001

beta_reverse_nova

Description:

The fraction of mass ejected in a nova explosion which is accreted back onto the companion star. Set to -1 to automatically calculate based on a geometric argument, or 0 or positive to set the value.

Parameter input type:

Float

Default value:

0

Macros:

[‘BETA_REVERSE_NOVAE_GEOMETRY = -1’]

WD_accretion_rate_nalgorithms

Description:

Wrapper to set all WD accretion rates in one go (see WD_accretion_rate_*).

Parameter input type:
Default value:

NULL

WD_accretion_rate_novae_upper_limit_hydrogen_donor

Description:

Upper limit of the stable mass transfer rate onto a white dwarf that leads to novae when the donor is hydrogen rich: above this rate the mass transfer leads to stable burning.

Parameter input type:

Float

Default value:

-1

Macros:

[‘DONOR_RATE_ALGORITHM_CLAEYS2014 = -1’, ‘DONOR_RATE_ALGORITHM_BSE = -2’, ‘DONOR_RATE_ALGORITHM_K2014_P2014 = -3’, ‘DONOR_RATE_ALGORITHM_WANGWU = -4’, ‘Random variation : log-spaced double between 1e-10 and 1e-05 ‘]

WD_accretion_rate_novae_upper_limit_helium_donor

Description:

Upper limit of the stable mass transfer rate onto a white dwarf that leads to novae when the donor is helium rich: above this rate the mass transfer leads to stable burning.

Parameter input type:

Float

Default value:

-1

Macros:

[‘DONOR_RATE_ALGORITHM_CLAEYS2014 = -1’, ‘DONOR_RATE_ALGORITHM_BSE = -2’, ‘DONOR_RATE_ALGORITHM_K2014_P2014 = -3’, ‘DONOR_RATE_ALGORITHM_WANGWU = -4’, ‘Random variation : log-spaced double between 1e-10 and 1e-05 ‘]

WD_accretion_rate_novae_upper_limit_other_donor

Description:

Upper limit of the stable mass transfer rate onto a white dwarf that leads to novae when the donor is neither hydrogen nor helium rich: above this rate the mass transfer leads to stable burning.

Parameter input type:

Float

Default value:

-1

Macros:

[‘DONOR_RATE_ALGORITHM_CLAEYS2014 = -1’, ‘DONOR_RATE_ALGORITHM_BSE = -2’, ‘DONOR_RATE_ALGORITHM_K2014_P2014 = -3’, ‘DONOR_RATE_ALGORITHM_WANGWU = -4’, ‘Random variation : log-spaced double between 1e-10 and 1e-05 ‘]

WD_accretion_rate_new_giant_envelope_lower_limit_hydrogen_donor

Description:

Lower limit of the mass transfer rate onto a white dwarf that leads to a the formation of a new giant envelope with a hydrogen-rich donor. Below this mass transfer leads to stable burning.

Parameter input type:

Float

Default value:

-1

Macros:

[‘DONOR_RATE_ALGORITHM_CLAEYS2014 = -1’, ‘DONOR_RATE_ALGORITHM_BSE = -2’, ‘DONOR_RATE_ALGORITHM_K2014_P2014 = -3’, ‘DONOR_RATE_ALGORITHM_WANGWU = -4’]

WD_accretion_rate_new_giant_envelope_lower_limit_helium_donor

Description:

Lower limit of the mass transfer rate onto a white dwarf that leads to a the formation of a new giant envelope with a helium-rich donor. Below this mass transfer leads to stable burning.

Parameter input type:

Float

Default value:

-1

Macros:

[‘DONOR_RATE_ALGORITHM_CLAEYS2014 = -1’, ‘DONOR_RATE_ALGORITHM_BSE = -2’, ‘DONOR_RATE_ALGORITHM_K2014_P2014 = -3’, ‘DONOR_RATE_ALGORITHM_WANGWU = -4’]

WD_accretion_rate_new_giant_envelope_lower_limit_other_donor

Description:

Lower limit of the mass transfer rate onto a white dwarf that leads to a the formation of a new giant envelope when the donor is neither hydrogen nor helium rich. Below this mass transfer leads to stable burning.

Parameter input type:

Float

Default value:

-1

Macros:

[‘DONOR_RATE_ALGORITHM_CLAEYS2014 = -1’, ‘DONOR_RATE_ALGORITHM_BSE = -2’, ‘DONOR_RATE_ALGORITHM_K2014_P2014 = -3’, ‘DONOR_RATE_ALGORITHM_WANGWU = -4’]

CRAP_parameter

Description:

Tidally enhanced mass loss parameter. See Tout and Eggleton’s paper on the subject. (Was the parameter bb).

Parameter input type:

Float

Default value:

0

individual_novae

Description:

If individual_novae is True, novae are resolved such that each explosion is performed separtaely.

Parameter input type:

True|False

Default value:

False

nova_timestep_accelerator_num

Description:

The nova timestep is accelerated if the nova number exceeds nova_timestep_accelerator_num. If zero or negative, acceleration is off. See also nova_timestep_accelerator_index and nova_timestep_accelerator_max. Only used if individual_novae is on.

Parameter input type:

Float

Default value:

100

nova_timestep_accelerator_index

Description:

The index at which the nova timestep is accelerated. A larger value gives longer timesteps. See also nova_timestep_accelerator_num and nova_timestep_accelerator_max. Only used if individual_novae is on.

Parameter input type:

Float

Default value:

0.5

nova_timestep_accelerator_max

Description:

The nova timestep is accelerated by a factor that is capped at nova_timestep_accelerator_max. This parameter is ignored if it is zero or negative. See also nova_timestep_accelerator_num and nova_timestep_accelerator_index. Only used if individual_novae is on.

Parameter input type:

Float

Default value:

-1

eccentric_RLOF_model

Description:

Chooses which model is used to handle eccentric RLOF. Default is RLOF_ECCENTRIC_AS_CIRCULAR, i.e. ignore the eccentricity. Note: requires force_corotation_of_primary_and_orbit to be FALSE.

Parameter input type:

Integer

Default value:

0

force_circularization_on_RLOF

Description:

If True forces circularization of stars and orbit when RLOF starts, this is as in the BSE algorithm. (True).

Parameter input type:

True|False

Default value:

True

resolve_stellar_type_changes

Description:

If True forces binary_c to resolve stellar type changes to within the minimum timestep. Warning: can be slow. (True).

Parameter input type:

True|False

Default value:

False

nonconservative_angmom_gamma

Description:

Mass lost from the system (but NOT from a stellar wind) takes a fraction gamma of the orbital angular momentum with it. Set to -1 to take the specific angular momentum of the donor star. Set to -2 to take super-Eddington, nova and disk-wind angular momenta as if a wind from the accretor.

Parameter input type:

Float

Default value:

-1

Macros:

[‘RLOF_NONCONSERVATIVE_GAMMA_DONOR = -1’, ‘RLOF_NONCONSERVATIVE_GAMMA_ISOTROPIC = -2’, ‘Random variation : int array or double scalar from [-2,1] or in range 0 to 10’]

max_stellar_angmom_change

Description:

Maxmimum fractional change in stellar angular momentum allowed before a timestep is rejected (0.05).

Parameter input type:

Float

Default value:

0.05

id_number

Description:

Set the integer ID number of the system.

Parameter input type:

Integer

Default value:

0

Extra:

0

decretion_disc_radius_algorithm

Description:

Decretion disc radius calculation algorithm.

Parameter input type:

Integer

Default value:

0

Macros:

[‘DECRETION_DISC_RADIUS_ZHANG_RIMULO = 0’, ‘DECRETION_DISC_RADIUS_COE_KIRK_2015 = 1’]

Extra:

0

system_list

Description:

Set the path to the a system_list file which has a list of arguments on each line.

Parameter input type:

String

Default value:

Extra:

core_collapse_energy

Description:

Set algorithm for computation of the energy of a core-collapse supernova. 0 = CORE_COLLAPSE_ENERGY_ONE_FOE = always 10^51 erg. 1 = CORE_COLLAPSE_ENERGY_MARASSI2019 use data from Marassi et al. (2019).

Parameter input type:

Integer

Default value:

0

Macros:

[‘CORE_COLLAPSE_ENERGY_ONE_FOE = 0’, ‘CORE_COLLAPSE_ENERGY_MARASSI2019 = 1’]

Extra:

0

baryonic_to_gravitational_remnant_mass_prescription

Description:

Baryonic to gravitational remnant mass prescription. During the collapse of a massive star the proto-compact object loses mass through neutrino emission.

Parameter input type:

Integer

Default value:

1

Macros:

[‘BARY_TO_GRAV_FRYER_2012 = 0’, ‘BARY_TO_GRAV_NEUTRINO_EMISSION_LIMITED = 1’]

Extra:

0

Section: nucsyn

nucsyn

Option

Description

third_dup

Description:

If True, enables third dredge up. Requires NUCSYN and NUCSYN_THIRD_DREDGE_UP.

Parameter input type:

True|False

Default value:

True

third_dup_multiplier

Description:

Usage: –third_dup_multiplier <i> <f>. Multiplies the abundance of element <i> by <f> during third dredge up.

Parameter input type:
Default value:

NULL

Extra:

1.0

NeNaMgAl

Description:

Enables NeNaMgAl reaction network. Requires NUCSYN and NUCSYN_HBB.

Parameter input type:

True|False

Default value:

NULL

Extra:

Ignore

nucsyn_network%d

Description:

Usage: –nucsyn_network%d <boolean>. Turn a nuclear network on or off.

Parameter input type:

Boolean(scanf)

Default value:

NULL

nucsyn_network_error%d

Description:

Usage: –nucsyn_network_error%d <f>. Threshold error in nuclear network solver for network %d.

Parameter input type:

Float(scanf)

Default value:

NULL

nucreacmult%d

Description:

Usage: –nucreacmult%d <f>. Multiply nuclear reaction given by the integer %d (integer) by f (float).

Parameter input type:

Float(scanf)

Default value:

NULL

nucsyn_metallicity

Description:

This sets the metallicity of the nucleosynthesis algorithms, i.e. the amount (by mass) of matter which is not hydrogen or helium. Usually you’d just set this with the metallicity parameter, but if you want the nucleosynthesis to be outside the range of the stellar evolution algorithm (e.g. Z=0 or Z=0.04) then you need to use nucsyn_metallicity. That said, it’s also outside the range of some of the nucleosynthesis algorithms as well, so you have been warned!.

Parameter input type:

Float

Default value:

NULL

Macros:

[‘DEFAULT_TO_METALLICITY = -1’, ‘Random variation : double between 0 and 0.05 ‘]

nucsyn_solver

Description:

Choose the solver used in nuclear burning. 0 = KAPS_RENTROP is a Kaps-Rentrop scheme (fast, not great for stiff problems), 1 = LSODA (Adams/BSF switcher), 2 = CVODE library (https://computing.llnl.gov/projects/sundials. Default 0.

Parameter input type:

Unsigned integer

Default value:

NULL

Extra:

0

initial_abundance_mix

Description:

Initial abundance mixture: 0=AG89, 1=Karakas 2002, 2=Lodders 2003, 3=Asplund 2005 (not available?), 4=Garcia Berro, 5=Grevesse Noels 1993.

Parameter input type:

Unsigned integer

Default value:

NULL

Extra:

0

init_abund

Description:

Usage: –init_abund <i> <X>. Sets the initial abundance of isotope number <i> to mass fraction <X>.

Parameter input type:
Default value:

NULL

Extra:

0.02

init_abund_mult

Description:

Usage: –init_abund_mult <i> <f>. Multiplies the initial abundance of isotope number <i> by <f>.

Parameter input type:
Default value:

NULL

Extra:

1.0

init_abund_dex

Description:

Usage: –init_abund_dex <i> <f>. Changes the initial abundance of isotope number <i> by <f> dex.

Parameter input type:
Default value:

NULL

Extra:

0.0

init_abunds_only

Description:

If True, outputs only the initial abundances, then exits.

Parameter input type:

True|False

Default value:

NULL

initial_abunds_only

Description:

If True, outputs only the initial abundances, then exits.

Parameter input type:

True|False

Default value:

NULL

no_thermohaline_mixing

Description:

If True, disables thermohaline mixing.

Parameter input type:

True|False

Default value:

NULL

lithium_GB_post_Heflash

Description:

Sets the lithium abundances after the helium flash. Requires NUCSYN and LITHIUM_TABLES.

Parameter input type:

Float

Default value:

NULL

lithium_GB_post_1DUP

Description:

Sets the lithium abundance after first dredge up. Requires NUCSYN and LITHIUM_TABLES.

Parameter input type:

Float

Default value:

NULL

lithium_hbb_multiplier

Description:

Multiplies the lithium abundances on the AGB during HBB (based on Karakas/Fishlock et al models).Requires NUCSYN and LITHIUM_TABLES.

Parameter input type:

Float

Default value:

NULL

angelou_lithium_decay_function

Description:

Functional form which describes Li7 decay. Requires NUCSYN and NUCSYN_ANGELOU_LITHIUM. Choices are: 0 expoential (see angelou_lithium_decay_time).

Parameter input type:

Integer

Default value:

NULL

angelou_lithium_LMMS_time

Description:

Time at which lithium manufacture is triggered in a low-mass (convective) main sequence (Myr). Requires NUCSYN and NUCSYN_ANGELOU_LITHIUM. Ignored if 0 (for the start, use 1e-6).

Parameter input type:

Float

Default value:

NULL

angelou_lithium_MS_time

Description:

Time at which lithium manufacture is triggered on the main sequence (Myr). Requires NUCSYN and NUCSYN_ANGELOU_LITHIUM. Ignored if 0 (for the start, use 1e-6).

Parameter input type:

Float

Default value:

NULL

angelou_lithium_HG_time

Description:

Time at which lithium manufacture is triggered on the Hertzsprung gap (Myr). Requires NUCSYN and NUCSYN_ANGELOU_LITHIUM. Ignored if 0 (for the start, use 1e-6).

Parameter input type:

Float

Default value:

NULL

angelou_lithium_GB_time

Description:

Time at which lithium manufacture is triggered on the giant branch (Myr). Requires NUCSYN and NUCSYN_ANGELOU_LITHIUM. Ignored if 0 (for the start, use 1e-6).

Parameter input type:

Float

Default value:

NULL

angelou_lithium_CHeB_time

Description:

Time at which lithium manufacture is triggered during core helium burning (Myr). Requires NUCSYN and NUCSYN_ANGELOU_LITHIUM. Ignored if 0 (for the start, use 1e-6).

Parameter input type:

Float

Default value:

NULL

angelou_lithium_EAGB_time

Description:

Time at which lithium manufacture is triggered on the early AGB (Myr). Requires NUCSYN and NUCSYN_ANGELOU_LITHIUM. Ignored if 0 (for the start, use 1e-6).

Parameter input type:

Float

Default value:

NULL

angelou_lithium_TPAGB_time

Description:

Time at which lithium manufacture is triggered on the thermally pulsing AGB (Myr). Requires NUCSYN and NUCSYN_ANGELOU_LITHIUM. Ignored if 0 (for the start, use 1e-6).

Parameter input type:

Float

Default value:

NULL

angelou_lithium_LMMS_decay_time

Description:

Decay time for surface lithium abundance during the low-mass (convective) main sequence (Myr). Requires NUCSYN and NUCSYN_ANGELOU_LITHIUM. Ignored if 0 (for the start, use 1e-6).

Parameter input type:

Float

Default value:

NULL

angelou_lithium_MS_decay_time

Description:

Decay time for surface lithium abundance on the main sequence (Myr). Requires NUCSYN and NUCSYN_ANGELOU_LITHIUM. Ignored if 0 (for the start, use 1e-6).

Parameter input type:

Float

Default value:

NULL

angelou_lithium_HG_decay_time

Description:

Decay time for surface lithium abundance on the Hertzsprung gap (Myr). Requires NUCSYN and NUCSYN_ANGELOU_LITHIUM. Ignored if 0 (for the start, use 1e-6).

Parameter input type:

Float

Default value:

NULL

angelou_lithium_GB_decay_time

Description:

Decay time for surface lithium abundance on the giant branch (Myr). Requires NUCSYN and NUCSYN_ANGELOU_LITHIUM. Ignored if 0 (for the start, use 1e-6).

Parameter input type:

Float

Default value:

NULL

angelou_lithium_CHeB_decay_time

Description:

Decay time for surface lithium abundance during core helium burning (Myr). Requires NUCSYN and NUCSYN_ANGELOU_LITHIUM. Ignored if 0 (for the start, use 1e-6).

Parameter input type:

Float

Default value:

NULL

angelou_lithium_EAGB_decay_time

Description:

Decay time for surface lithium abundance on the early AGB (Myr). Requires NUCSYN and NUCSYN_ANGELOU_LITHIUM. Ignored if 0 (for the start, use 1e-6).

Parameter input type:

Float

Default value:

NULL

angelou_lithium_TPAGB_decay_time

Description:

Decay time for surface lithium abundance on the thermally pulsing AGB (Myr). Requires NUCSYN and NUCSYN_ANGELOU_LITHIUM. Ignored if 0 (for the start, use 1e-6).

Parameter input type:

Float

Default value:

NULL

angelou_lithium_LMMS_massfrac

Description:

Lithium mass fraction when its manufacture is triggered during the low-mass (convective) main sequence (Myr). Requires NUCSYN and NUCSYN_ANGELOU_LITHIUM. Ignored if 0 (for the start, use 1e-6).

Parameter input type:

Float

Default value:

NULL

angelou_lithium_MS_massfrac

Description:

Lithium mass fraction when its manufacture is triggered on the main sequence (Myr). Requires NUCSYN and NUCSYN_ANGELOU_LITHIUM. Ignored if 0 (for the start, use 1e-6).

Parameter input type:

Float

Default value:

NULL

angelou_lithium_HG_massfrac

Description:

Lithium mass fraction when its manufacture is triggered on the Hertzsprung gap (Myr). Requires NUCSYN and NUCSYN_ANGELOU_LITHIUM. Ignored if 0 (for the start, use 1e-6).

Parameter input type:

Float

Default value:

NULL

angelou_lithium_GB_massfrac

Description:

Lithium mass fraction when its manufacture is triggered on the giant branch (Myr). Requires NUCSYN and NUCSYN_ANGELOU_LITHIUM. Ignored if 0 (for the start, use 1e-6).

Parameter input type:

Float

Default value:

NULL

angelou_lithium_CHeB_massfrac

Description:

Lithium mass fraction when its manufacture is triggered during core helium burning (Myr). Requires NUCSYN and NUCSYN_ANGELOU_LITHIUM. Ignored if 0 (for the start, use 1e-6).

Parameter input type:

Float

Default value:

NULL

angelou_lithium_EAGB_massfrac

Description:

Lithium mass fraction when its manufacture is triggered on the early AGB (Myr). Requires NUCSYN and NUCSYN_ANGELOU_LITHIUM. Ignored if 0 (for the start, use 1e-6).

Parameter input type:

Float

Default value:

NULL

angelou_lithium_TPAGB_massfrac

Description:

Lithium mass fraction when its manufacture is triggered on the thermally pulsing AGB (Myr). Requires NUCSYN and NUCSYN_ANGELOU_LITHIUM. Ignored if 0 (for the start, use 1e-6).

Parameter input type:

Float

Default value:

NULL

angelou_lithium_vrot_trigger

Description:

Equatorial rotational velocity at which lithium manufacture is triggered (km/s). Requires NUCSYN and NUCSYN_ANGELOU_LITHIUM. Ignored if 0.

Parameter input type:

Float

Default value:

NULL

angelou_lithium_vrotfrac_trigger

Description:

Fraction of Keplerian (breakup) equatorial rotational velocity at which lithium manufacture is triggered (must be <1, ignored if 0). Requires NUCSYN and NUCSYN_ANGELOU_LITHIUM. Ignored if 0.

Parameter input type:

Float

Default value:

NULL

Section: output

output

Option

Description

cf_amanda_log

Description:

Enable logging to compare to Amanda’s models.

Parameter input type:

True|False

Default value:

NULL

float_overflow_checks

Description:

Turn on to enable floating-point overflow checks at the end of each timestep, if they are available. 0=off, 1=warn (stderr) on failure, 2=exit on failure (0).

Parameter input type:

Integer

Default value:

0

save_pre_events_stardata

Description:

Enable this to save a copy of stardata to stardata->pre_events_stardata just before an event.

Parameter input type:

True|False

Default value:

False

disable_end_logging

Description:

Disable the logging that happens at the end of the evolution.

Parameter input type:

True|False

Default value:

False

ensemble

Description:

Turn on ensemble calculations and output.

Parameter input type:

True|False

Default value:

False

ensemble_filters_off

Description:

Sets all ensemble filters to be off (FALSE) - these can then be enabled one-by-one with –ensemble_filter_[…] TRUE.

Parameter input type:

True|False

Default value:

False

ensemble_filter_%d

Description:

Turn on or off ensemble filter <n>. For a list of filters, see ensemble_macros.h.

Parameter input type:

Boolean(scanf)

Default value:

NULL

ensemble_filter_SCALARS

Description:

Turn on or off ensemble filter SCALARS. For a list of filters, see ensemble_macros.h.

ensemble_filter_HRD

Description:

Turn on or off ensemble filter HRD. For a list of filters, see ensemble_macros.h.

ensemble_filter_SUPERNOVAE

Description:

Turn on or off ensemble filter SUPERNOVAE. For a list of filters, see ensemble_macros.h.

ensemble_filter_CHEMICALLY_PECULIAR

Description:

Turn on or off ensemble filter CHEMICALLY_PECULIAR. For a list of filters, see ensemble_macros.h.

ensemble_filter_SPECTRAL_TYPES

Description:

Turn on or off ensemble filter SPECTRAL_TYPES. For a list of filters, see ensemble_macros.h.

ensemble_filter_INITIAL_DISTRIBUTIONS

Description:

Turn on or off ensemble filter INITIAL_DISTRIBUTIONS. For a list of filters, see ensemble_macros.h.

ensemble_filter_TEST

Description:

Turn on or off ensemble filter TEST. For a list of filters, see ensemble_macros.h.

ensemble_filter_ORBIT

Description:

Turn on or off ensemble filter ORBIT. For a list of filters, see ensemble_macros.h.

ensemble_filter_MASS_FUNCTIONS

Description:

Turn on or off ensemble filter MASS_FUNCTIONS. For a list of filters, see ensemble_macros.h.

ensemble_filter_LUMINOSITY_FUNCTIONS

Description:

Turn on or off ensemble filter LUMINOSITY_FUNCTIONS. For a list of filters, see ensemble_macros.h.

ensemble_filter_VEQ_FUNCTIONS

Description:

Turn on or off ensemble filter VEQ_FUNCTIONS. For a list of filters, see ensemble_macros.h.

ensemble_filter_MERGED

Description:

Turn on or off ensemble filter MERGED. For a list of filters, see ensemble_macros.h.

ensemble_filter_STELLAR_TYPE_COUNTS

Description:

Turn on or off ensemble filter STELLAR_TYPE_COUNTS. For a list of filters, see ensemble_macros.h.

ensemble_filter_CHEMICAL_YIELDS

Description:

Turn on or off ensemble filter CHEMICAL_YIELDS. For a list of filters, see ensemble_macros.h.

ensemble_filter_EMP

Description:

Turn on or off ensemble filter EMP. For a list of filters, see ensemble_macros.h.

ensemble_filter_CBDISCS

Description:

Turn on or off ensemble filter CBDISCS. For a list of filters, see ensemble_macros.h.

ensemble_filter_HRD_PERIOD_DISTRIBUTIONS

Description:

Turn on or off ensemble filter HRD_PERIOD_DISTRIBUTIONS. For a list of filters, see ensemble_macros.h.

ensemble_filter_HRD_TIME_SLICES

Description:

Turn on or off ensemble filter HRD_TIME_SLICES. For a list of filters, see ensemble_macros.h.

ensemble_filter_HRD_PERIOD_DISTRIBUTIONS_TIME_SLICES

Description:

Turn on or off ensemble filter HRD_PERIOD_DISTRIBUTIONS_TIME_SLICES. For a list of filters, see ensemble_macros.h.

ensemble_filter_HRD_MASS_DISTRIBUTIONS

Description:

Turn on or off ensemble filter HRD_MASS_DISTRIBUTIONS. For a list of filters, see ensemble_macros.h.

ensemble_filter_GAIA_CMD

Description:

Turn on or off ensemble filter GAIA_CMD. For a list of filters, see ensemble_macros.h.

ensemble_filter_GAIA_CMD_TIME_SLICES

Description:

Turn on or off ensemble filter GAIA_CMD_TIME_SLICES. For a list of filters, see ensemble_macros.h.

ensemble_filter_GAIA_CMD_PERIOD_DISTRIBUTIONS

Description:

Turn on or off ensemble filter GAIA_CMD_PERIOD_DISTRIBUTIONS. For a list of filters, see ensemble_macros.h.

ensemble_filter_GAIA_CMD_PERIOD_DISTRIBUTIONS_TIME_SLICES

Description:

Turn on or off ensemble filter GAIA_CMD_PERIOD_DISTRIBUTIONS_TIME_SLICES. For a list of filters, see ensemble_macros.h.

ensemble_filter_GAIA_CMD_MASS_DISTRIBUTIONS

Description:

Turn on or off ensemble filter GAIA_CMD_MASS_DISTRIBUTIONS. For a list of filters, see ensemble_macros.h.

ensemble_filter_PULSATORS

Description:

Turn on or off ensemble filter PULSATORS. For a list of filters, see ensemble_macros.h.

ensemble_filter_TIDES

Description:

Turn on or off ensemble filter TIDES. For a list of filters, see ensemble_macros.h.

ensemble_filter_PNE

Description:

Turn on or off ensemble filter PNE. For a list of filters, see ensemble_macros.h.

ensemble_filter_RRLYRAE

Description:

Turn on or off ensemble filter RRLYRAE. For a list of filters, see ensemble_macros.h.

ensemble_filter_TRANSIENTS

Description:

Turn on or off ensemble filter TRANSIENTS. For a list of filters, see ensemble_macros.h.

ensemble_filter_COMENV

Description:

Turn on or off ensemble filter COMENV. For a list of filters, see ensemble_macros.h.

ensemble_filter_HRD_BOKEH

Description:

Turn on or off ensemble filter HRD_BOKEH. For a list of filters, see ensemble_macros.h.

ensemble_filter_CMD_BOKEH

Description:

Turn on or off ensemble filter CMD_BOKEH. For a list of filters, see ensemble_macros.h.

ensemble_filter_IFMR

Description:

Turn on or off ensemble filter IFMR. For a list of filters, see ensemble_macros.h.

ensemble_filter_BeXRB

Description:

Turn on or off ensemble filter BeXRB. For a list of filters, see ensemble_macros.h.

ensemble_filter_NUMBER

Description:

Turn on or off ensemble filter NUMBER. For a list of filters, see ensemble_macros.h.

ensemble_legacy_ensemble

Description:

Turn on ensemble legacy population output.

Parameter input type:

True|False

Default value:

False

legacy_yields

Description:

Turn on ensemble legacy yield output.

Parameter input type:

True|False

Default value:

NULL

ensemble_defer

Description:

Defer ensemble output.

Parameter input type:

True|False

Default value:

False

ensemble_dt

Description:

When doing ensemble calculations, data are stored and/or output every ensemble_dt Myr. See also ensemble, ensemble_logdt, ensemble_startlogtime.

Parameter input type:

Float

Default value:

1

ensemble_logdt

Description:

When doing ensemble calculations, and when logensembletimes is set, the ensemble is stored/output every ensemble_logdt Myr. See also ensemble, ensemble_dt, ensemble_startlogtime.

Parameter input type:

Float

Default value:

0.1

ensemble_startlogtime

Description:

Start log ensemble data storage/calculations/output at ensemble_startlogtime. See also ensemble, ensemble_dt, ensemble_startlogtime.

Parameter input type:

Float

Default value:

0.1

ensemble_logtimes

Description:

When doing ensemble calculations/output, set this to act at log times rather than linear times.

Parameter input type:

True|False

Default value:

False

postagb_legacy_logging

Description:

Turn on post-AGB legacy logging.

Parameter input type:

True|False

Default value:

False

disc_legacy_logging

Description:

Turn on disc legacy logging.

Parameter input type:

True|False

Default value:

False

EMP_logg_maximum

Description:

Maximum logg that EMP stars are allowed to have. See Izzard et al 2009. See also CEMP_cfe_minimum, NEMP_nfe_minimum, EMP_minimum_age.

Parameter input type:

Float

Default value:

NULL

EMP_minimum_age

Description:

Minimum age that EMP stars are required to have. See Izzard et al 2009. See also CEMP_cfe_minimum, NEMP_nfe_minimum, EMP_logg_maximum.

Parameter input type:

Float

Default value:

NULL

EMP_feh_maximum

Description:

Maximum [Fe/H] that an EMP stars may have. See Izzard et al 2009. See also CEMP_cfe_minimum, NEMP_nfe_minimum, EMP_logg_maximum, EMP_minimum_age. Default -2.0.

Parameter input type:

Float

Default value:

NULL

CEMP_cfe_minimum

Description:

Minimum [C/Fe] that CEMP stars are required to have. See Izzard et al 2009. See also NEMP_cfe_minimum, EMP_logg_maximum, EMP_minimum_age. Default 0.7.

Parameter input type:

Float

Default value:

NULL

NEMP_cfe_minimum

Description:

Minimum [N/Fe] that NEMP stars are required to have. See Izzard et al 2009, Pols et al. 2012. See also CEMP_cfe_minimum, EMP_logg_maximum, EMP_minimum_age. Default 1.0.

Parameter input type:

Float

Default value:

NULL

thick_disc_start_age

Description:

Lookback time for the start of the thick disc star formation, e.g. 13e3 Myr. Units = Myr.

Parameter input type:

Float

Default value:

NULL

thick_disc_end_age

Description:

Lookback time for the end of the thick disc star formation, e.g. 4e3 Myr. Units = Myr.

Parameter input type:

Float

Default value:

NULL

thick_disc_logg_min

Description:

Minimum logg for thick disc giants to be logged.

Parameter input type:

Float

Default value:

NULL

thick_disc_logg_max

Description:

Maximum logg for thick disc giants to be logged.

Parameter input type:

Float

Default value:

NULL

escape_velocity

Description:

A parameter used in constructing galactic chemical evolution (GCE) models. If the stellar wind velocity exceeds this value, any chemical yield from the wind is ignored, i.e. assumed lost. (km/s) Requires NUCSYN_GCE_OUTFLOW_CHECKS. Default 1e9 km/s. See also escape_fraction.

Parameter input type:

Float

Default value:

NULL

escape_fraction

Description:

A parameter used in constructing galactic chemical evolution (GCE) models. If the stellar wind velocity exceeds this value, any chemical yield from the wind is ignored, i.e. assumed lost. (km/s) Requires NUCSYN_GCE_OUTFLOW_CHECKS. Default 0.0. See also escape_velocity.

Parameter input type:

Float

Default value:

NULL

colour_log

Description:

If set to True, thelog is coloured with ANSI colour formatting. Requires FILE_LOG to be defined. (If clean_log is True, colours are never used.).

Parameter input type:

True|False

Default value:

False

Extra:

log_legacy_stellar_types

Description:

If set to True, we output the stellar types in the log as numbers rather than (more human-readable) letters.

Parameter input type:

True|False

Default value:

False

Extra:

log_filename

Description:

Location of the output logging filename. If set to “/dev/null” then there is no logging.

Parameter input type:

String

Default value:

/tmp/c_log.dat

Extra:

log_separator

Description:

Separator used in the standard log. By default this is an empty string.

Parameter input type:

String

Default value:

Extra:

log_period_unit

Description:

Define the unit used to log orbital periods. If ‘auto’ or 0, the default, automatically choses the unit.

Parameter input type:

Integer

Default value:

0

Macros:

[‘auto = 0’, ‘tplanck = 1’, ‘s = 2’, ‘day = 3’, ‘d = 4’, ‘hr = 5’, ‘yr = 6’, ‘y = 7’, ‘Kyr = 8’, ‘Ky = 9’, ‘Myr = 10’, ‘My = 11’, ‘Gyr = 12’, ‘Gy = 13’]

clean_log

Description:

If TRUE, all other log formatting is ignored, arrows are off, so that the log is parseable by 3rd-party applications in a fixed form. (FALSE).

Parameter input type:

True|False

Default value:

False

log_arrows

Description:

Add arrows to the output log to show whether values are increasing or decreasing. (If clean_log is True, arrows are never output.).

Parameter input type:

True|False

Default value:

False

Extra:

stopfile

Description:

File which, when it exists, will stop the current binary_c repeat run.

Parameter input type:

String

Default value:

Extra:

stardata_dump_filename

Description:

Location of the stardata dump file.

Parameter input type:

String

Default value:

Extra:

stardata_load_filename

Description:

Location of the stardata file to load.

Parameter input type:

String

Default value:

Extra:

api_log_filename_prefix

Description:

Location of the output logging filename prefix for the API. If set to “/dev/null” then there is no logging.

Parameter input type:

String

Default value:

Extra:

0

hrdiag_output

Description:

Set to True to output high time-resolution Hertzstrpung-Russell diagram information. Requires HRDIAG.

Parameter input type:

True|False

Default value:

NULL

internal_buffering

Description:

Experimental. Set to non-zero values to implement internal buffering prior to output. For use with binary_grid, you shouldn’t really be playing with this.

Parameter input type:

Integer

Default value:

2

Macros:

[‘INTERNAL_BUFFERING_OFF = 0’, ‘INTERNAL_BUFFERING_PRINT = 1’, ‘INTERNAL_BUFFERING_STORE = 2’]

wtts_log

Description:

If True, enables log file output for WTTS2.

Parameter input type:

True|False

Default value:

False

fabian_imf_log_time

Description:

Time at which to output for Fabian Schneider’s IMF project. Requires FABIAN_IMF_LOG.

Parameter input type:

Float

Default value:

NULL

Extra:

Ignore

fabian_imf_log_timestep

Description:

Timestep for Fabian Schneider’s IMF project logging. Requires FABIAN_IMF_LOG.

Parameter input type:

Float

Default value:

NULL

Extra:

Ignore

version

Description:

Display binary_c version and build information. Also performs timing tests.

Parameter input type:
Default value:

NULL

Extra:

Ignore

dumpversion

Description:

Display binary_c version number (short format).

Parameter input type:
Default value:

NULL

Extra:

Ignore

version_only

Description:

Display binary_c version number and build information, but do not perform timing tests or anything that requires stardata to be non-NULL.

Parameter input type:
Default value:

NULL

Extra:

Ignore

tides_diagnosis_log

Description:

Enable logging to test MINT tides. Requires MINT. Choices are: 0 disabled, 1 enable lambda test.

Parameter input type:

Integer

Default value:

0

Extra:

Ignore

YBC_path

Description:

Path to the YBC bolometric correction database (git clone https://gitlab.com/cycyustc/ybc_tables).

Parameter input type:

String

Default value:

Extra:

0

YBC_listfile

Description:

File to use as YBC’s .list. If omitted, this is constructed for you.

Parameter input type:

String

Default value:

Extra:

0

YBC_instruments

Description:

Comma-separated list of the YBC instruments to be used when computing magnitudes, e.g. “GAIA,SLOAN”.

Parameter input type:

String

Default value:

Extra:

0

YBC_all_instruments

Description:

If True, we set all magnitudes associated with all known and defined instruments in the YBC library (see ybc.h for the YBC_INSTRUMENTS_LIST).

Parameter input type:

True|False

Default value:

False

Extra:

0

event_based_logging_%d

Description:

Enable logging of event type <n> (e.g. SN, RLOF, DCO). The event logstrings will first be stored in the stardata, and at the end of the evolution they will all be printed. This deals with the evol-splitting better.

Parameter input type:

Boolean(scanf)

Default value:

NULL

Extra:

Ignore

event_based_logging_SN

Description:

Enable logging of event type SN (e.g. SN, RLOF, DCO). The event logstrings will first be stored in the stardata, and at the end of the evolution they will all be printed. This deals with the evol-splitting better.

Extra:

Ignore

event_based_logging_RLOF

Description:

Enable logging of event type RLOF (e.g. SN, RLOF, DCO). The event logstrings will first be stored in the stardata, and at the end of the evolution they will all be printed. This deals with the evol-splitting better.

Extra:

Ignore

event_based_logging_DCO

Description:

Enable logging of event type DCO (e.g. SN, RLOF, DCO). The event logstrings will first be stored in the stardata, and at the end of the evolution they will all be printed. This deals with the evol-splitting better.

Extra:

Ignore

Section: input

input

Option

Description

MINT_dir

Description:

Location of MINT algorithm data.

Parameter input type:

String

Default value:

Extra:

MINT_load_state_file

Description:

File from which the MINT state should be loaded.

Parameter input type:

String

Default value:

Extra:

MINT_save_state_file

Description:

File to which the MINT state should be saved.

Parameter input type:

String

Default value:

Extra:

MINT_data_cleanup

Description:

Activate checks on incoming data to try to account for problems. Will make data-loading slower, but may fix a few things.

Parameter input type:

True|False

Default value:

False

Extra:

MINT_MS_rejuvenation

Description:

Turn on or off (hydrogen) main-sequence rejuvenation.

Parameter input type:

True|False

Default value:

True

Extra:

MINT_remesh

Description:

Turn on or off MINT’s remeshing.

Parameter input type:

True|False

Default value:

True

Extra:

MINT_filename_vb

Description:

Turn on or off verbose logging of MINT filename searches. Handy for debugging when you can’t get your data file to load.

Parameter input type:

True|False

Default value:

False

Extra:

MINT_use_ZAMS_profiles

Description:

Use chemical profiles at the ZAMS if MINT_use_ZAMS_profiles is TRUE, otherwise set homogeneous abundances. (Default is TRUE, so we use the profiles if they are available.).

Parameter input type:

True|False

Default value:

True

Extra:

MINT_fallback_to_test_data

Description:

If TRUE, use the MINT test_data directory as a fallback when data is unavailable. (FALSE).

Parameter input type:

True|False

Default value:

False

Extra:

MINT_use_fallback_comenv

Description:

If TRUE, use the BSE common-envelope prescription as a fallback when MINT is not available.

Parameter input type:

True|False

Default value:

False

Extra:

MINT_disable_grid_load_warnings

Description:

Use this to disable MINT’s warnings when loading a grid with, e.g., missing or too much data.

Parameter input type:

True|False

Default value:

False

Extra:

MINT_disable_warnings

Description:

Use this to disable all MINT’s warnings.

Parameter input type:

True|False

Default value:

False

Extra:

MINT_Kippenhahn

Description:

Turn on or off MINT’s Kippenhahn diagrams. If 0, off, if 1, output star 1 (index 0), if 2 output star 2 (index 1). Default 0.

Parameter input type:

Integer

Default value:

0

Extra:

MINT_nshells

Description:

Set the initial number of shells MINT uses in each star when doing nuclear burning. Note: remeshing can change this. If MINT_nshells is 0, shellular burning and other routines that require shells will not be available. (200).

Parameter input type:

Integer

Default value:

200

Extra:

MINT_maximum_nshells

Description:

Set the maximum number of shells MINT uses in each star when doing nuclear burning. Note that this will be limited to MINT_HARD_MAX_NSHELLS. (1000).

Parameter input type:

Integer

Default value:

1000

Extra:

MINT_minimum_nshells

Description:

Set the minimum number of shells MINT uses in each star when doing nuclear burning. Note that this will be greater than or equal to MINT_HARD_MIN_NSHELLS, which is 0 by default. (0).

Parameter input type:

Integer

Default value:

10

Extra:

MINT_Kippenhahn_stellar_type

Description:

Stellar type selector for Kippenhahn plots. Set to -1 to ignore, otherwise the stellar type number for which Kippenhahn plot data should be output.

Parameter input type:

Integer

Default value:

-1

Macros:

[‘LOW_MASS_MS = 0’, ‘MS = 1’, ‘HG = 2’, ‘GIANT_BRANCH = 3’, ‘CHeB = 4’, ‘EAGB = 5’, ‘TPAGB = 6’, ‘HeMS = 7’, ‘HeHG = 8’, ‘HeGB = 9’, ‘HeWD = 10’, ‘COWD = 11’, ‘ONeWD = 12’, ‘NS = 13’, ‘BH = 14’, ‘MASSLESS_REMNANT = 15’]

Extra:

MINT_Kippenhahn_companion_stellar_type

Description:

Companion stellar type selector for Kippenhahn plots. Set to -1 to ignore, otherwise the stellar type number for the companion for which Kippenhahn plot data should be output.

Parameter input type:

Integer

Default value:

-1

Macros:

[‘LOW_MASS_MS = 0’, ‘MS = 1’, ‘HG = 2’, ‘GIANT_BRANCH = 3’, ‘CHeB = 4’, ‘EAGB = 5’, ‘TPAGB = 6’, ‘HeMS = 7’, ‘HeHG = 8’, ‘HeGB = 9’, ‘HeWD = 10’, ‘COWD = 11’, ‘ONeWD = 12’, ‘NS = 13’, ‘BH = 14’, ‘MASSLESS_REMNANT = 15’]

Extra:

MINT_nuclear_burning

Description:

Turn on or off MINT’s nuclear burning algorithm.

Parameter input type:

True|False

Default value:

False

Extra:

MINT_minimum_shell_mass

Description:

Minimum shell mass in MINT’s nuclear burning routines.

Parameter input type:

Float

Default value:

1e-06

Extra:

MINT_maximum_shell_mass

Description:

Maximum shell mass in MINT’s nuclear burning routines. :.

Parameter input type:

Float

Default value:

0.1

Extra:

Section: i/o

i/o

Option

Description

go

Description:

Batchmode control command.

Extra:

Ignore

gogo

Description:

Batchmode control command.

Extra:

Ignore

reset_stars

Description:

Reset the star structures. Used in batchmode.

Extra:

Ignore

reset_stars_defaults

Description:

Reset the star structures and set defaults. Used in batchmode.

Extra:

Ignore

defaults

Description:

Reset all defaults. Used in batchmode.

Extra:

Ignore

echo

Description:

Activate batchmode command echoing, i.e. when you enter a command, binary_c repeats the command then executes it.

Extra:

Ignore

noecho

Description:

Deactivate batchmode command echoing. See ‘echo’.

Extra:

Ignore

noechonow

Description:

Deactivate batchmode command echoing. See ‘echo’.

Extra:

Ignore

bye

Description:

Quit binary_c. Used in batchmode.

Extra:

Ignore

fin

Description:

Batchmode control command.

Extra:

Ignore

reset_prefs

Description:

Reset preferences struct. Used in batchmode.

Extra:

Ignore

status

Description:

Output batchmode status information.

Extra:

Ignore

Section: algorithms

algorithms

Option

Description

repeat

Description:

If > 1, repeats the system as many times as required. Handy if you’re using Monte-Carlo kicks and want to sample the parameter space well. Also, if you are running speed tests this is good to give a statistically more reasonable result. (See e.g. ‘tbse pgo’).

Parameter input type:

Integer

Default value:

1

pause_after_repeat

Description:

If true, pauses after a repeated system is run to wait for someone to press RETURN on stdin.

Parameter input type:

True|False

Default value:

False

random_systems

Description:

Experimental. Use this to apply random initial system parameters (masses, separations, etc.). Useful for testing only.

Parameter input type:

Integer

Default value:

0

Section: misc

misc

Option

Description

defaults_set

Description:

Choose the set of defaults to be used by binary_c.

Parameter input type:
Default value:

NULL

Macros:

[‘DEFAULTS_SET_2022a = 1’]

skip_bad_args

Description:

If not set to SKIP_BAD_ARGS_NONE, skip unmatched arguments assuming they are in the form “–x y”. Can be SKIP_BAD_ARGS_WITH_WARNING (will warn to the output buffer) or SKIP_BAD_ARGS_WITH_NO_WARNING (stays quiet). (SKIP_BAD_ARGS_NONE).

Parameter input type:

Integer

Default value:

0

Macros:

[‘SKIP_BAD_ARGS_NONE = 0’, ‘SKIP_BAD_ARGS_WITH_WARNING = 1’, ‘SKIP_BAD_ARGS_NO_WARNING = 2’]

random_seed

Description:

Random number seed, usually a (possibly negative) integer. Useful for exactly reproducing the evolution of a system which involves a kick (which is a Monte-Carlo, i.e. pseudorandom, process).

Parameter input type:

Integer

Default value:

0

random_systems_seed

Description:

Random number seed for the generation of random systems. See random_systems and random_seed.

Parameter input type:

Integer

Default value:

0

random_skip

Description:

Skip the first <random_seed> random numbers that are generated. Usually this is 0 so they are all used.

Parameter input type:

Integer

Default value:

0

idum

Description:

[NB: deprecated, please use ‘random_seed’ instead.] Random number seed, usually a (possibly negative) integer. Useful for exactly reproducing the evolution of a system which involves a kick (which is a Monte-Carlo, i.e. pseudorandom, process).

Parameter input type:

Integer

Default value:

0

reverse_time

Description:

Make time go backwards. To be considered very experimental!.

Parameter input type:

True|False

Default value:

NULL

start_time

Description:

Start time for the simulation.

Parameter input type:

Float

Default value:

0

warmup_cpu

Description:

Uses the CPU at maximum power the given number of seconds, prior to running normal stellar evolution.

Parameter input type:
Default value:

NULL

Extra:

Ignore

help

Description:

Display help pages. Usage: –help <help topic>.

Parameter input type:
Default value:

NULL

Extra:

Ignore

argopts

Description:

Display argument options. Usage: –argopts <argument>.

Parameter input type:
Default value:

NULL

Extra:

Ignore

help_all

Description:

Display all help pages.

Parameter input type:
Default value:

NULL

Extra:

Ignore

list_args

Description:

Display list of arguments with their default values. Useful for batchmode.

Parameter input type:
Default value:

NULL

Extra:

Ignore

bjorn

Description:

Usage: –bjorn … shows an ASCII-art picture of Björn (requires jp2a to be installed).

Parameter input type:
Default value:

NULL

Extra:

0

logo

Description:

Usage: –logo … shows an ASCII-art picture of the binary_c logo (requires jp2a to be installed).

Parameter input type:
Default value:

NULL

Extra:

0

nanchecks

Description:

Turn nanchecks on or off. Requires SOFT_NANCHECK to be defined. True by default if so, otherwise ignored.

Parameter input type:

True|False

Default value:

False

random_system_list

Description:

Output a list of nrepeat random systems (see also “repeat”) and exit.

Parameter input type:
Default value:

NULL